make_irls_warmup function

Create starting values for iterative reweighted least squares