--Fix bug in hessian function, which would fail to compute the Hessian if evaluation
point was int (instead of float).
--hessian and jacobian can take lists of numbers as evaluation points; they will
be converted to numpy arrays.
--Simplified notation in implementation of golden search method.
--Added Gauss-Jacobi and Gauss-Seidel solvers for linear equations.