来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/05/15 06:40:16
求用matlab解下面这个方程并绘图
其中,epsilon=[0.01,0.001]
参考代码:
E=[0.01,0.001];
clf,hold on
for e = E
dy = @(t,y)[y(2); e*y(2)^3-y(1)];
[t,y] = ode45(dy,[0 100],[1; 0]);
plot(t,y(:,1),'color',rand(1,3))
end
legend(arrayfun(@(e){sprintf('\\epsilon = %g',e)},E),0)