The conservative force term of the Euler-Lagrange equation is therefore:
∂s∂L=(rsummit−∫0ssinθ(σ)dσ)2GMmsinθ(s). Hence the functional derivative of S with respect to s is therefore:
δsδS=ms¨−(rsummit−∫0ssinθ(σ)dσ)2GMmsinθ(s). The time derivative of generalized momentum is the same as above. As for Qs it is:
Qs=−(rsummit−∫0ssinθ(σ)dσ)2GMmCrrcosθ(s)−2ρCDAs˙2. Substituting Equations (1) and (2) into Equation (??) yields:
ms¨−(rsummit−∫0ssinθ(σ)dσ)2GMmsinθ(s)=−(rsummit−∫0ssinθ(σ)dσ)2GMmCrrcosθ(s)−2ρCDAs˙2. Rearranging and factorizing:
ms¨=(rsummit−∫0ssinθ(σ)dσ)2GMm[sinθ(s)−Crrcosθ(s)]−2ρCDAs˙2. Dividing by m:
s¨=(rsummit−∫0ssinθ(σ)dσ)2GM[sinθ(s)−Crrcosθ(s)]−2mρCDAs˙2.