anchor.diff:          h(1)=half*(tt*tt-tt)
           h(2)=half*(tt*tt+tt)
           h(3)=one-tt*tt
  -        dq(1,1)=-fourth*((one-ss)*(-rr-ss-one)+(one-rr)*(one-ss))
  -        dq(2,1)= fourth*((one-ss)*( rr-ss-one)+(one+rr)*(one-ss))
  -        dq(3,1)= fourth*((one+ss)*( rr+ss-one)+(one+rr)*(one+ss))
  -        dq(4,1)=-fourth*((one+ss)*(-rr+ss-one)+(one-rr)*(one+ss))
  -        dq(5,1)=-rr*(one-ss)
  -        dq(6,1)= half*(one-ss*ss)
  -        dq(7,1)=-rr*(one+ss)
  -        dq(8,1)=-half*(one-ss*ss)
  -        dq(1,2)=-fourth*((one-rr)*(-rr-ss-one)+(one-rr)*(one-ss))
  -        dq(2,2)=-fourth*((one+rr)*( rr-ss-one)+(one+rr)*(one-ss))
  -        dq(3,2)= fourth*((one+rr)*( rr+ss-one)+(one+rr)*(one+ss))
  -        dq(4,2)= fourth*((one-rr)*(-rr+ss-one)+(one-rr)*(one+ss))
  -        dq(5,2)=-half*(one-rr*rr)
  -        dq(6,2)=-ss*(one+rr)
  -        dq(7,2)= half*(one-rr*rr)
  -        dq(8,2)=-ss*(one-rr)
  +        dq(1,1)=(-fourth)*((one-ss)*(-rr-ss-one)+(one-rr)*(one-ss))
  +        dq(2,1)=( fourth)*((one-ss)*( rr-ss-one)+(one+rr)*(one-ss))
  +        dq(3,1)=( fourth)*((one+ss)*( rr+ss-one)+(one+rr)*(one+ss))
  +        dq(4,1)=(-fourth)*((one+ss)*(-rr+ss-one)+(one-rr)*(one+ss))
  +        dq(5,1)=(-rr)*(one-ss)
  +        dq(6,1)=( half)*(one-ss*ss)
  +        dq(7,1)=(-rr)*(one+ss)
  +        dq(8,1)=(-half)*(one-ss*ss)
  +        dq(1,2)=(-fourth)*((one-rr)*(-rr-ss-one)+(one-rr)*(one-ss))
  +        dq(2,2)=(-fourth)*((one+rr)*( rr-ss-one)+(one+rr)*(one-ss))
  +        dq(3,2)=( fourth)*((one+rr)*( rr+ss-one)+(one+rr)*(one+ss))
  +        dq(4,2)=( fourth)*((one-rr)*(-rr+ss-one)+(one-rr)*(one+ss))
  +        dq(5,2)=(-half)*(one-rr*rr)
  +        dq(6,2)=(-ss)*(one+rr)
  +        dq(7,2)=( half)*(one-rr*rr)
  +        dq(8,2)=(-ss)*(one-rr)
           dq(9,1)=zero
           dq(9,2)=zero
           dh(1)=half*(two*tt-one)
           dh(2)=half*(two*tt+one)
  -        dh(3)=-two*tt
  +        dh(3)=(-two)*tt
           do i=1,nen
             sh(4,i,l)=q(iq(i))*h(ih(i))
             sh(1,i,l)=dq(iq(i),1)*h(ih(i))
  
positive.code: c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  c
  c
        subroutine pquadpyr(sh,shj,gauss,ngauss,nen,nsd,nenmax,ngaussmax,
       & intord)
  c
  c... Subroutine to compute shape functions in natural coordinates,
  c    integration points, and weights for a quadratic pyramid.
  c
        include "implicit.inc"
  c
  c...  subroutine arguments
  c
        integer nsd,nenmax,ngaussmax,intord
        integer ngauss,nen
        double precision sh(nsd+1,nenmax,ngaussmax)
        double precision shj(nsd+1,nenmax,ngaussmax)
        double precision gauss(nsd+1,ngaussmax)
  c
  c...  defined constants
  c
        include "nconsts.inc"
        include "rconsts.inc"
  c
        double precision c1,c2,c3,c4,c5
        parameter(c1 =  15.0d0,
  
negative.code:       end
  c
  c version
  c $Id: pquadpyr.f,v 1.2 2004/07/07 19:48:45 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  