请教牛顿迭代法
给你一个模板吧
#include "math。h"
main()
{double x,x1,f,f1;
double xx1(),yy(),yy1(),pd();
x=0;x1=500;
for(;pd(x,x1)>=1e-5;)
{x=x1;
x1=xx1(x);
}
printf("%f",x1);
}
double yy(double x)
{double f;
f=59。 7622-8013。69/x-5。08*log(x)-log(1。013e5);
return(f);
}
double yy1(double x)
{double f1;
f1=8013。 69/(x*x)...全部
给你一个模板吧
#include "math。h"
main()
{double x,x1,f,f1;
double xx1(),yy(),yy1(),pd();
x=0;x1=500;
for(;pd(x,x1)>=1e-5;)
{x=x1;
x1=xx1(x);
}
printf("%f",x1);
}
double yy(double x)
{double f;
f=59。
7622-8013。69/x-5。08*log(x)-log(1。013e5);
return(f);
}
double yy1(double x)
{double f1;
f1=8013。
69/(x*x)-5。
08/x; return(f1);
}
double xx1(double x)
{double x1;
x1=x-yy(x) 1(x); return(x1);
}
double pd(double x,double x1)
{double y;
y=abs(x-x1);
return(y);
}。收起