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

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (4000)
>虫友互识 (565)
>导师招生 (401)
>文献求助 (182)
>考博 (126)
>学术会议 (95)
>硕博家园 (75)
>公派出国 (69)
>博后之家 (67)
>论文投稿 (62)
>招聘信息布告栏 (60)
>论文道贺祈福 (46)
>考研 (46)
>基金申请 (45)
>找工作 (40)
>休闲灌水 (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;
然后怎么运行得出结果呢？

回复此楼

» 猜你喜欢

售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急已经有5人回复
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急已经有7人回复
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急已经有11人回复
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急已经有8人回复
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急已经有17人回复
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急已经有16人回复
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急已经有4人回复
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急已经有5人回复
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急已经有6人回复
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有5人回复

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

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

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

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

laizhianzhi

铁虫 (初入文坛)

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

【答案】应助回帖

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

赞一下