d(sqrt((u/M)^2-(P/(M^2))^2))/du