function F=fase_functie(b,omega_u,m,omega_0);

F=atan(b.*omega_u./(m.*(omega_0.^2-omega_u.^2)))/(2*pi)*360;
