yexl 发表于 2006-12-5 09:23

请问怎么可以调出并察看matlab中ode23或ode45的源程序

这是我自己编的一段关于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=;

happy 发表于 2006-12-5 10:55

http://forum.vibunion.com/forum/thread-17615-1-1.html

参考这里吧
页: [1]
查看完整版本: 请问怎么可以调出并察看matlab中ode23或ode45的源程序