close all;
%
xlim([-2.5 2]);
ylim([-1.5 1.5]);
plot(-2.5:0.05:2,(-2.5:0.05:2)-(-2.5:0.05:2).^3/3);
hold on
plot(-2.5:0.05:2,((-2.5:0.05:2)+0.7)./0.8);
%}
%
for i = 0:0.1:1
x0 = [-2, -1.1+i];
[t, phi] = ode45(@(t,x)fnode2(t,x),[0 50],x0);
plot(phi(:,1),phi(:,2),'r');
end
%}
