This form uses the classical theory of beams in simple bending as represented by the equation M=-EJy'' adapted for beams with circular axis.
The solution is based onto a polynomial-sinusoidal formula (exact) obtained from strain energy U minimization (hit the 'Formulas' button to see, relationships enforcing boundary conditions not shown).