anchor.diff:     ry = rho*Iyy;
      po = rho*Jzz*L;   
   
  -   M(1:12,1:12) = 0
  +   M(1:12,1:12) = 0.0_ReKi
         
  -   M( 9,  9) = t/3.0
  -   M( 7,  7) = 13.0*t/35.0 + 6.0*ry/(5.0*L)
  -   M( 8,  8) = 13.0*t/35.0 + 6.0*rx/(5.0*L)
  -   M(12, 12) = po/3.0
  -   M(10, 10) = t*L*L/105.0 + 2.0*L*rx/15.0
  -   M(11, 11) = t*L*L/105.0 + 2.0*L*ry/15.0
  -   M( 2,  4) = -11.0*t*L/210.0 - rx/10.0
  -   M( 1,  5) =  11.0*t*L/210.0 + ry/10.0
  -   M( 3,  9) = t/6.0
  -   M( 5,  7) =  13.*t*L/420. - ry/10.
  -   M( 4,  8) = -13.*t*L/420. + rx/10. 
  -   M( 6, 12) = po/6.
  -   M( 2, 10) =  13.*t*L/420. - rx/10. 
  -   M( 1, 11) = -13.*t*L/420. + ry/10.
  -   M( 8, 10) =  11.*t*L/210. + rx/10.
  -   M( 7, 11) = -11.*t*L/210. - ry/10. 
  -   M( 1,  7) =  9.*t/70. - 6.*ry/(5.*L)
  -   M( 2,  8) =  9.*t/70. - 6.*rx/(5.*L)
  -   M( 4, 10) = -L*L*t/140. - rx*L/30. 
  -   M( 5, 11) = -L*L*t/140. - ry*L/30.
  +   M( 9,  9) = t/3.0_ReKi
  +   M( 7,  7) = 13.0_ReKi*t/35.0_ReKi + 6.0_ReKi*ry/(5.0_ReKi*L)
  +   M( 8,  8) = 13.0_ReKi*t/35.0_ReKi + 6.0_ReKi*rx/(5.0_ReKi*L)
  +   M(12, 12) = po/3.0_ReKi
  +   M(10, 10) = t*L*L/105.0_ReKi + 2.0_ReKi*L*rx/15.0_ReKi
  +   M(11, 11) = t*L*L/105.0_ReKi + 2.0_ReKi*L*ry/15.0_ReKi
  +   M( 2,  4) = -11.0_ReKi*t*L/210.0_ReKi - rx/10.0_ReKi
  +   M( 1,  5) =  11.0_ReKi*t*L/210.0_ReKi + ry/10.0_ReKi
  +   M( 3,  9) = t/6.0_ReKi
  +   M( 5,  7) =  13._ReKi*t*L/420._ReKi - ry/10._ReKi
  +   M( 4,  8) = -13._ReKi*t*L/420._ReKi + rx/10._ReKi
  +   M( 6, 12) = po/6._ReKi
  +   M( 2, 10) =  13._ReKi*t*L/420._ReKi - rx/10._ReKi
  +   M( 1, 11) = -13._ReKi*t*L/420._ReKi + ry/10._ReKi
  +   M( 8, 10) =  11._ReKi*t*L/210._ReKi + rx/10._ReKi
  +   M( 7, 11) = -11._ReKi*t*L/210._ReKi - ry/10._ReKi
  +   M( 1,  7) =  9._ReKi*t/70._ReKi - 6._ReKi*ry/(5._ReKi*L)
  +   M( 2,  8) =  9._ReKi*t/70._ReKi - 6._ReKi*rx/(5._ReKi*L)
  +   M( 4, 10) = -L*L*t/140._ReKi - rx*L/30._ReKi 
  +   M( 5, 11) = -L*L*t/140._ReKi - ry*L/30._ReKi
      
      M( 3,  3) = M( 9,  9)
      M( 1,  1) = M( 7,  7)
  
positive.code:    M( 8,  2) =  9._ReKi/70._ReKi * t
     M( 9,  3) = t/6.0_ReKi
     
     DC = 0
     DC( 1: 3,  1: 3) = DirCos
     DC( 4: 6,  4: 6) = DirCos
     DC( 7: 9,  7: 9) = DirCos
  
negative.code:    Ax = kappa*A
     Ay = kappa*A
     
     K(1:12,1:12) = 0.0_ReKi
     
     IF (Shear) THEN
        Kx = 12.0_ReKi*E*Iyy / (G*Ax*L*L)
        Ky = 12.0_ReKi*E*Ixx / (G*Ay*L*L)
     ELSE
        Kx = 0.0_ReKi
        Ky = 0.0_ReKi
     ENDIF
        
     K( 9,  9) = E*A/L
     K( 7,  7) = 12.0_ReKi*E*Iyy/( L*L*L*(1.0_ReKi + Kx) )
     K( 8,  8) = 12.0_ReKi*E*Ixx/( L*L*L*(1.0_ReKi + Ky) )
     K(12, 12) = G*Jzz/L
     K(10, 10) = (4.0_ReKi + Ky)*E*Ixx / ( L*(1.0_ReKi+Ky) )  
     K(11, 11) = (4.0_ReKi + Kx)*E*Iyy / ( L*(1.0_ReKi+Kx) )
     K( 2,  4) = -6._ReKi*E*Ixx / ( L*L*(1.0_ReKi+Ky) )
     K( 1,  5) =  6._ReKi*E*Iyy / ( L*L*(1.0_ReKi+Kx) )
     K( 4, 10) = (2.0_ReKi-Ky)*E*Ixx / ( L*(1.0_ReKi+Ky) )
     K( 5, 11) = (2.0_ReKi-Kx)*E*Iyy / ( L*(1.0_ReKi+Kx) )
     
     K( 3,  3)  = K(9,9)
     K( 1,  1)  = K(7,7)
  