基本信息
源码名称：四阶四级龙格库塔法数值求解常微分方程
源码大小：1.07KB
文件格式：.m
开发语言：MATLAB
更新时间：2020-08-15
友情提示：（无需注册或充值，赞助后即可获取资源下载链接）
嘿，亲！知识可是无价之宝呢，但咱这精心整理的资料也耗费了不少心血呀。小小地破费一下，绝对物超所值哦！如有下载和支付问题，请联系我们QQ(微信同号)：813200300
本次赞助数额为： 2 元×
微信扫码支付：2 元
×
请留下您的邮箱，我们将在2小时内将文件发到您的邮箱
源码介绍
应用四阶四级龙格库塔法求解常微分方程的数值解
for x=0:h:5-h
k1=dyx(x,y1(i-1));
k2=dyx(x h/2,y1(i-1) h/2*k1);
k3=dyx(x h/2,y1(i-1) h/2*k2);
k4=dyx(x h,y1(i-1) h*k3);
y1(i)=y1(i-1) h/6*(k1 2*k2 2*k3 k4);
i=i 1;
end %公式
i=2;
for x=0:-h:-5 h
k1=dyx(x,y2(i-1));
k2=dyx(x-h/2,y2(i-1)-h/2*k1);
k3=dyx(x-h/2,y2(i-1)-h/2*k2);
k4=dyx(x-h,y2(i-1)-h*k3);
y2(i)=y2(i-1)-h/6*(k1 2*k2 2*k3 k4);
i=i 1;
end
应用四阶四级龙格库塔法求解常微分方程的数值解
for x=0:h:5-h
k1=dyx(x,y1(i-1));
k2=dyx(x h/2,y1(i-1) h/2*k1);
k3=dyx(x h/2,y1(i-1) h/2*k2);
k4=dyx(x h,y1(i-1) h*k3);
y1(i)=y1(i-1) h/6*(k1 2*k2 2*k3 k4);
i=i 1;
end %公式
i=2;
for x=0:-h:-5 h
k1=dyx(x,y2(i-1));
k2=dyx(x-h/2,y2(i-1)-h/2*k1);
k3=dyx(x-h/2,y2(i-1)-h/2*k2);
k4=dyx(x-h,y2(i-1)-h*k3);
y2(i)=y2(i-1)-h/6*(k1 2*k2 2*k3 k4);
i=i 1;
end