var_main now has an option include_squared_trend to include the square of a trend variable when a trend variable is included for a model. var_main now only includes squared residuals when exogenous_max_iterations = 3.