版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

>论坛更新日志 (1462)
>论文投稿 (108)
>学术会议 (100)
>导师招生 (31)
>虫友互识 (30)
>硕博家园 (26)
>考博 (25)
>基金申请 (22)
>考研 (13)
>公派出国 (12)
>博后之家 (10)
>文献求助 (10)
>教师之家 (8)
>找工作 (8)
>论文道贺祈福 (7)
>休闲灌水 (7)

返回列表

茗茗暖阳

铜虫 (小有名气)

应助: 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;
然后怎么运行得出结果呢？

回复此楼

» 猜你喜欢

澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域已经有56人回复
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人已经有15人回复
欢迎发来filecode的Mz6后的代码验证其规律已经有29人回复
关于Filecode分析方法已经有12人回复
咱们一起用铁证分析2026国家社科基金中标与否已经有11人回复
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊已经有10人回复
阿姨已经有4人回复
一作与独作在应聘高校教师时区别大吗已经有4人回复
静等基金结果已经有21人回复
FileCode能看出啥？已经有28人回复

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

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

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

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

laizhianzhi

铁虫 (初入文坛)

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

【答案】应助回帖

你把这段代码复制进matlab脚本中，可以运行，在主界面显示结果。

赞一下