Todo List
void chiSquared( double[][] $x, double[] $a, $y, $s, object $f, double[] $y,, double[] $s,)
Chi-squared is the distribution of squared Gaussian errors, thus the name.
the solve( x $x, y $a, a $y, vary $s, s2 $vary, lambda $f, termepsilon $lambda, maxiter $termepsilon, verbose $maxiter, $verbose)
Note that LMfunc implements the value and gradient of f(x,a), NOT the value and gradient of E with respect to a!