function cap = polfunction(caplag,prodlag,shocks,iii)
global decision capstst

cap =decision(1,iii)+decision(2,iii)*(caplag-capstst) ...	
    +decision(3,iii)*prodlag + decision(4,iii)*shocks;
