syms x
A=[40;50;60;70;80;90;]
for ii = 1:length(A)
[t,x]=ode45(@odefun,[0 1],0.2,A(ii));
tr{ii} = t;
xr{ii} = x;
end
% 绘图的代码还需要自己修改
plot(t,x)
xlabel('t')
ylabel('x')
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!