This should now be faster; we are not computing every possible fit, but only the ones needed to reach the maxerror threshold
Added Vs, Vn in MA2; also, fixed order of x3 and x4 terms