#include <stdio.h>
#include <stdlib.h>
#include <stdint.h>
#include <vector>
#include "aritmetica.hpp"
#include "util.hpp"
using namespace std;
int prec;
numero b_i(int j)
{
numero ret(prec, 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]);
// Vamos a hacer una lista con los terminos para ir sumando del ultimo
// al primero
vector<numero> terminos;
/* Ponemos el primer termino de la seie */
numero a_i(prec, 1);
terminos.push_back(a_i);
/* Loop de la serie, nterm - 1 términos (el 1 lo acabamos de pushear
* arriba) */
for (int i = 0; i < nterm - 1; i++) {
a_i = a_i * b_i(i);
terminos.push_back(a_i);
/* 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 generando los parciales, sumando siempre de atras para
* adelante */
for (int i = 0; i < nterm; i++) {
numero n(prec, 0);
for (int j = 0; j <= i; j++) {
n = n + terminos[i - j];
}
print_result(i + 1, n);
}
return 0;
}