您的当前位置:首页正文

牛顿迭代法c++程序

来源:个人技术集锦


#include

#include

#include

#define EPS 1e-6

#define DELTA 1e-6

#define N 100

float f(float);

float fd(float);

float f(float x)

{

return x*x*x+x*x-3*x-3;

}

float fd(float x)

{

return 3*x*x+2*x-3;

}

main()

{

int i=0,n=0,c;

float x[N];

printf(\"ÇëÊäÈëx0\\n\");

scanf(\"%f\

do{

x[i+1]=x[i]-f(x[i])/fd(x[i]);

if(fd((n++>N)||x[i])==0)

{

printf(\"Ë㷨ʧ°Ü\");

break;

}

x[i]=x[i+1];

i++;

}

while(fabs(x[i]-x[i-1])>EPS||fabs(f(x[i]))>DELTA);

for(i=0;i<=n;i++)

{printf(\"x[%d]=%f\\

}

}

因篇幅问题不能全部显示,请点此查看更多更全内容