function [ x_p ] = HW5_Solution_CentralPotential(X,A)
%CentralPotential first derivatives

x_p(1)=X(4);
x_p(2)=X(5);
x_p(3)=X(6);

r=sqrt( sum((X(1:3,:)).^2) );

x_p(4)=A*X(1)/r^3;
x_p(5)=A*X(2)/r^3;
x_p(6)=A*X(3)/r^3;

x_p=x_p';


end

