用户名  找回密码
 注册
查看: 1413|回复: 0

经典龙格-库塔法

[复制链接]
发表于 2004-7-1 07:42:46 | 显示全部楼层 |阅读模式

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

您需要 登录 才可以下载或查看,没有账号?注册

x
已知
function R=rk4(f,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 j=0:m
    k1=feval(f,T(j),Y(j));
    k2=feval(f,T(j)+h/2,Y(j)+k1/2);
    k3=feval(f,T(j)+h/2,Y(j)+k2/2);
    k4=feval(f,T(j)+h,Y(j)+k3);
    Y(j+1)=Y(j)+h*(k1+2*k2+2*k3+k4)/6;
end
R=[T'Y']
如何求:
f=50*T(j)*T(j)+2*T(j)-50*Y(j)
谢谢
您需要登录后才可以回帖 登录 | 注册

本版积分规则

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