function c=course(dx,dy);

if dx>=0 & dy>0
c=atan(dx/dy)*180/pi;
elseif dy<0 
c=atan(dx/dy)*180/pi+180;
elseif dx>0 & dy==0
c=90;
elseif dx<0 & dy==0
c=270;
elseif dx==0 & dy==0
c=0;
else
c=atan(dx/dy)*180/pi+360;
end



