subroutine gen_M()
use work
! Genereaza matricea de masa M(ngl, ngl)
implicit real(8) (a-h, o-z)

! ngl =3

M =0.

m(1,1) =1.
m(2,2) =1.
m(3,3) =1

M =M*40_8

END subroutine gen_M
