untitled6.m

来自「matlab-M文件库 可用于学习Matlab 标准化编程」· M 代码 · 共 9 行

y0=[1;1;0.5;1;1];
dt=0.1;t=0;

y=y0;
for n=1:5
t=t+dt;
f=[1 - t*y(1), exp(-t) - 3*y(2), t^2-y(3), -y*abs(y), sin(t)-sqrt(abs(y))];
y=y+dt*f'
end