close all
clear all
figure(1)
hold on 
for theta=[0:50]*pi/25
    x0=1e-2*[cos(theta);sin(theta)];
    [t,x]=ode45(@dxdt1,[0 10],x0);
    plot(x(:,1),x(:,2))
end


figure(2)
hold on 
for theta=[-50:50]/10
    x0=[theta 2];
    [t,x]=ode45(@dxdt1,[0 20],x0);
    plot(x(:,1),x(:,2))
    x0=[theta -2];
    [t,x]=ode45(@dxdt1,[0 20],x0);
    plot(x(:,1),x(:,2))
end
axis([-10 10 -2 2])