anchor.diff:       1             pnt(2,1),pnt(2,2),pnt(2,3),
        2             pnt(3,1),pnt(3,2),pnt(3,3),
        3             eig(1),eig(2),eig(3)
  -      mxx = eig(1)*pnt(1,1)*pnt(1,1) + eig(2)*pnt(2,1)*pnt(2,1) + 
  +      mxx = eig(1)*pnt(1,1)*pnt(1,1) + eig(2)*pnt(2,1)*pnt(2,1) +
        1                                          eig(3)*pnt(3,1)*pnt(3,1)
         myy = eig(1)*pnt(1,2)*pnt(1,2) + eig(2)*pnt(2,2)*pnt(2,2) +
        1                                          eig(3)*pnt(3,2)*pnt(3,2)
  
positive.code:       t(1) = a(1,3)
        t(2) = a(2,3)
        t(3) = a(3,3)
  C Write P, N, and T vectors
        write(soln,1001) p(1),p(2),p(3),o(1),o(2),o(3),t(1),t(2),t(3),
       1                 w(1),w(2),w(3)
   1001 format(9F7.3,3(1P1E14.6))
        RETURN
        END
  
negative.code: C----
        IMPLICIT none
        CHARACTER*180 line,soln
  #ifdef USELAPACK
        INTEGER n,lda,lwmax
        PARAMETER (n=3,lda=n,lwmax=250)
        INTEGER info,lwork
        REAL*8 a(lda,n),w(n),work(lwmax)
        REAL*8 m(6)
        REAL*8 p(3),o(3),t(3)
  #endif
  
  C There seems to be some confusion in my mind as to the definitions of r, t, and p
  
  