function A=A_functie(F,m,omega_u,omega_0,b);

A=F./sqrt(m.^2.*(omega_u.^2-omega_0.^2).^2 + b.^2.*omega_u.^2);
