function ang = quadc(traj_x, traj_y)
ang = [];
for i =1:length(traj_x)
    temp_a = atand(abs(traj_y(i))/abs(traj_x(i)));
    if traj_x(i)> 0 & traj_y(i) > 0;
        ang = [ang temp_a];
    elseif traj_x(i)< 0 & traj_y(i) < 0;
        ang = [ang 360-temp_a];
    elseif traj_x(i)< 0 & traj_y(i) > 0;
        ang = [ang 180-temp_a];
    else   traj_x(i)> 0 & traj_y(i) < 0;
        ang = [ang 180+temp_a];
    end
end
ang = ang';