clc;clear;close;
[x,y]=meshgrid(linspace(-2,2));
h=streamslice(x,y,y,-2.*x-3.*y);
% h=streamslicetitle('相图')`
xlabel('x','Color','r');
ylabel('y','Color','r');
xlim([-2,2]);ylim([-2,2]);
set(h,'Color','k')
axis equal
hold on
% theta=0:pi/30:2*pi;
% x1=cos(theta);y1=sin(theta);
% plot(x1,y1,'r--')
% f = @(t,x)[x(2)+x(2).*(1-x(1).^2).*(x(2).^2-(x(1).^2).*(1-0.5*(1-0.5*x(1).^2)));
% x(1).*(1-x(1).^2)-x(2).*(x(2).^2-(x(1).^2).*(1-0.5*(1-0.5*x(1).^2)))]
% [t,x]=ode45(f,[0 100],[1 1]);
% plot(t,x(:,1),t,x(:,2)),