这个非线性方程怎么解啊？

yaoyaolan

大家帮忙看看：
function dz=yaoyao(t,z)
% z(1)=x1;z(2)=x2;z(3)=x3;z(4)=y1;z(5)=y2;z(6)=y3;
alpha=10;beta=16;b=-0.143;
dz=zeros(6,1);
e(1)=z(4)-z(1);
e(2)=z(5)-z(2);
e(3)=z(6)-z(3);
U1=e(3)-(beta*beta+alpha+1)*e(2)-(beta-alpha*b+1)*e(1)+alpha*e(1)*(e(1).^2+3*z(1)*e(1)+3*z(1).^2);
dz(1)=alpha*(z(2)-b*z(1)-z(1).^3);
dz(2)=z(1)-z(2)+z(3);
dz(3)=-beta*z(2);
dz(4)=alpha*(z(5)-b*z(4)-z(4).^3)+U1;
dz(5)=z(4)-z(5)+z(6);
dz(6)=-beta*z(5);

clear;
[T,z]=ode45('yaoyao',[0 30],[1 2 1 10 5 5]);
n=length(T);
for i=1:n
    E1(i)=z(i,4)-z(i,1);
    E2(i)=z(i,5)-z(i,2);
    E3(i)=z(i,6)-z(i,3);
end
figure
subplot(2,3,1)
plot(T,z(:,1));xlabel('T/s');ylabel('x1');grid;
subplot(2,3,2)
plot(T,z(:,2));xlabel('T/s');ylabel('x2');grid;
subplot(2,3,3)
plot(T,z(:,3));xlabel('T/s');ylabel('x3');grid;
subplot(2,3,4)
plot(T,z(:,4));xlabel('T/s');ylabel('y1');grid;
subplot(2,3,5)
plot(T,z(:,5));xlabel('T/s');ylabel('y2');grid;
subplot(2,3,6)
plot(T,z(:,6));xlabel('T/s');ylabel('y3');grid;
figure
subplot(2,2,1)
plot(T,E1);xlabel('T/s');ylabel('e1');grid;
subplot(2,2,2)
plot(T,E2);xlabel('T/s');ylabel('e2');grid;
subplot(2,2,3)
plot(T,E3);xlabel('T/s');ylabel('e3');grid;
这是以下方程的程序：

dx(1)/dt=alpha(x(2)-bx(1)-x(1).^3),
dx(2)/dt=x(1)-x(2)-x(3),
dx(3)/dt=-betax(2);

dy(1)=alpha(y(2)-by(1)-y(1).^3)+U(1),
dy(2)/dt=y(1)-y(2)-y(3),
dy(3)/dt=-betay(3).
其中，
U(1)=y(3)-x(3)-(beta^2+alpha+1)(y(2)-x(2))-(beta-alpha*b+1)(y(1)-x(1))+alpha*(y(1)-x(1))((y(1)-x(1))^2+3x(1)(y(1)-x(1))+x(1)^2).

我想要画出y(1)-x(1), y(2)-x(2), y(3)-x(3)的图像

yaoyaolan

是y(1)减去x(1)，y(2)减去x(2),y(3)减去x(3)的图像