% RBC2BMS    computes the VARMA(1,1) matrices B1, M, Sig:
%
%               Y[t] = (B1+M) Y[t-1] + v[t] - M v[t-1],  Ev[t]v[t]'=Sig
%
%            for an economy with households solving
%
%               max E sum beta^t log (c[t]) + psi log(1-h[t])
%
%               st.  c[t]+x[t] = r[t]k[t]+w[t]h[t] - taxes
%                    k[t+1]    = (1-delta) k[t] + x[t]
%                    taxes     = tauk[t](r[t]-delta)k[t]
%                                  +tauh[t]w[t]h[t]
%                                  +taux[t]x[t]+tr[t]
%
%            and resource constraint c[t]+x[t]+g[t] = F(k[t],z[t]h[t]).
%            The user can choose elements of the data vector Y.
%

%            Ellen McGrattan, 4-27-06

or, give code to get log(l), log(y) etc.

then form matrices D0, D1, etc.

then get C1,C2 from D's

also need alpha

B1   = C1;
M    = C2*inv(C1)-C1;
Sig  = D0*Q*Q'*D0;


