The subroutines elmstfMRT, elmstfDKT, sm3mb, and sm3mh come from
the file assm325.f of the STADYN finite element program 
(http://roger.ecn.purdue.edu/~jfdoyle/code.html) by James F. Doyle.

The subroutines are included with permission of the author.
