###
#
#octave_control_bode_phase.m
#https://octave-online.blogspot.jp/search/label/control_bode_phase
#
###
function out=octave_control_bode_phase( num, den, min, max, n )
 w = logspace( min, max, n  );
 s = polyval(num, complex(0, w))./polyval(den, complex(0, w));
 phase = angle(s)/3.14159*180;
 semilogx(w, phase);
 out = "";
endfunction

