function Xprime=twodspring(t,X)

     global k;
     global m;
     Xprime=zeros(4,1);

     Xprime(1)= X(2); 
     Xprime(2)= (k/m)*( (1./sqrt(X(1).^2+X(3).^2)) - 1 ).*X(1);
     Xprime(3)= X(4); 
     Xprime(4)= -9.8 + (k/m)*( (1./sqrt(X(1).^2+X(3).^2)) - 1 ).*X(3);