《改进欧拉法求解初值问题.docx》由会员分享,可在线阅读,更多相关《改进欧拉法求解初值问题.docx(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份;.