peak frequencies
fr=freson(h)
:h syslin list : :fr vector of peak frequencies in Hz :
returns the vector of peak frequencies in Hz for the SISO plant h
h=`syslin`_('c',-1+%s,(3+2*%s+%s^2)*(50+0.1*%s+%s^2))
fr=freson(h)
`bode`_(h)
g=20*`log`_(`abs`_(`repfreq`_(h,fr)))/`log`_(10)