I am a new comer to R circle. I am trying to simulate a VAR estimation. The problem is given above in the image. Assume that the errors are iid and normally distributed. Here the number of observations x1=x2=64.
I've written the code below. Not sure if it is correct.