program assgin8_1
!
use denseop; use sparsem; use sparseop
integer,parameter :: n=3
integer :: i,j
type (densem) :: x
type (sparse_hashm)::xx
type (sparse_ija) :: y

real(8)::a(n,n)=reshape( (/36,30,18,30,41,23,18,23,14/),(/n,n/) ),&
b(n)=(/150,181,106/)

call printmat(a,'A matrix')

call zerom(x,n)
call zerom(xx,n)
do i=1,n
do j=1,n
call addm(a(i,j),i,j,x)
enddo
enddo

call printm(x)

end program assgin8_1