solve u= (-m*sqrt(r^2-u^2)-m*u)/(m^2+1) for m