改进欧拉法求解初值问题.doc

《改进欧拉法求解初值问题.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、#include"stdio.h"#include"math.h"#include"conio.h"#defineN10floatf(floatx,floaty){floats;s=y-2*x/y;returns;}main(){floath,T1,T2,x=0.0,y=1.0,a=0.0,b=1.0;inti;h=(b-a)/N;printf("OUTPUT:N");for(i=0;i<=N-1;i++){printf("k=%dtx%d=%fty%d=%f",i,i,x,i,y);T

2、1=y+h*f(x,y);x+=h;T2=y+h*f(x,T1);y=(T1+T2)/2;}printf("k=%dtx%d=%fty%d=%f",i,i,x,i,y);}求解方程为的数值解。其中区间为N=10份