figure
axis equal off
hold on
t=linspace(0,2*pi);
planet1 = struct('A',10,'P',2,'phi', pi/8); 
planet2 = struct('A', 4,'P',1,'phi',-pi/7);
pLoc1=orbit(t,planet1,'b'); pt1=orbit(3,planet1); % Radan 1 piste, kun t=3
pLoc2=orbit(t,planet2,'r'); pt2=orbit(1,planet2); % Radan 2 piste, kun t=1
% Katsotaan:
pt1
pt1.x

plot([pt1.x pt2.x],[pt1.y pt2.y],'o');
plot([pt1.x pt2.x],[pt1.y pt2.y],'g--');
plot([0 0],[0,0],'*y')
text(0,-.2,'Aurinko')
hold off

% Tässä on mahdollisuus kokeilla parametrien A,P ja phi vaikutusta rataan.
