最速下降法.txt

该压缩文件有一些经典插值c语言版的源程序

#include<stdio.h>
#include<math.h>

main()
{
float x1,x2,T1,T2,f,g,M=0.001,t,s,p;
clrscr();
x1=1;
x2=1;
T1=2*x1;
T2=8*x2;
t=(T1*T1+T2*T2)/(2*T1*T1+8*T2*T2);
x1=x1-t*T1;
x2=x2-t*T2;
T1=2*x1;
T2=8*x2;
s=T1*T1+T2*T2;
p=sqrt(s);
while(p>M)
{
T1=2*x1;
T2=8*x2;
t=(T1*T1+T2*T2)/(2*T1*T1+8*T2*T2);
x1=x1-t*T1;
x2=x2-t*T2;
T1=2*x1;
T2=8*x2;
s=T1*T1+T2*T2;
p=sqrt(s);
printf("x1=%f\t x2=%f\t  s=%f\n",x1,x2,p);
}
f=x1*x1+4*x2*x2;
printf("the best result is:\n");
printf("x1=%f\t x2=%f\t  f=%f",x1,x2,f);
}