改进欧拉法求解初值问题.docx
.# 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份;.