24小时热门版块排行榜    

查看: 443  |  回复: 1

茗茗暖阳

铜虫 (小有名气)

[求助] 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;
然后怎么运行得出结果呢?
回复此楼

» 猜你喜欢

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

幸运从来不会光顾到没有奋斗的人身上!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

laizhianzhi

铁虫 (初入文坛)

【答案】应助回帖

你把这段代码复制进matlab脚本中,可以运行,在主界面显示结果。
2楼2015-08-27 02:05:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 茗茗暖阳 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 295材料求调剂,一志愿武汉理工085601专硕 +4 Charlieyq 2026-03-19 4/200 2026-03-20 14:26 by 无懈可击111
[考研] 本人考085602 化学工程 专硕 +18 不知道叫什么! 2026-03-15 20/1000 2026-03-20 13:52 by danney002
[考研] 265求调剂 +9 梁梁校校 2026-03-19 9/450 2026-03-20 12:33 by lature00
[考研] 279分求调剂 一志愿211 +9 chaojifeixia 2026-03-19 10/500 2026-03-20 12:29 by lature00
[考研] 274求调剂 +8 S.H1 2026-03-18 8/400 2026-03-20 11:53 by 学员8dgXkO
[考研] 工科材料085601 279求调剂 +6 困于星晨 2026-03-17 8/400 2026-03-20 11:24 by kkcoco25
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[考研] 330求调剂 +3 小材化本科 2026-03-18 3/150 2026-03-18 21:55 by 无懈可击111
[考研] 295求调剂 +3 一志愿京区211 2026-03-18 5/250 2026-03-18 17:03 by zhaoqian0518
[考研] 302求调剂 +10 呼呼呼。。。。 2026-03-17 10/500 2026-03-18 12:45 by Linda Hu
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
信息提示
请填处理意见