Ejemplo013.c

Go to the documentation of this file.

 
#include <stdio.h>
 
int main(int argc, char *argv[])
{
    (void)argc;
    (void)argv;
 
    int n, i;
    float x;        /* valor para ln(x) */
    float lnx;      /* acumulador de la serie */
    float fct;      /* f = (x-1)/(x+1) */
    float fct2;     /* f^2 */
    float suma;     /* término actual t_i = f^{2i+1}/(2i+1) */
 
    /* Leer n >= 1 */
    do {
        printf("Ingres el numero de terminos: ");
        scanf("%d", &n);
    } while (n < 1);
 
    /* Leer x > 0 */
    do {
        printf("Ingrese el valor de x: ");
        scanf("%f", &x);
    } while (x <= 0);
 
    /*
      Inicialización:
        f = (x-1)/(x+1)
        f2 = f^2
        t0 = f^(1)/1 = f
        lnx = 0
    */
    for (i = 0, fct = (x - 1) / (x + 1), lnx = 0, fct2 = fct * fct, suma = fct;
         i < n;
         i++)
    {
        /* Acumular término t_i */
        lnx += suma;
 
        /*
          Actualizar al siguiente término:
            t_{i+1} = t_i * f^2 * (2i+1)/(2i+3)
          (evita recalcular potencias desde cero)
        */
        suma *= (fct2 * (2*i + 1)) / (2*i + 3);
    }
 
    /* ln(x) = 2 * serie */
    lnx *= 2;
 
    printf("ln(%f) = %f\n", x, lnx);
    return 0;
}