A method of solution of linear differential equations is derived which requires no special starting procedure, and which is, in principle at least, capable of extension to any order of accuracy. Further advantages of the method are that boundary conditions are easy to satisfy, and that the solution of families of equations differing only in the input function requires little extra computation. The method requires the reciprocation of a matrix of order equal to the number of points for which a solution is to be obtained, but this can be much simplified by suitably partitioning the matrix, and with most automatic high-speed computers the reciprocation of a matrix is in any case an easily programmed operation.