www.pudn.com > fxxfcz.rar > DEMO0302.C

#include  
#include  
#include  
#define MN 51 
double f(double x){ return(x*x-5.0);} 
double g(double x){ return(2.0*x);} 
static   int     RN=1; 
static   double  X[MN],Y[MN],Z[MN],T[MN],E[MN]; 
static   double  A0=2.00,B0=3.0; 
static   double  EPS=1.0e-15; 
double solution; 
Operation() 
{  int K; 
   double x; 
   if(f(A0)>0) X[0]=A0; 
       else    X[0]=B0; 
   Y[0]=f(X[0]); 
   Z[0]=g(X[0]); 
   T[0]=X[0]-Y[0]/Z[0]; 
   E[0]=fabs(T[0]-X[0]); 
   for(K=1;K0) x=a; 
     else    x=b; 
  for(k=0;k<100;k++) 
  { y =f(x); 
    z=g(x); 
    t= x-y/z; 
    if(fabs(x-t)