#include #include /* newton.c : 1999 2 28 */ double newton(double a, int n); /* main */ main() { double x; x=newton(2.0,10); printf(" %f\n",x); } /* newton method */ double newton(double a, int n) { double x,x0; int i; x0=a/2.0; for(i=1; i<=n; i++){ x=(x0+a/x0)/2.0; if (fabs(x-x0)<=0.0001) break; x0=x; } return x; }