New implementation of the FD-method for Sturm–Liouville problems with Dirichlet–Neumann boundary conditions

A new algorithm of the FD-method is proposed for solving Sturm–Liouville problems on an interval with Dirichlet–Neumann boundary conditions where the potential is a polynomial. A software implementation of the algorithm using a computer algebra software package shows its high efficiency.