[顯示] [隱藏]

#include <cstdlib>
#include <cstdio>
#include <iostream>
#include <cmath>
/*
利用對分法求非線性方程式
f(x)=x^6-5x^5+x^3-7x^2+7x-20=0
*/
using namespace std;
  int dhrt(double a,double b,double h,double eps,double *x,int m,double (*f)(double))
//  int m;
//  double a,b,h,eps,x[],(*f)();
  { int n,js;
    double z,y,z1,y1,z0,y0;
    n=0; z=a; y=(*f)(z);
    while ((z<=b+h/2.0)&&(n!=m))
      { if (fabs(y)<eps)
          { n=n+1; x[n-1]=z;
            z=z+h/2.0; y=(*f)(z);
          }
        else
          { z1=z+h; y1=(*f)(z1);
            if (fabs(y1)<eps)
              { n=n+1; x[n-1]=z1;
                z=z1+h/2.0; y=(*f)(z);
              }
            else if (y*y1>0.0)
              { y=y1; z=z1;}
            else
              { js=0;
                while (js==0)
                  { if (fabs(z1-z)<eps)
                      { n=n+1; x[n-1]=(z1+z)/2.0;
                        z=z1+h/2.0; y=(*f)(z);
                        js=1;
                      }
                    else
                      { z0=(z1+z)/2.0; y0=(*f)(z0);
                        if (fabs(y0)<eps)
                          { x[n]=z0; n=n+1; js=1;
                            z=z0+h/2.0; y=(*f)(z);
                          }
                        else if ((y*y0)<0.0)
                          { z1=z0; y1=y0;}
                        else { z=z0; y=y0;}
                      }
                  }
              }
          }
      }
    return(n);
  }
  double dhrtf(double x)
//  double x;
  { double z;
    z=(((((x-5.0)*x+3.0)*x+1.0)*x-7.0)*x+7.0)*x-20.0;
    return(z);
  }
int main(int argc, char *argv[])
{
    int i,n;
    static int m=6;
    static double x[6];
    double dhrtf(double);
    n=dhrt(-2.0,5.0,0.2,0.000001,x,m,dhrtf);
    printf("M=%d\n",n);
    for (i=0; i<=n-1; i++)
      printf("x(%d)=%13.7e\n",i,x);
    printf("\n");
    system("PAUSE");
    return EXIT_SUCCESS;
}