We use Legendre function expansion methods to find numerical solutions of linear ordinary differential equations on the interval and functional differential equations and where and are given, and
