This function might be useful, e.g. for fitting in different norms. e.g. [see this discourse discussion](https://discourse.julialang.org/t/multivariate-polynomial-regression-of-discrete-data-in-l-infinity-norm/125369/7?u=stevengj)