matlab解微分方程组

解微分方程组：
x’=-x^3-y,x(0)=1
y’=x-y^3,y(0)=0.5 0<t<30

编辑器窗口：
%M函数eg6_3fun.m
function f=eg6_3fun(t,x)
f(1)=-x(1)^3-x(2);
f(2)=x(1)-x(2)^3;
f=f(😃;%注意要保证f为列向量

指令窗口：

[t,x]=ode45(@eg6_3fun,[0 30],[1;0.5] %1和0.5是初始值

解应该是x（t），y（t）两个函数
所以
把x和y都写成一个x
x表示一个向量
x----x(1)
y----x(2）

算法本身要求f是列向量，
所以，f=f（:）