/* complete this simple linear model */
decl ct, mx, vbeta, veps, vy;
ct = 4;                 // sample size T
mx = ...                // create $X$, where $x_t=t$
vbeta = ...             // 2 x 1 matrix with ones
veps = 0.1 * rann(ct, 1);
vy = ...                // create $Y$ 
print("y = X * beta + eps\n",
    "X", mx, "beta", vbeta, "eps", veps, "y", vy,
    "b_hat", (1/(mx'mx)) * (mx'vy) );
