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


M =0.

do i =1,ngl
m(i,i) =1._8
enddo


END subroutine gen_M
