function traj = traje(t,x,flag,rad,omega)
% Computes time derivatives for use with script coriolis.m.
traj = zeros(4,1);
ind = 1;
traj(1) = (2.*omega+ind*x(1)/(rad*cos(x(4))))*sin(x(4))*x(2);
traj(2) = -(2.*omega+ind*x(1)/(rad*cos(x(4))))*sin(x(4))*x(1);
traj(3) = x(1)/(rad*cos(x(4)));
traj(4) = x(2)/rad;