% main1.m 
y0 = [10 0]';
[t y] = ode45('f',[0 100],y0);
clf; hold on;
axis([-6 10 -6 10]);
axis('square');
plot(y(:,1),y(:,2));

t = 0:0.05:6.3
[X, Xs, Y, Ys] = child(t);
plot(X,Y,':')
hold off;

