function plawn_faza;
a=[1 1];
b=[0 1];
[h,w]=freqs(b,a);
phi=angle(h);
phi=unwrap(phi);
plot(w,phi*180/pi);
