stopR = 10;
initx1 = 0;   inity1 = 0;
initvx1 = 0;   initvy1 = 0;   

initx2 = 1;    inity2 = 0;
initvx2 = 0;   initvy2 = 0; 

initx3 = 1.7;    inity3 = 0.7;
initvx3 = 0;   initvy3 = 0; 
sim('ThreeBody', 70);

%r12 = sqrt((x1 - x2) .^ 2 + (y1 - y2) .^ 2);
%r13 = sqrt((x1 - x3) .^ 2 + (y1 - y3) .^ 2);
%r32 = sqrt((x3 - x2) .^ 2 + (y3 - y2) .^ 2);

%plot(tout, r12, tout, r13, tout, r32)
%figure
%plot3(tout, x1, y1)
%xlabel('time')
%ylabel('x')
%zlabel('y')
%figure
%plot3(tout, x2, y2)
%figure
%plot3(tout, x3, y3)