《歐拉及改進的歐拉法求解常微分方程》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫。

1、生物信息技術(shù)0801徐聰U200812594#include#includevoidfl(double*y,double*x,double*yy){y[o]=2.o；x[0]=0.0;yy[0]=2.0;for(inti=l;i<=9;i++){x[i]=x[i-l]+0.2;y[i]=y[i-l]+0.2*(y[i-l]-x[i-l]);yy[i]=x[i]+l+exp(x[i]);printf(”若x=%f,計算值是%f，真實值是%f，截斷誤差是%fx[i],y[i],yy[i],y[i]-yy[i]);}}；void

2、f2(double*y,double*x,double*yy){y[o]=i.o；x[0]=0.0;yy[0]=1.0;for(inti=l;i<=9;i++){x[i]=x[i-l]+0.2;y[i]=y[i-l]+0.2*(2*y[i-l]+x[i-l]*x[i-l]);yy[i]=-0.5*(x[i]*x[i]+x[i]+0.5)+1.25*exp(2*x[i]);printf(”若x=%f,計算值是%f，真實值是%f，截斷誤差是%fx[i],y[i],yy[i],y[i]-yy[i]);}}；voidf3(double*y,double*x,dou

3、ble*yy,double*y0){y[o]=2.o；x[0]=0.0;yy[0]=2.0;for(inti=l;i<=9;i++){x[i]=x[i-l]+0.2;y0[i]=y[i-l]+0.2*(y[i-l]-x[i-l]);y[i]=y[i-l]+0.1*(y[i-l]-x[i-l]+y0[i-l]-x[i-l]);yy[i]=x[i]+l+exp(x[j]);printf('喏x=%f,計算值是%f，真實值是%f，截斷誤差是%fn,x[i]zy[i],yy[i],y[i]-yy[i]);voidf4(double*y,double*x,doubl

4、e*yy,double*y0){y[0]=i.0；x[0]=0.0;yy[0]=1.0;for(inti=l;i<=9;i++){x[i]=x[i-l]+0.2;y0[i]=y[i-l]+0.2*(2*y[i-l]+x[i-l]*x[i-l]);y[i]=y[i-l]+0.1*(2*y0[i-l]+x[i-l]*x[i-l]+2*y[i-l]+x[i-l]*x[i-l]);yy[i]=-0.5*(x[i]*x[i]+x[i]+0.5)+1.25*exp(2*x[i]);printfC若x=%f,計算值是％f，真實值是%f，截斷誤差是%f?,,x[i],y[

5、i],yy[i]/y[i]-yy[i]);}}；voidmain(){doubley[10],x[10],yy[10]zy0[10];printf("常微分方程組是y-y-xandy(0)=2采用歐拉法求解得”)；fl(y,x,yy);printf("常微分方程組是y'=x*x+2yandy(0)=l釆用歐拉法求解得“);f2(y,x,yy);printf("常微分方程組是y-y-xandy(0)=2采用歐拉改進法求解得rT);fl(y,x,yy);printf("常微分方程組是y*=x*x+2yandy(0)=l采用歐拉改進法求解得

6、“);fl(y,x,yy);c：b'BACKUPl我的文檔I計算方法I計算方法C作業(yè)Qebug'歐拉公式解方程2心『常微分方程組是yit——acaaaaa0.400000,VJ0.600000,TI0.800000,VJ1.000000,VJ1.200000,yJ1.400080,y-

7、1?600000,7

8、1.800000,VIgx=0.200000,tj-右X-右X右X右X右X右x右x右X=y_x■丿古曰■丿古曰■丿古曰■丿古曰■丿古曰

9、

10、得03腿4fe1P2習(xí)4法2.拉E是是是是是是是是歐詐值值值?:?:值值值磊實實實實實實實實=2，0,0.0.0

11、.4,」7.0.>0000288180000639187<008385399V-042788895083841989d4n-23345567a2E是是是是是是是是最差差差差差差差差八折折折折圻折折hlnuhlnuhlnuhlnuhlmJhlMuhly鐵截截截截截截截截-0.021403-0.051825-0.094119-0.151941-0.229962-0.334133-0.472019-0.653215-0.889867常微分方程組是=x*x+2yandy<0>=l采用歐拉法求解若xW.200000"計算值是1.400000,直實值是1?494281」―

12、ifWr曰，'nS右x右X右X右x右x