function [q]=ptp(beta_e,beta_a)

%Parameter
 v_max=double(10);
 a_max=double(5);

 global params
 params  = [beta_e; beta_a; v_max; a_max;];

 delta_q=beta_a-beta_e;
 t_b=v_max/a_max;
 t_e=delta_q/v_max+v_max/a_max;
 t_v=(t_e-t_b);

%  q=step_konstant(t);
q=v_max*t-0.5*v_max^2/a_max;


end