求数值分析打靶法的C语言程序！高悬赏！

求数值分析打靶法的C语言程序！高悬赏！

#include
#include
#include

int main()
{
int i,cc;
double h,epsilon,ma,mb,mc,mm,yfa,yfb,xx,y0,yf,panduan;
double yy[11];
h=0.1;
cc=0;
epsilon=0.0002;
ma=5;mm=ma;
y0=1;yf=exp(1);
yy[0]=y0;
do
{
yy[1]=h*mm+yy[0];
for(i=2;i<=10;i++)
{
xx=i*h;
yy[i]=2*yy[i-1]-yy[i-2]+pow(h, 2)*((2+9*xx)*yy[i-1]-(1+9*xx)*exp(xx));
}
panduan=(yy[10]-yf)*(yy[10]-yf)-epsilon*epsilon;
if(panduan>0&&cc==0)
{
yfa=yy[10];
mb=yf*ma/yfa;
mm=mb;
}
else
{
yfb=yy[10];
mc=ma+(mb-ma)*(yf-yfa)/(yfb-yfa);
mm=mc;
ma=mb;yfa=yfb;
}
cc=cc+1;
}
while(panduan>0);

printf("%f\n",yy[0]);
printf("%f\n",yy[1]);
printf("%f\n",yy[2]);
printf("%f\n",yy[3]);
printf("%f\n",yy[4]);
printf("%f\n",yy[5]);
printf("%f\n",yy[6]);
printf("%f\n",yy[7]);
printf("%f\n",yy[8]);
printf("%f\n",yy[9]);
printf("%f\n",yy[10]);
printf("%d\n",cc);
printf("%f\n",panduan);
system("PAUSE");
return 0;
}
最简单的一个打靶法程序，计算y''=-(2+9x)y=-(1+9x)e^x,y(0)=1,y(1)=e

这么点分。懒的做

确实如此，这个该给RMB了
如果只有线性的，你给200分，我可以考虑一试
把题目写详细点

如以上回答内容为低俗、色情、不良、暴力、侵权、涉及违法等信息，可以点下面链接进行举报！

点此我要举报以上问答信息