#include <stdio.h>
#include <stdlib.h>
#include <stdint.h>
#include "aritmetica.hpp"
#include "util.hpp"
int prec;
numero b_i(int j)
{
numero ret(prec, 1);
// (1/2 -j) / (j + 1) = [(1 - 2J) / (j+1)] / 2
numero num(prec, 2);
num = num * j; //2*J
num = 1 - num ; //1-2J
numero den(prec, 1);
den = den +j;
ret = ret * (num / den) / 2;
//printf(" %d %.60f 0x%" print_64t "x\n", j, ret.to_double(),
// to_u64(ret));
return ret;
}
int main(int argc, char *argv[])
{
if (argc != 3) {
printf("Uso: %s <prec> <nterm>\n", argv[0]);
return 1;
}
prec = atoi(argv[1]);
int nterm = atoi(argv[2]);
/* Valor inicial, el del primer término */
numero n(prec, 1);
print_result(1, n);
numero a_i(prec, 1);
/* Loop de la serie, nterm - 1 términos (el primero, por el cual
* arrancamos, ya esta fijo arriba) */
for (int i = 0; i < nterm - 1; i++) {
numero antes = n;
a_i = a_i * b_i(i);
//printf("%.60f 0x%" print_64t "x\n", a_i.to_double(), to_u64(a_i));
n = n + a_i;
if (n == antes)
{
//printf("Sale por no sumar\n");
break;
}
/* XXX DEBUG - SACAR */
#if 0
printf(" %.60f 0x%" print_64t "x\n",
numerador.to_double(), to_u64(numerador));
printf(" %.60f 0x%" print_64t "x\n",
denominador.to_double(), to_u64(denominador));
printf(" ---\n");
#endif
//printf(" %d %.60f 0x%" print_64t "x\n", i, n.to_double(), to_u64(n));
/* Vamos emitiendo los resultados a medida que los tenemos */
print_result(i + 2, n);
}
return 0;
}