《改进欧拉法求解初值问题.doc》由会员分享,可在线阅读,更多相关《改进欧拉法求解初值问题.doc(1页珍藏版)》请在三一文库上搜索。
# include stdio.h# include math.h# include conio.h# define N 10float f(float x, float y)float s;s=y-2*x/y;return s;main ()float h,T1,T2,x=0.0,y=1.0,a=0.0,b=1.0;int i;h=(b-a)/N;printf(OUTPUT:N);for (i=0;i=N-1;i+)printf(k=%dtx%d=%fty%d=%fn,i,i,x,i,y);T1=y+h*f(x,y);x+=h; T2=y+h*f(x,T1);y=(T1+T2)/2;printf(k=%dtx%d=%fty%d=%fn,i,i,x,i,y);求解方程为的数值解。其中区间为N=10份