This form uses the classical theory of beams in simple bending as represented by the equation M=-EJy''
The solution is based onto a polynomial (exact) formula obtained from strain energy V minimization (choose 'Show formulae' under 'Options' to see, relationships enforcing boundary conditions not shown).