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