Oops, correct my faulty interpretation of leastsq's 'diag' argument in fit.py