newton.c

C语言程 反对萨芬记录卡上的交付离开的撒娇发了啊

/*牛顿迭代法*/
#include<stdio.h>
#include<conio.h>
#include<math.h>
#define EPS 1e-6
#define ETA 1e-8
#define N 100
void main()
{ float f(float);
  float f1(float);
  float x0,y0;
   float Newton(float(*)(float),float(*)(float),float) ;
   clrscr();
   printf("Please input x0\n");
   scanf("%f",&x0);
   printf("x(0)=%f\n",x0);
   y0=Newton(f,f1,x0);
   printf("\nThe root of the equation is x=%f\n",y0);
   getch();
 }
float Newton(float(*f)(float),float(*f1)(float),float x0)
  { float x1,d;
     int k=0;
     do { x1=x0-f(x0)/f1(x0);
          if(k++>N||(fabs(f1(x1))<EPS))
             { printf("\nNewton method failed!");
              getch();
              exit();
              }
      d=(fabs(x1)<1?(x1-x0):(x1-x0)/x1);
      x0=x1;
      printf("x(%d)=%f\t",k,x0);
      }
     while(fabs(d)>EPS&&fabs(f(x1))>ETA) ;
     return x1;
     }
  float f(float x)
     { return 3*x+sin(x)-exp(x);
     }
  float f1(float x)
     {  return 3+cos(x)-exp(x);
     }