%  
N=1;
s0=1000*N;    % 
v0=10*N;        %    
a0=0.1*N;       %       
%   
c=(v0^2)/a0/2-s0;
b=2*v0;
a=a0;
%   
t1=(-b+sqrt(b^2-4*a*c))/(2*a)
t2=(v0+a0*t1)/a0
t0=t1+t2
% 
s1=v0*t1+(a0*t1^2)/2
s2=(a0*t2^2)/2
v=a0*t2
s = s1+s2
