I am currently working on my big honours project and need some help on a mainly coding problem I may be having. So basically in the images attached is the Lagrangian density for my soliton that I must solve numerically, initial I thought my best approach would be to implement spectral methods using Fourier transforms as that is one of the more accurate methods of solving BVPs with smooth solutions on periodic domains. Now when discussing with my supervisor he told me I was simply overcomplicating it, because we were working in one space dimension and because of the property of solitons being constant with respect to time, we can do some “Math trickier” to reduce the second order ODE to a simple integral and then find the inverse of the function gained from that integral.
Lagrangian density and converting Euler Lagrange equation to an integral
enter image description here

Reduced differential transform method