function Yp=kakskpl(t,Y)
% Kahden kappaleen ongelma. Gravitaatiovakio = 1
global m1 m2
R1=Y(1:2);
V1=Y(3:4);
R2=Y(5:6);
V2=Y(7:8);

r=norm(R2-R1);
Yp=[V1;
    m2*(R2-R1)/r^3;
    V2;
    m1*(R1-R2)/r^3];

