24小时热门版块排行榜

>论坛更新日志 (4762)
>休闲灌水 (1256)
>导师招生 (1039)
>虫友互识 (527)
>文献求助 (412)
>博后之家 (282)
>招聘信息布告栏 (201)
>论文投稿 (163)
>考博 (147)
>硕博家园 (125)
>找工作 (110)
>第一性原理 (94)
>绿色求助(高悬赏) (42)
>基金申请 (42)
>生物科学 (38)
>考研 (36)

返回列表

茗茗暖阳

铜虫 (小有名气)

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

回复此楼