| emplt=0.94, pt=1, piwt=0, pict=0, xt=0 |
| mu=0.05, bpiw=0.5, bpic=1.5, b1=−1.14, |
| b2=5.53, b3=3.68, bx=0.5, ax=1, sx=1.2, |
| sg=4, sp=2, spi=100, g1=0.3, g2=0.8, |
| g3=0.8, g4=0.5, nn=0.05, barp=1 |
| inflt |
| c 0.01 |
| Ut = (1−emplt)100 |
| inflt = (b1+b2/Ut+b3/Ut2)/100+piwt |
| ggt = nn−g1(pt−barp)−g2(inflt−(mu−nn)) |
| + g3(pict−(mu−nn))−g4(piwt−(mu−nn)) |
| emplt1=emplt(ggt−nn) |
| pt1=ptinflt |
| piwt1 = bpiw(((1+xt)/2)inflt−piwt) |
| pict1 = bpic(inflt−pict) |
| xt1 = bx((1−xt)math.exp(ax(sxxt+sg(ggt−nn) |
| − sp(pt−barp)2−spi(pict−(mu−nn))2)) |
| − (1+xt)math.exp(−ax(sxxt+sg(ggt−nn)-sp(pt−barp)2 |
| − spi(pict−(mu−nn))2))) |
| emplt=emplt1 |
| pt=pt1 |
| piwt=piwt1 |
| pict=pict1 |
| xt=xt1 |
