kill(all);

amp(f_ratio,zeta):=1/sqrt((1-f_ratio^2)^2+(2*zeta*f_ratio)^2);

wxplot2d([amp(f_ratio,0),1],[f_ratio,0,2],[y,0,4]);

solve(amp(f_ratio,0)^2=1.1^2,f_ratio);
%[2],numer;

solve(amp(f_ratio,0)^2=1.05^2,f_ratio);
%[2],numer;

solve(amp(f_ratio,0)^2=1.02^2,f_ratio);
%[4],numer;

solve(amp(f_ratio,0)^2=1.01^2,f_ratio);
%[2],numer;