anchor.diff:                  
       END SELECT
   
  -end subroutine HexTwnSeven
  +end subroutine HexTwn
   
   
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  
positive.code:     integer :: NumberOfElements, NumberOfMaterials, ElementGroupSize, QuadratureOrder=2
      integer :: N(11) !Pointers
      
      NPAR(5) = 27
      NPAR(4) = 0
      
      NumberOfElements = NPAR(2)
  
negative.code:     implicit none
      integer ::  ElementShapeNodes, i
      real(8) ::  NMatrix(3, 3*ElementShapeNodes), Transformed(3), N(ElementShapeNodes)
      real(8) :: xi, eta, zta
      
      xi = Transformed(1)
      eta = Transformed(2)
      zta = Transformed(3)
      select case (ElementShapeNodes)
      case (20)
      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>WORKING PROGRESS
          N(1) = (1-xi)*(1-eta)*(1-zta)*(-xi-eta-zta-2)/8
          N(2) = (1+xi)*(1-eta)*(1-zta)*(xi-eta-zta-2)/8
          N(3) = (1+xi)*(1+eta)*(1-zta)*(+xi+eta-zta-2)/8
          N(4) = (1-xi)*(1+eta)*(1-zta)*(-xi+eta-zta-2)/8
          N(5) = (1-xi)*(1-eta)*(1+zta)*(-xi-eta+zta-2)/8
          N(6) = (1+xi)*(1-eta)*(1+zta)*(xi-eta+zta-2)/8
          N(7) = (1+xi)*(1+eta)*(1+zta)*(xi+eta+zta-2)/8
          N(8) = (1-xi)*(1+eta)*(1+zta)*(-xi+eta+zta-2)/8
          N(9) = (1-xi**2)*(1-eta)*(1-zta)/4
          N(10) = (1-xi**2)*(1+eta)*(1-zta)/4
          N(11) = (1-xi**2)*(1+eta)*(1+zta)/4
          N(12) = (1-xi**2)*(1-eta)*(1+zta)/4
          N(13) = (1-eta**2)*(1-xi)*(1-zta)/4
          N(14) = (1-eta**2)*(1+xi)*(1-zta)/4
          N(15) = (1-eta**2)*(1+xi)*(1+zta)/4
          N(16) = (1-eta**2)*(1-xi)*(1+zta)/4
          N(17) = (1-zta**2)*(1-xi)*(1-eta)/4
          N(18) = (1-zta**2)*(1+xi)*(1-eta)/4
          N(19) = (1-zta**2)*(1+xi)*(1+eta)/4
          N(20) = (1-zta**2)*(1-xi)*(1+eta)/4
          NMatrix = reshape((/((/N(i), 0D0, 0D0, 0D0, N(i), 0D0, 0D0, 0D0, N(i)/),i=1,ElementShapeNodes)/), &
                            shape(NMatrix))
      end select
  