求这个二元微分方程组的MATLAB程序

假设x(1)=x1,x(2)=x1',x(3)=x2;把高阶导数换算成低阶倒数；得到3 个方程联立的方程组1：x(1)'=x(2) 2：x(2)'=-2*x(2)-2*x(3) 3：x(3)'=-x(2)-2x(1)+x(3)。求解程序如下：clear allclcf=@(t,x)[x(2);-2*x(2)-2*x(3);-x(2)-2*x(1)+x(3)];x0=[1;-1;0];[t,x]=ode45(f,[0,5],x0);plot(t,x)figure;plot(x(:,1),x(:,3))可以的到其时间曲线和相平面曲线。