|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
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)
谢谢 |
|