The numerical solution of boundary value problems for linear systems of first order equations with a regular singular point at one endpoint is considered. The standard procedure of expanding about the singularity to get a nonsingular problem over a reduced interval is justified in some detail. Quite general boundary conditions are included which permit unbounded solutions. Error estimates are given and some numerical calculations are presented to check the theory.