Question

matlab求解以下微分方程组!

Answer 1

clear

clc

[t,x]=ode45(@(t,x)[-x(1)^2-x(2);x(1)-x(2)^3],[0 5],[1 0.5]);

plot(t,x(:,1),'r',t,x(:,2),'b');

legend('x(t)','y(t)');

axis([0,5,-10,2])