原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
#includestdio.h
#includemath.h
#define f(x) x+20/x
#define M 0.618
void main()
{
double y1,y2,x1,x2,x,a,b,e;
int n;
n=1;
printf(请输入收敛精度e=);
scanf(%lf,e);
printf(请输入区间左值a=);
scanf(%lf,a);
printf(请输入区间右值b=);
scanf(%lf,b);
printf(n a b x1 x2 y1 y2 b-a\n);
x1=b-M*(b-a);
x2=a+M*(b-a);
y1=f(x1);
y2=f(x2);
printf(%d %.4lf %.4lf %.4lf %.4lf %.4lf %.4lf\n,n,a,b,x1,x2,y1,y2,b-a);
n=n++;
do
{
if(y1y2)
{
b=x2;
x2=x1;
y2=y1;
x1=b-M*(b-a);
y1=f(x1);
printf(%d %.4lf %.4lf %.4lf %.4lf %.4lf %.4lf %.4lf\n,n,a,b,x1,x2,y1,y2,b-a);
n=n++;
}
else
{
a=x1;
x1=x2;
y1=y2;
x2=a+M*(b-a);
y2=f(x2);
printf(%d %.4lf %.4lf %.4lf %.4lf %.4lf %.4lf %.4lf\n,n,a,b,x1,x2,y1,y2,b-a);
n=n++;
}
}
while(fabs((b-a))=e);
x=(a+b)*0.5;
printf(x=%.5lf f(x)=%.5lf\n,x,f(x));
getchar();
}
#includestdio.h
#includestdlib.h
#includemath.h
#includeconio.h
#define TT 0.08
#define FF 0.0001
#define AC #define AD #define N 2
float ia;
float fny(float *x)
{
float f;
f=(float)(10*pow((x[0]+x[1]-5),2)+pow((x[0]-x[1]),2));
return(f);
}
float *iterate(float *x,float a, float *s)
{
float *x1;
x1=(float *)malloc(N*sizeof(float));
for(int i=0; iN; i++)
x1[i]=x[i]+a*s[i];
return(x1);
}
float func(float *x, float a, float*s)
{
float *x1;
x1=iterate(x,a,s);
float f=fny(x1);
return(f);
}
void finding(float a[3], float f[3], float *xk, float *s)
{
float t=(float)TT;
float a1,f1;
a[0]=0; f[0]=func(xk,a[0],s);
for(int i=0;;i++)
{
a[1]=a[0]+t; f[1]=func(xk,a[1],s);
if(f[1]f[0]) break;
if(fabs(f[1]-f[0])=AD)
{
t=-t; a[0]=a[1]; f[0]=f[1];
}
else
{
if(ia==1) return;
t=t/2; ia=1;
}
}
for(i=0;;i++)
{
a[2]=a[1]+t; f[2]=func(xk,a[2],s);
if(f[2]f[1])break;
t=2*t;
a[0]=a[1]; f[0]=f[1];
a[1]=a[2]; f[1]=f[2];
}
if(a[0]a[2])
{
a1=a[0]; f1=f[0]; a[0]=a[2]; f[0]=f[2]; a[2]=a1; f[2]=f1;
}
return;
}
float lagrange(float *xk,float*ft,float *s)
{
float a[3],f[3];
float b,c,d,aa;
finding(a,f,xk,s);
for(int i
文档评论（0）