// Batch code for SYS( 1)
module("PcGive");
package("PcGive", "Multiple-equation");
usedata("DHSY.xlsx");
system
{
    Y = c, y, D4p;
    Z = c_1, y_1, D4p_1, Trend;
    U = Constant, Seasonal, Seasonal_1, Seasonal_2;
}
estimate("OLS", 1959, 2, 1976, 2, 40);
forecast(40, 1, 1, 0, 1, 0, 0);
