  We apply a new calculation scheme of a finite element method (FEM) for
solving an elliptic boundary-value problem describing a quadrupole vibration
collective nuclear model with tetrahedral symmetry. We use of shape functions
constructed with interpolation Lagrange polynomials on a triangle finite
element grid and compare the FEM results with obtained early by a finite
difference method.
