A numerical method for solving systems of higher order linear functional differential equations