The whole programming frame is efficient for the numerical solution of many P.D.E. problems (the stiffness matrix is symmetrical and positive) using Finite Element Method.