改进欧拉法求解初值问题
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
# include ""
# include ""
# include ""
# define N 10
float f(float x, float y)
{
float s;
s=y-2*x/y;
return s;
}
main ()
{
float h,T1,T2,x=,y=,a=,b=;
int i;
h=(b-a)/N;
printf("OUTPUT:\N");
for (i=0;i<=N-1;i++)
{
printf("k=%d\tx%d=%f\ty%d=%f\n",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=%d\tx%d=%f\ty%d=%f\n",i,i,x,i,y); }
求解方程为⎩⎨
⎧=<=<=-='1)0()10(/2y x y x y y 的数值解。其中区间为N=10份