
clear all

ModelName='centereq'; 
options=odeset('e','off','BDF','on','AbsTol',1e-8,'RelTol',1e-8); % for solve ode, use =[] for default.

tran=8;
tend=30;

X0=[-1 1];
param=[1 1];

y0=[1,.667,.333,0,-1,-.667,-.333];

figure(10)
hold off

plot([-1.2 1.2], [0 0], '--','color',[.7 .7 .7])
hold on
plot([0 0], [-1.3 1.3], '--','color',[.7 .7 .7])

plot([-.4 1],[0,0],'r-','linewidth',2);

plot([0 0],[-1.2,1.2],'b-','linewidth',2);

plot([0], [0], 'ko','markerfacecolor','w','markersize',10)

for i=1:length(y0)
xinit=[-1,y0(i)];
[~,X]=ode15s(ModelName, [0 tran], xinit, options, param);
plot(X(:,1),X(:,2),'r-','linewidth',2);
end
text(1.2,-.1,'x','fontsize',14)
text(.1,1.2,'y','fontsize',14)
text(.1,-.15,'0','fontsize',14)
xlabel('x','fontsize',16)
ylabel('y','fontsize',16)

axis off

%%
