newton法、美しいアルゴリズム

newton.c

#include<stdio.h>
#include<stdlib.h>

double newton(double c)
{
    double x;
    int i;

    x = ( c + c / c ) / 2;

    for(i=0;i < 1000000; i++)
    {
        x = ( x + c /x ) / 2;

    }

    return x;

}

int
main(int argc, char **argv)
{

    double c =2*2*2*2*2*2*2*2*2*2*2*2*2*2*2*2*2*2*2-1,
           t;
    t = newton(c);
    printf("%10.100lf\n",t);



}