second derivative writeup; removed cruft from *-fns.R
switched "beta" to "exp(beta)".
can now compute derivatives and Gram matrix to identify collinear parameters