% compute equilibrium values
global beta w r L K eta sumK pensshare N expenditures;

Lshare=L/sum(L);
P=(eta*pensshare+(1-eta)*Lshare)*ones(1,N);

[w,K]=fixk(w,sumK);
exp_nohouse= L.*w+ (eta*pensshare+(1-eta)*Lshare)*sum(r*K);
expenditures=inv(eye(N)-(1-beta)*P)*exp_nohouse;