anchor.diff:      !ftype = 2: eps
       !ftype = 4: png
       ftype = 2
  -    
  +
       n = hecMAT%N
       ndof = 3*n
       nnz = 9*n + 9*2*hecMAT%indexL(n)
  
positive.code:     rnum = (7.21d0+0.01*dlog10(dble(hecMAT%N)))*10.0d0/dble(hecMAT%N)
      cond = 1.0d0
      !rnum = (7.25d0)*10.0d0/dble(hecMAT%N)
      
      open(fio,file='nonzero.dat',status='replace')
        write(fio,"(a,f12.5,i)")"##magic number 10 : 7.2, ",rnum,hecMAT%N
        do i= 1, n
          jS= hecMAT%indexL(i-1) + 1
          jE= hecMAT%indexL(i  )
          write(fio,"(i,a,i)")i,"    ",i
          do j= jS, jE
            in = hecMAT%itemL(j)
            write(fio,"(i,a,i)")i,"    ",in
            write(fio,"(i,a,i)")in,"    ",i
          enddo
        enddo
      close(fio)
      
      open(fio,file='nonzero.plt',status='replace')
        if(ftype == 4)then
          write(fio,"(a)")'set terminal png'
  
negative.code:     implicit none
      type (hecmwST_local_mesh) :: hecMESH
      type (hecmwST_matrix)     :: hecMAT
  
      integer(kind=kint) :: i, j, in, jS, jE, ftype, n, ndof, nnz, fio
      real(kind=kreal) :: rnum, dens, cond
  
  