function ddt=RHS2(y,x)

y1=y(1);
y2=y(2);
y3=y(3);
y4=y(4);

ddt=[y2; y3; y4; 1-x^2];