unsigned long int i, j;
 
float sumatoria, buff_sumatoria, buff_ciclo_superior, buff_ciclo_inferior, primera_parte, buff_ciclo_inferior2, factorial_superior, factorial_inferior;
 
sumatoria=1103;
buff_sumatoria=0;
buff_ciclo_superior=0;
buff_ciclo_inferior=0;
primera_parte=0;
buff_ciclo_inferior2=0;
factorial_superior=1;
factorial_inferior=1;
 
//sqrt(2)
primera_parte=sqrt(2);
 
//sqrt(2) * 2
primera_parte*=2;
 
//sqrt(2)*2 / 9801
primera_parte/=9801
 
for (i = 1; i <= digitos; i++) {//sabiendo que ramanujan produce 8 digitos con un margen de error de 1% por ciclo, según mis pruebas
        for (j = (i - 1)*4 + 1; j <= i * 4; j++) {
            factorial_superior*=*j;//ciclo acumulador del factorial (4k)!
        }
        //29690 * k
        buff_ciclo_superior= 26390;
        buff_ciclo_superior *= i;
 
        //(1103 + 29690*k)
        buff_ciclo_superior+=1103;
 
        //(4k)! * (1103+29690*k)
        buff_ciclo_superior*=factorial_superior;
 
        //(k!)^4
        factorial_inferior*=i;
        buff_ciclo_inferior = pow(factorial_inferior, 4);
 
        //396^4k
        buff_ciclo_inferior2 = 396;
        buff_ciclo_inferior2 = pow(buff_ciclo_inferior2, 4 * i);
 
        //(k!)^4 * 396^4
        buff_ciclo_inferior*=buff_ciclo_inferior2;
 
        // (4k)!*(1103+29690*k) / (k!)^4*396^4
        buff_sumatoria = buff_ciclo_superior / buff_ciclo_inferior;
 
        sumatoria+=buff_sumatoria;//la sumatoria...
 
    }
 
    //sqrt(2)*2/9801 * SUMATORIA (4k)!*(1103+29690*k) / (k!)^4*396^4
    sumatoria*=primera_parte;
 
    // el resultado original es 1/pi, entonces con esto hago la inversión
    sumatoria = 1 / sumatoria;
 
    //convierte a texto para poder imprimirlo
    return mpf_get_str(NULL, &exponente, 10, 0, sumatoria);