err_steinmetz is a function. [err, Wfit, c] = err_steinmetz(X, fpoints, bpoints, wpoints, include_excess, normalize)