function error = check(theta,t,t_0,vel)

t1  = 1/2 + 2*(t-t_0)*cos(theta) - vel*t; 
t2 = 2*(t-t_0)*sin(theta) - 1; 

error = [t1; t2];
