y0=[-4 -3 5 -6]; 
x0=0;  %起始時間點
xf=200;%終止時間點
[t,ym]=ode45('HyChen',x0:0.001:xf,y0);
figure(1)
subplot(221)
plot(ym(:,1),ym(:,2));
grid on
xlabel('y1');ylabel('y2');
subplot(222)
plot(ym(:,1),ym(:,3));
grid on
xlabel('y1');ylabel('y3');
subplot(223)
plot(ym(:,1),ym(:,4));
grid on
xlabel('y1');ylabel('y4');
subplot(224)
plot(ym(:,2),ym(:,3));
grid on
xlabel('y2');ylabel('y3');
figure(2)
subplot(221)
plot(ym(:,2),ym(:,4));
grid on
xlabel('y2');ylabel('y4');

subplot(222)
plot(ym(:,3),ym(:,4));
grid on
xlabel('y3');ylabel('y4');

subplot(223)
plot3(ym(:,1),ym(:,2),ym(:,3));
grid on
xlabel('y1');ylabel('y2');zlabel('y3')

subplot(224)
plot3(ym(:,1),ym(:,2),ym(:,4));
grid on
xlabel('y1');ylabel('y2');zlabel('y4')