برنامه بدست آوردن ژیزمان

clc
clear all
format long g
%barname gizman%
xa=input('xa=');ya=input('ya=');
xb=input('xb=');yb=input('yb=');
dxab=xb-xa;dyab=yb-ya;
vab=atand(abs(dxab/dyab));
if dxab>0&dyab>0
    gab=vab
elseif dxab>0&dyab<0
    gab=180-vab
elseif dxab<0&dyab<0
    gab=180+vab
elseif dxab<0&dyab>0
    gab=360-vab
end
if gab<=180
    gba=360-(180-gab)
elseif gab>180
    gba=180-(360-gab)
end