VARXFit function

Fit a VAR or VARX model by least squares