求非线性方程一个实根的牛顿法

#include "stdio.h"

#include "math.h"

#include "equation.h"

/*牛顿迭代法：求非线性方程一个实根*/

int newton_iterative(ITEP np,NLFP fp)

{

   int l;

  double eps,y[2],d,p,x0,x1;

  eps=np->eps;                            /*控制精度*/

  x0=np->x;                                  /*np->x为迭代初值*/

   nlf2(fp,x0,y);

  l=np->iter; d=eps+1.0;             /*牛顿迭代法*/

while((d>=eps)&&(l>0)){

   if(fabs(y[1])<0.0000001)  {

          printf("error\n");  return (-1);

   }

     x1=x0-y[0]/y[1];

   nlf2(fp,x1,y);

  d=fabs(x1-x0);

  p=fabs(y[0]);

   if(p>d)  d=p;

  x0=x1;

  l--;

}

  np->x=x1;   /*np->x返回时为迭代终值*/

  return(np->iter-1);

}

main() {

  double a[]={-1.0,0.0,-2.0,1.0};    /*多项式系数*/

  NLF   fa={3,a};                           /*多项式描述结构*/

 ITE   nit ={50,0.000001,2.5};

int k;

  if((k=newton_interative(&nit,&fa))>=0)

       printf("\nk=%d  x=%13.7e\n",k,nit.x);         /*输出迭代次数和根*/

   printf("\n");

getch();

} /*程序完毕*/

本文出自 “王小毛-爱你！” 博客，转载请与作者联系！