声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1945|回复: 1

[综合讨论] 请问怎么可以调出并察看matlab中ode23或ode45的源程序

[复制链接]
发表于 2006-12-5 09:23 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
这是我自己编的一段关于3阶Runge-Kutta方法的matlab程序,中间好像有错误,请求指导,谢谢
function E=rk3(a,b,ya,m)
% a,b区间两端点值
% ya为给定某一点的函数值
% m最大迭代次数
h=(b-a)/m;
T=zeros(1,m+1);
Y=zeros(1,m+1);
T=a:h:b;
Y(1)=ya;
for y=1:m;
    k1=f(T(j),Y(j));
    k2=f(T(j)+0.5*h,Y(j)+0.5*h*k1);
    k3=f(T(j)+h,Y(j)-h*k1+2*h*k2);
    Y(j+1)=Y(j)+1/6*h*(k1+4*k2+k3);
end
E=[T',Y'];
回复
分享到:

使用道具 举报

发表于 2006-12-5 10:55 | 显示全部楼层
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-9-24 23:30 , Processed in 0.056730 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表