Thus D={C,pϕ} is conserved without re-solving the motion. Explicitly D=∂ϕC=pRpϕcosϕ−(pϕ2/R+A)sinϕ. Squaring and adding the two rotated components gives
Consequently the four quantities are not functionally independent: generically there are three independent constants. One may choose, locally, H,pϕ and the orientation encoded by C,D; their magnitude is fixed by the displayed relation. This also matches the maximum of three independent trajectory-labeling constants for a nonstationary autonomous flow in four-dimensional phase space. On singular subsets the local functional rank can be lower: for example, at A=0 and pϕ=0, the differentials of C,D are proportional to dpϕ, leaving at most two independent differentials together with dH; “three” is the generic count, not an assertion about every singular configuration.