求解方程组，画相平面图

081270053

function dy=eq(t,y)
e=0.1;
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=-y(1)+e*(1-y(1)^2)*y(2);
end
建立m函数文件，y(1)是x，y（2）是y。
[t,y]=ode45('eq',[0 0.00001],[0 0]);
subplot(121);
plot(t,y(1,)
subplot(122);
plot(y(:,2),y(:,1),'*')
ε被写在函数里了。
大致是这样，个人不是特明白，仅会做而已。