24小时热门版块排行榜

返回列表

茗茗暖阳

铜虫 (小有名气)

应助: 6 (幼儿园)
金币: 396.1
散金: 43
红花: 1
帖子: 127
在线: 47.2小时
虫号: 3729976
注册: 2015-03-11
性别: MM
专业: 钢铁冶金

[求助] matlab编的4阶龙格库塔法解微分方程的程序已有1人参与

求matlab编的4阶龙格库塔法解微分方程的程序。。谢谢大家。最好是有具体的过程的。
m文件中写了如下内容：
function varargout=saxplaxliu(varargin)
clc,clear
x0=0;xn=1.2;y0=1;h=0.1;
[y,x]=lgkt4j(x0,xn,y0,h);
n=length(x);
fprintf(' i   x(i)   y(i)\n');
for i=1:n
    fprintf('%2d %4.4f %4.4f\n',i,x(i),y(i));
end
function z=f(x,y)
   z=-2*x*y^2;
function [y,x]=lgkt4j(x0,xn,y0,h)
x=x0:h:xn;
n=length(x);
y1=x;
y1(1)=y0;
for i=1:n-1
    K1=f(x(i),y1(i));
    K2=f(x(i)+h/2,y1(i)+h/2*K1);
    K3=f(x(i)+h/2,y1(i)+h/2*K2);
    K4=f(x(i)+h,y1(i)+h*K3);
    y1(i+1)=y1(i)+h/6*(K1+2*K2+2*K3+K4);
end
   y=y1;
然后怎么运行得出结果呢？

回复此楼

» 猜你喜欢

谈谈两天一夜的“延安行” 已经有15人回复
EST投稿状态问题已经有6人回复
职称评审没过，求安慰已经有15人回复
垃圾破二本职称评审标准已经有11人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有16人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
聘U V热熔胶研究人员已经有10人回复
求助文献已经有3人回复
投稿返修后收到这样的回复，还有希望吗已经有8人回复
三无产品还有机会吗已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

幸运从来不会光顾到没有奋斗的人身上！

1楼2015-03-19 20:07:00

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

laizhianzhi

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 243.7
红花: 1
帖子: 23
在线: 13.9小时
虫号: 1519959
注册: 2011-12-02
专业: 光学

【答案】应助回帖

相关版块跳转我要订阅楼主茗茗暖阳的主题更新

返回列表