function omega_u_ar=omega_uar_functie(omega_0,b,m);

omega_u_ar=sqrt(omega_0^2-(b^2/(2*m^2)));
