Abstract: The main features of the method proposed in this paper are orthogonal, very efficient, easy control of the mesh spacing both in and directions and the upper and lower boundaries can be discontinuous.Several numerical examples are performed on IBM/PC-XT mini computer. The typical run time for a field with 39×17 grid points is 8 min, which is several times faster than that with the method in [2].