“…The solution of such implicit schemes in most instances is extremely difficult, as it reduces ultimately to the solution of simultaneous equations involving a large number of unknowns interrelated in a manner quite unsuitable for easy determination. If the system is linear, there are several classical methods by which it may be treated; most notably, by direct inversion of matrices of large size [1], by iterative procedures [2] and by relaxation [3]. It is the purpose of this paper to introduce a finite difference technique especially adapted to the solution of boundary value problems.…”