A pascal program for the solution of sparse linear systems – Part I