Wrote a newton-raphson solver skeleton.  It could be cleaned up, but not much.  Also wrote vectorized initializers for gradient and inverse step size.  Still need to give each model type N_parameters.