clear all tic randn('state', sum(100*rand*clock)); rand('state', sum(100*rand*clock)); run=40; time=100; %number of periods mxfx=23; %maximum number of fixed periods phi=.05; %this is the necessary raise for sticky wages Ai=100; Aj=100; Ak=100; tauij=zeros(1,run); tauk=zeros(1,run); winfav=zeros(1,run); menu=zeros(1,run); mktall=zeros(1,run); mkti=zeros(1,run); mktj=zeros(1,run); mktk=zeros(1,run); bav=zeros(1,run); bvl=zeros(1,run); bi=zeros(1,run); bj=zeros(1,run); bk=zeros(1,run); pvlav=zeros(1,run); pthruavgi=zeros(1,run); pthruavgj=zeros(1,run); pthruavgk=zeros(1,run); pthrurngi=zeros(1,run); pthrurngj=zeros(1,run); pthrurngk=zeros(1,run); pthruij=zeros(1,run); pthruik=zeros(1,run); pthruji=zeros(1,run); pthrujk=zeros(1,run); pthruki=zeros(1,run); pthrukj=zeros(1,run); risk=zeros(1,run); eiklnvar=zeros(1,run); ejklnvar=zeros(1,run); loop=zeros(1,run); wage=zeros(run,time); wgshdw=zeros(run,time); for aaa=1:run tauij(1,aaa)=.01*rand; %transaction cost between country i and j tauk(1,aaa)=rand*.008;%transaction cost with country k is random but could be more than tauij winfav(1,aaa)=(.005/12)+(rand*.080/12); %monthly average inflation is random between .5-8.5prct annually menu(1,aaa)=.01+(rand*.06); %random menu cost between .01 and .07 mkti(1,aaa)=rand; mktj(1,aaa)=rand; mktk(1,aaa)=rand; mktall(1,aaa)=mkti(1,aaa)+mktj(1,aaa)+mktk(1,aaa); bav(1,aaa)=4+(rand*6); %average b for demand (A-bp)=Q bvl(1,aaa)=rand*2; %range of b around the average bi(1,aaa)=bav(1,aaa)+(rand*bvl(1,aaa)); bj(1,aaa)=bav(1,aaa)+(rand*bvl(1,aaa)); bk(1,aaa)=bav(1,aaa)+(rand*bvl(1,aaa)); pvlav(1,aaa)=rand*(.023); %the average var of random movement in underlying price is roughly 1.5%StdDev=.23 or less risk(1,aaa)=.9+(rand*.1); %setting firm risk aversion eiklnvar(1,aaa)=rand*.0005; %setting volatility of two ex rates ejklnvar(1,aaa)=rand*.0005; loop(1,aaa)=.02+(rand*.3); wage(aaa,1)=4+(6*rand); wgshdw(aaa,1)=wage(aaa,1); end leik=log(zeros(run,time)+.5); %starting points for ex rate lejk=log(zeros(run,time)+.5); leij=log(zeros(run,time)+1); eik=zeros(run,time)+.5; %starting points for reg ex. rate ejk=zeros(run,time)+.5; eij=zeros(run,time)+1; infl=zeros(run,time); lpi=log(zeros(run,time)+2); lpj=log(zeros(run,time)+2); lpk=log(zeros(run,time)+4); pi=zeros(run,time)+2; pj=zeros(run,time)+2; pk=zeros(run,time)+4; pj_ink=zeros(time,1)+4; pi_ink=zeros(time,1)+4; pk_ink=zeros(time,1)+4; corrpipjh=zeros(2,2); corrpipkh=zeros(2,2); corrpjpkh=zeros(2,2); corrpipj=zeros(run,1); corrpipk=zeros(run,1); corrpjpk=zeros(run,1); for zrn=1:run pj_ink=zeros(time,1)+4; pi_ink=zeros(time,1)+4; pk_ink=zeros(time,1)+4; corrpipjh=zeros(2,2); corrpipkh=zeros(2,2); corrpjpkh=zeros(2,2); for ztm=2:time leik(zrn,ztm)=leik(zrn,ztm-1)+((eiklnvar(1,zrn)^.5)*randn); %setting exrates and no arbitrage lejk(zrn,ztm)=lejk(zrn,ztm-1)+((ejklnvar(1,zrn)^.5)*randn); leij(zrn,ztm)=leik(zrn,ztm)-lejk(zrn,ztm)-tauij(1,zrn)-(2*tauk(1,zrn)); infl(zrn,ztm)=winfav(1,zrn)+randn*.5*winfav(1,aaa); wgshdw(zrn,ztm)=(1+infl(zrn,ztm))*wgshdw(zrn,ztm-1); wage(zrn,ztm)=wage(zrn,ztm-1); if wgshdw(zrn,ztm)>=wage(zrn,ztm)*(1+phi) wage(zrn,ztm)=wgshdw(zrn,ztm); %setting wage inflation else wage(zrn,ztm)=wage(zrn,ztm-1); end lpi(zrn,ztm)=lpi(zrn,ztm-1)+(pvlav(zrn)*randn); lpj(zrn,ztm)=lpj(zrn,ztm-1)+(pvlav(zrn)*randn); lpk(zrn,ztm)=lpk(zrn,ztm-1)+(pvlav(zrn)*randn); %you may need a no arbitrage contraint here eik(zrn,ztm)=2.718^leik(zrn,ztm); ejk(zrn,ztm)=2.718^lejk(zrn,ztm); eij(zrn,ztm)=2.718^leij(zrn,ztm); pi(zrn,ztm)=2.718^lpi(zrn,ztm); pj(zrn,ztm)=2.718^lpj(zrn,ztm); pk(zrn,ztm)=2.718^lpk(zrn,ztm); if (pj(zrn,ztm))*(1/ejk(zrn,ztm))*(1-loop(zrn))>=(pk(zrn,ztm)) pj(zrn,ztm)=pk(zrn,ztm)*(ejk(zrn,ztm)); lpj(zrn,ztm)=log(pj(zrn,ztm)); end if (pj(zrn,ztm))*(1/ejk(zrn,ztm))*(1+loop(zrn))<=(pk(zrn,ztm)) pj(zrn,ztm)=pk(zrn,ztm)*(ejk(zrn,ztm)); lpj(zrn,ztm)=log(pj(zrn,ztm)); end if (pi(zrn,ztm))*(1/eik(zrn,ztm))*(1-loop(zrn))>=(pk(zrn,ztm)) pi(zrn,ztm)=pk(zrn,ztm)*(eik(zrn,ztm)); lpi(zrn,ztm)=log(pi(zrn,ztm)); end if (pi(zrn,ztm))*(1/eik(zrn,ztm))*(1+loop(zrn))<=(pk(zrn,ztm)) pi(zrn,ztm)=pk(zrn,ztm)*(eik(zrn,ztm)); lpi(zrn,ztm)=log(pi(zrn,ztm)); end pj_ink(ztm,1)=pj(zrn,ztm)*(1/ejk(zrn,ztm)); pi_ink(ztm,1)=pi(zrn,ztm)*(1/eik(zrn,ztm)); pk_ink(ztm,1)=pk(zrn,ztm); end corrpipjh=corrcoef(pi_ink,pj_ink); corrpipkh=corrcoef(pi_ink,pk_ink); corrpjpkh=corrcoef(pj_ink,pk_ink); corrpipj(zrn,1)=corrpipjh(1,2); corrpipk(zrn,1)=corrpipkh(1,2); corrpjpk(zrn,1)=corrpjpkh(1,2); end varlpi=var(lpi'); varlpj=var(lpj'); varlpk=var(lpk'); varleik=var(leik'); varlejk=var(lejk'); varleij=var(leij'); %now to calculate prices and profits pricegv=ones(run*mxfx,time); pricegd=ones(run*mxfx,time); pricedindf=ones(run*mxfx,time); pricevinv=ones(run*mxfx,time); pricedind=ones(run*mxfx,time); pricevinfv=ones(run*mxfx,time); pricefindf=ones(run*mxfx,time); pricedindv=ones(run*mxfx,time); pricefinf=ones(run*mxfx,time); pricevindv=ones(run*mxfx,time); pricevinf3p=ones(run*mxfx,time); pricevind3p=ones(run*mxfx,time); pricedinf3p=ones(run*mxfx,time); pricedinv3p=ones(run*mxfx,time); pricefind3p=ones(run*mxfx,time); pricefinv3p=ones(run*mxfx,time); pricevinv3p=ones(run*mxfx,time); pricedind3p=ones(run*mxfx,time); pricefinf3p=ones(run*mxfx,time); pgvd=0; pgdd=0; pdindfd=0; pvinvd=0; pdindd=0; pvinfvd=0; pfindfd=0; pdindvd=0; pfinfd=0; pvindvd=0; pvinfd=0; pvindd=0; pdinfd=0; pdinvd=0; pfindd=0; pfinvd=0; pgvn=0; pgdn=0; pdindfn=0; pvinvn=0; pdindn=0; pvinfvn=0; pfindfn=0; pdindvn=0; pfinfn=0; pvindvn=0; pvinfn=0; pvindn=0; pdinfn=0; pdinvn=0; pfindn=0; pfinvn=0; %here's a list of the various profit functions profitgv=zeros(run*mxfx,time); profitgd=zeros(run*mxfx,time); profitddv=zeros(run*mxfx,time); profitdvv=zeros(run*mxfx,time); profitffv=zeros(run*mxfx,time); profitdfd=zeros(run*mxfx,time); profitvfv=zeros(run*mxfx,time); profitdfv=zeros(run*mxfx,time); profitgv3p=zeros(run*mxfx,time); profitgd3p=zeros(run*mxfx,time); profitddv3p=zeros(run*mxfx,time); profitdvv3p=zeros(run*mxfx,time); profitffv3p=zeros(run*mxfx,time); profitdfd3p=zeros(run*mxfx,time); profitvfv3p=zeros(run*mxfx,time); lastup=zeros(1,mxfx)+24; for yrn=1:run lastup=zeros(1,mxfx)+24; for ytm=24:time pgvd=0; pgdd=0; pdindfd=0; pvinvd=0; pdindd=0; pvinfvd=0; pfindfd=0; pdindvd=0; pfinfd=0; pvindvd=0; pvinfd=0; pvindd=0; pdinfd=0; pdinvd=0; pfindd=0; pfinvd=0; pgvn=0; pgdn=0; pdindfn=0; pvinvn=0; pdindn=0; pvinfvn=0; pfindfn=0; pdindvn=0; pfinfn=0; pvindvn=0; pvinfn=0; pvindn=0; pdinfn=0; pdinvn=0; pfindn=0; pfinvn=0; for ymx=1:mxfx if ytm>=lastup(ymx)+ymx pricegv((yrn*mxfx)+ymx-mxfx,ytm)=(pgvn+((1-menu(yrn))*((mkti(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ai)+(mktj(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Aj)+(mktk(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ak)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*eik(yrn,ytm)*(1+tauk(yrn))*(1/pi(yrn,ytm)))+(bj(yrn)*mktj(yrn)*ejk(yrn,ytm)*(1+tauk(yrn))*(1/pj(yrn,ytm)))+(bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm)))))))))/(pgvd+((1-menu(yrn))*((((1-tauk(yrn))*eik(yrn,ytm))*((2*bi(yrn)*mkti(yrn)*(1+tauk(yrn))*eik(yrn,ytm)*(1/pi(yrn,ytm)))+(2*bj(yrn)*mktj(yrn)*(1+tauk(yrn))*ejk(yrn,ytm)*(1/pj(yrn,ytm)))+(2*bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm)))))))); pricegd((yrn*mxfx)+ymx-mxfx,ytm)=(pgdn+((1-menu(yrn))*((mkti(yrn)*Ai)+(mktj(yrn)*Aj)+(mktk(yrn)*Ak)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(bj(yrn)*mktj(yrn)*(1/eij(yrn,ytm))*(1+tauij(yrn))*(1/pj(yrn,ytm)))+(bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm)))))))))/(pgdd+((1-menu(yrn))*(((2*bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(2*bj(yrn)*mktj(yrn)*(1+tauij(yrn))*(1/eij(yrn,ytm))*(1/pj(yrn,ytm)))+(2*bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm))))))); pricedindf((yrn*mxfx)+ymx-mxfx,ytm)=(pdindfn+(((1-menu(yrn))^2)*((mkti(yrn)*Ai)+(mktj(yrn)*Aj)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(bj(yrn)*mktj(yrn)*(1/eij(yrn,ytm))*(1+tauij(yrn))*(1/pj(yrn,ytm)))))))))/(pdindfd+(((1-menu(yrn))^2)*(((2*bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(2*bj(yrn)*mktj(yrn)*(1+tauij(yrn))*(1/eij(yrn,ytm))*(1/pj(yrn,ytm))))))); pricevinv((yrn*mxfx)+ymx-mxfx,ytm)=(pvinvn+(((1-menu(yrn))^2)*((mktk(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ak)+((wage(yrn,ytm)*((bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm)))))))))/(pvinvd+(((1-menu(yrn))^2)*(((1-tauk(yrn))*eik(yrn,ytm))*((2*bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm))))))); pricedind((yrn*mxfx)+ymx-mxfx,ytm)=(pdindn+(((1-menu(yrn))^2)*((mkti(yrn)*Ai)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))))))))/(pdindd+(((1-menu(yrn))^2)*(((2*bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm))))))); pricevinfv((yrn*mxfx)+ymx-mxfx,ytm)=(pvinfvn+(((1-menu(yrn))^2)*((mktj(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Aj)+(mktk(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ak)+((wage(yrn,ytm)*((bj(yrn)*mktj(yrn)*ejk(yrn,ytm)*(1+tauk(yrn))*(1/pj(yrn,ytm)))+(bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm)))))))))/(pvinfvd+(((1-menu(yrn))^2)*(((1-tauk(yrn))*eik(yrn,ytm))*((2*bj(yrn)*mktj(yrn)*(1+tauk(yrn))*ejk(yrn,ytm)*(1/pj(yrn,ytm)))+(2*bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm))))))); pricefindf((yrn*mxfx)+ymx-mxfx,ytm)=(pfindfn+(((1-menu(yrn))^2)*((mkti(yrn)*(1-tauij(yrn))*eij(yrn,ytm)*Ai)+(mktj(yrn)*(1-tauij(yrn))*eij(yrn,ytm)*Aj)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*eij(yrn,ytm)*(1+tauij(yrn))*(1/pi(yrn,ytm)))+(bj(yrn)*mktj(yrn)*(1/pj(yrn,ytm)))))))))/(pfindfd+(((1-menu(yrn))^2)*(((1-tauij(yrn))*eij(yrn,ytm))*((2*bi(yrn)*mkti(yrn)*(1+tauij(yrn))*eij(yrn,ytm)*(1/pi(yrn,ytm)))+(2*bj(yrn)*mktj(yrn)*(1/pj(yrn,ytm))))))); pricedindv((yrn*mxfx)+ymx-mxfx,ytm)=(pdindvn+(((1-menu(yrn))^2)*((mkti(yrn)*Ai)+(mktk(yrn)*Ak)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm)))))))))/(pdindvd+(((1-menu(yrn))^2)*(((2*bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(2*bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm))))))); pricefinf((yrn*mxfx)+ymx-mxfx,ytm)=(pfinfn+(((1-menu(yrn))^2)*((mktj(yrn)*(1-tauij(yrn))*eij(yrn,ytm)*Aj)+((wage(yrn,ytm)*((bj(yrn)*mktj(yrn)*(1/pj(yrn,ytm)))))))))/(pfinfd+(((1-menu(yrn))^2)*(((1-tauij(yrn))*eij(yrn,ytm))*((2*bi(yrn)*mkti(yrn)*(1+tauij(yrn))*eij(yrn,ytm)*(1/pi(yrn,ytm))))))); pricevindv((yrn*mxfx)+ymx-mxfx,ytm)=(pdindvn+(((1-menu(yrn))^2)*((mkti(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ai)+(mktk(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ak)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*eik(yrn,ytm)*(1+tauk(yrn))*(1/pi(yrn,ytm)))+(bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm)))))))))/(pdindvd+(((1-menu(yrn))^2)*(((2*bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(2*bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm))))))); pricevinf3p((yrn*mxfx)+ymx-mxfx,ytm)=(pvinfn+(((1-menu(yrn))^3)*((mktj(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Aj)+((wage(yrn,ytm)*((bj(yrn)*mktj(yrn)*ejk(yrn,ytm)*(1+tauk(yrn))*(1/pj(yrn,ytm)))))))))/(pvinfd+(((1-menu(yrn))^3)*(((1-tauk(yrn))*eik(yrn,ytm))*((2*bj(yrn)*mktj(yrn)*(1+tauk(yrn))*ejk(yrn,ytm)*(1/pj(yrn,ytm))))))); pricevind3p((yrn*mxfx)+ymx-mxfx,ytm)=(pvindn+(((1-menu(yrn))^3)*((mkti(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ai)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*eik(yrn,ytm)*(1+tauk(yrn))*(1/pi(yrn,ytm)))))))))/(pvindd+(((1-menu(yrn))^3)*(((1-tauk(yrn))*eik(yrn,ytm))*((2*bi(yrn)*mkti(yrn)*(1+tauk(yrn))*eik(yrn,ytm)*(1/pi(yrn,ytm))))))); pricedinf3p((yrn*mxfx)+ymx-mxfx,ytm)=(pdinfn+(((1-menu(yrn))^3)*((mktj(yrn)*Aj)+((wage(yrn,ytm)*((bj(yrn)*mktj(yrn)*(1/eij(yrn,ytm))*(1+tauij(yrn))*(1/pj(yrn,ytm)))))))))/(pdinfd+(((1-menu(yrn))^3)*(((2*bj(yrn)*mktj(yrn)*(1+tauij(yrn))*(1/eij(yrn,ytm))*(1/pj(yrn,ytm))))))); pricefind3p((yrn*mxfx)+ymx-mxfx,ytm)=(pfindn+(((1-menu(yrn))^3)*((mkti(yrn)*(1-tauij(yrn))*eij(yrn,ytm)*Ai)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*eij(yrn,ytm)*(1+tauij(yrn))*(1/pi(yrn,ytm)))))))))/(pfindd+(((1-menu(yrn))^3)*(((1-tauij(yrn))*eij(yrn,ytm))*((2*bi(yrn)*mkti(yrn)*(1+tauij(yrn))*eij(yrn,ytm)*(1/pi(yrn,ytm))))))); pricevinv3p((yrn*mxfx)+ymx-mxfx,ytm)=(pvinvn+(((1-menu(yrn))^3)*((mktk(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ak)+((wage(yrn,ytm)*((bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm)))))))))/(pvinvd+(((1-menu(yrn))^3)*(((1-tauk(yrn))*eik(yrn,ytm))*((2*bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm))))))); pricedind3p((yrn*mxfx)+ymx-mxfx,ytm)=(pdindn+(((1-menu(yrn))^3)*((mkti(yrn)*Ai)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))))))))/(pdindd+(((1-menu(yrn))^3)*(((2*bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm))))))); pricefinf3p((yrn*mxfx)+ymx-mxfx,ytm)=(pfinfn+(((1-menu(yrn))^3)*((mktj(yrn)*(1-tauij(yrn))*eij(yrn,ytm)*Aj)+((wage(yrn,ytm)*((bj(yrn)*mktj(yrn)*(1/pj(yrn,ytm)))))))))/(pfinfd+(((1-menu(yrn))^3)*(((1-tauij(yrn))*eij(yrn,ytm))*((2*bj(yrn)*mktj(yrn)*(1/pj(yrn,ytm))))))); pricedinv3p((yrn*mxfx)+ymx-mxfx,ytm)=(pdinvn+(((1-menu(yrn))^3)*((mktk(yrn)*Ak)+((wage(yrn,ytm)*((bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm)))))))))/(pdinvd+(((1-menu(yrn))^3)*(((2*bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm))))))); lastup(ymx)=ytm; else pricegv((yrn*mxfx)+ymx-mxfx,ytm)=pricegv((yrn*mxfx)+ymx-mxfx,ytm-1); pricegd((yrn*mxfx)+ymx-mxfx,ytm)=pricegd((yrn*mxfx)+ymx-mxfx,ytm-1); pricedindf((yrn*mxfx)+ymx-mxfx,ytm)=pricedindf((yrn*mxfx)+ymx-mxfx,ytm-1); pricevinv((yrn*mxfx)+ymx-mxfx,ytm)=pricevinv((yrn*mxfx)+ymx-mxfx,ytm-1); pricedind((yrn*mxfx)+ymx-mxfx,ytm)=pricedind((yrn*mxfx)+ymx-mxfx,ytm-1); pricevinfv((yrn*mxfx)+ymx-mxfx,ytm)=pricevinfv((yrn*mxfx)+ymx-mxfx,ytm-1); pricefindf((yrn*mxfx)+ymx-mxfx,ytm)=pricefindf((yrn*mxfx)+ymx-mxfx,ytm-1); pricedindv((yrn*mxfx)+ymx-mxfx,ytm)=pricedindv((yrn*mxfx)+ymx-mxfx,ytm-1); pricefinf((yrn*mxfx)+ymx-mxfx,ytm)=pricefinf((yrn*mxfx)+ymx-mxfx,ytm-1); pricevindv((yrn*mxfx)+ymx-mxfx,ytm)=pricevindv((yrn*mxfx)+ymx-mxfx,ytm-1); pricevinf3p((yrn*mxfx)+ymx-mxfx,ytm)=pricevinf3p((yrn*mxfx)+ymx-mxfx,ytm-1); pricevind3p((yrn*mxfx)+ymx-mxfx,ytm)=pricevind3p((yrn*mxfx)+ymx-mxfx,ytm-1); pricedinf3p((yrn*mxfx)+ymx-mxfx,ytm)=pricedinf3p((yrn*mxfx)+ymx-mxfx,ytm-1); pricedinv3p((yrn*mxfx)+ymx-mxfx,ytm)=pricedinv3p((yrn*mxfx)+ymx-mxfx,ytm-1); pricefind3p((yrn*mxfx)+ymx-mxfx,ytm)=pricefind3p((yrn*mxfx)+ymx-mxfx,ytm-1); pricefinv3p((yrn*mxfx)+ymx-mxfx,ytm)=pricefinv3p((yrn*mxfx)+ymx-mxfx,ytm-1); pricevinv3p((yrn*mxfx)+ymx-mxfx,ytm)=pricevinv3p((yrn*mxfx)+ymx-mxfx,ytm-1); pricedind3p((yrn*mxfx)+ymx-mxfx,ytm)=pricedind3p((yrn*mxfx)+ymx-mxfx,ytm-1); pricefinf3p((yrn*mxfx)+ymx-mxfx,ytm)=pricefinf3p((yrn*mxfx)+ymx-mxfx,ytm-1); end if ymx==1 pricegv((yrn*mxfx)+ymx-mxfx,ytm)=(pgvn+((1-menu(yrn))*((mkti(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ai)+(mktj(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Aj)+(mktk(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ak)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*eik(yrn,ytm)*(1+tauk(yrn))*(1/pi(yrn,ytm)))+(bj(yrn)*mktj(yrn)*ejk(yrn,ytm)*(1+tauk(yrn))*(1/pj(yrn,ytm)))+(bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm)))))))))/(pgvd+((1-menu(yrn))*((((1-tauk(yrn))*eik(yrn,ytm))*((2*bi(yrn)*mkti(yrn)*(1+tauk(yrn))*eik(yrn,ytm)*(1/pi(yrn,ytm)))+(2*bj(yrn)*mktj(yrn)*(1+tauk(yrn))*ejk(yrn,ytm)*(1/pj(yrn,ytm)))+(2*bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm)))))))); pricegd((yrn*mxfx)+ymx-mxfx,ytm)=(pgdn+((1-menu(yrn))*((mkti(yrn)*Ai)+(mktj(yrn)*Aj)+(mktk(yrn)*Ak)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(bj(yrn)*mktj(yrn)*(1/eij(yrn,ytm))*(1+tauij(yrn))*(1/pj(yrn,ytm)))+(bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm)))))))))/(pgdd+((1-menu(yrn))*(((2*bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(2*bj(yrn)*mktj(yrn)*(1+tauij(yrn))*(1/eij(yrn,ytm))*(1/pj(yrn,ytm)))+(2*bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm))))))); pricedindf((yrn*mxfx)+ymx-mxfx,ytm)=(pdindfn+(((1-menu(yrn))^2)*((mkti(yrn)*Ai)+(mktj(yrn)*Aj)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(bj(yrn)*mktj(yrn)*(1/eij(yrn,ytm))*(1+tauij(yrn))*(1/pj(yrn,ytm)))))))))/(pdindfd+(((1-menu(yrn))^2)*(((2*bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(2*bj(yrn)*mktj(yrn)*(1+tauij(yrn))*(1/eij(yrn,ytm))*(1/pj(yrn,ytm))))))); pricevinv((yrn*mxfx)+ymx-mxfx,ytm)=(pvinvn+(((1-menu(yrn))^2)*((mktk(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ak)+((wage(yrn,ytm)*((bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm)))))))))/(pvinvd+(((1-menu(yrn))^2)*(((1-tauk(yrn))*eik(yrn,ytm))*((2*bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm))))))); pricedind((yrn*mxfx)+ymx-mxfx,ytm)=(pdindn+(((1-menu(yrn))^2)*((mkti(yrn)*Ai)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))))))))/(pdindd+(((1-menu(yrn))^2)*(((2*bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm))))))); pricevinfv((yrn*mxfx)+ymx-mxfx,ytm)=(pvinfvn+(((1-menu(yrn))^2)*((mktj(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Aj)+(mktk(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ak)+((wage(yrn,ytm)*((bj(yrn)*mktj(yrn)*ejk(yrn,ytm)*(1+tauk(yrn))*(1/pj(yrn,ytm)))+(bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm)))))))))/(pvinfvd+(((1-menu(yrn))^2)*(((1-tauk(yrn))*eik(yrn,ytm))*((2*bj(yrn)*mktj(yrn)*(1+tauk(yrn))*ejk(yrn,ytm)*(1/pj(yrn,ytm)))+(2*bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm))))))); pricefindf((yrn*mxfx)+ymx-mxfx,ytm)=(pfindfn+(((1-menu(yrn))^2)*((mkti(yrn)*(1-tauij(yrn))*eij(yrn,ytm)*Ai)+(mktj(yrn)*(1-tauij(yrn))*eij(yrn,ytm)*Aj)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*eij(yrn,ytm)*(1+tauij(yrn))*(1/pi(yrn,ytm)))+(bj(yrn)*mktj(yrn)*(1/pj(yrn,ytm)))))))))/(pfindfd+(((1-menu(yrn))^2)*(((1-tauij(yrn))*eij(yrn,ytm))*((2*bi(yrn)*mkti(yrn)*(1+tauij(yrn))*eij(yrn,ytm)*(1/pi(yrn,ytm)))+(2*bj(yrn)*mktj(yrn)*(1/pj(yrn,ytm))))))); pricedindv((yrn*mxfx)+ymx-mxfx,ytm)=(pdindvn+(((1-menu(yrn))^2)*((mkti(yrn)*Ai)+(mktk(yrn)*Ak)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm)))))))))/(pdindvd+(((1-menu(yrn))^2)*(((2*bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(2*bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm))))))); pricefinf((yrn*mxfx)+ymx-mxfx,ytm)=(pfinfn+(((1-menu(yrn))^2)*((mktj(yrn)*(1-tauij(yrn))*eij(yrn,ytm)*Aj)+((wage(yrn,ytm)*((bj(yrn)*mktj(yrn)*(1/pj(yrn,ytm)))))))))/(pfinfd+(((1-menu(yrn))^2)*(((1-tauij(yrn))*eij(yrn,ytm))*((2*bi(yrn)*mkti(yrn)*(1+tauij(yrn))*eij(yrn,ytm)*(1/pi(yrn,ytm))))))); pricevindv((yrn*mxfx)+ymx-mxfx,ytm)=(pdindvn+(((1-menu(yrn))^2)*((mkti(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ai)+(mktk(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ak)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*eik(yrn,ytm)*(1+tauk(yrn))*(1/pi(yrn,ytm)))+(bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm)))))))))/(pdindvd+(((1-menu(yrn))^2)*(((2*bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(2*bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm))))))); pricevinf3p((yrn*mxfx)+ymx-mxfx,ytm)=(pvinfn+(((1-menu(yrn))^3)*((mktj(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Aj)+((wage(yrn,ytm)*((bj(yrn)*mktj(yrn)*ejk(yrn,ytm)*(1+tauk(yrn))*(1/pj(yrn,ytm)))))))))/(pvinfd+(((1-menu(yrn))^3)*(((1-tauk(yrn))*eik(yrn,ytm))*((2*bj(yrn)*mktj(yrn)*(1+tauk(yrn))*ejk(yrn,ytm)*(1/pj(yrn,ytm))))))); pricevind3p((yrn*mxfx)+ymx-mxfx,ytm)=(pvindn+(((1-menu(yrn))^3)*((mkti(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ai)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*eik(yrn,ytm)*(1+tauk(yrn))*(1/pi(yrn,ytm)))))))))/(pvindd+(((1-menu(yrn))^3)*(((1-tauk(yrn))*eik(yrn,ytm))*((2*bi(yrn)*mkti(yrn)*(1+tauk(yrn))*eik(yrn,ytm)*(1/pi(yrn,ytm))))))); pricedinf3p((yrn*mxfx)+ymx-mxfx,ytm)=(pdinfn+(((1-menu(yrn))^3)*((mktj(yrn)*Aj)+((wage(yrn,ytm)*((bj(yrn)*mktj(yrn)*(1/eij(yrn,ytm))*(1+tauij(yrn))*(1/pj(yrn,ytm)))))))))/(pdinfd+(((1-menu(yrn))^3)*(((2*bj(yrn)*mktj(yrn)*(1+tauij(yrn))*(1/eij(yrn,ytm))*(1/pj(yrn,ytm))))))); pricefind3p((yrn*mxfx)+ymx-mxfx,ytm)=(pfindn+(((1-menu(yrn))^3)*((mkti(yrn)*(1-tauij(yrn))*eij(yrn,ytm)*Ai)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*eij(yrn,ytm)*(1+tauij(yrn))*(1/pi(yrn,ytm)))))))))/(pfindd+(((1-menu(yrn))^3)*(((1-tauij(yrn))*eij(yrn,ytm))*((2*bi(yrn)*mkti(yrn)*(1+tauij(yrn))*eij(yrn,ytm)*(1/pi(yrn,ytm))))))); pricevinv3p((yrn*mxfx)+ymx-mxfx,ytm)=(pvinvn+(((1-menu(yrn))^3)*((mktk(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ak)+((wage(yrn,ytm)*((bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm)))))))))/(pvinvd+(((1-menu(yrn))^3)*(((1-tauk(yrn))*eik(yrn,ytm))*((2*bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm))))))); pricedind3p((yrn*mxfx)+ymx-mxfx,ytm)=(pdindn+(((1-menu(yrn))^3)*((mkti(yrn)*Ai)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))))))))/(pdindd+(((1-menu(yrn))^3)*(((2*bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm))))))); pricefinf3p((yrn*mxfx)+ymx-mxfx,ytm)=(pfinfn+(((1-menu(yrn))^3)*((mktj(yrn)*(1-tauij(yrn))*eij(yrn,ytm)*Aj)+((wage(yrn,ytm)*((bj(yrn)*mktj(yrn)*(1/pj(yrn,ytm)))))))))/(pfinfd+(((1-menu(yrn))^3)*(((1-tauij(yrn))*eij(yrn,ytm))*((2*bj(yrn)*mktj(yrn)*(1/pj(yrn,ytm))))))); pricedinv3p((yrn*mxfx)+ymx-mxfx,ytm)=(pdinvn+(((1-menu(yrn))^3)*((mktk(yrn)*Ak)+((wage(yrn,ytm)*((bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm)))))))))/(pdinvd+(((1-menu(yrn))^3)*(((2*bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm))))))); lastup(ymx)=ytm; end if ytm==24 pricegv((yrn*mxfx)+ymx-mxfx,ytm)=(pgvn+((1-menu(yrn))*((mkti(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ai)+(mktj(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Aj)+(mktk(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ak)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*eik(yrn,ytm)*(1+tauk(yrn))*(1/pi(yrn,ytm)))+(bj(yrn)*mktj(yrn)*ejk(yrn,ytm)*(1+tauk(yrn))*(1/pj(yrn,ytm)))+(bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm)))))))))/(pgvd+((1-menu(yrn))*((((1-tauk(yrn))*eik(yrn,ytm))*((2*bi(yrn)*mkti(yrn)*(1+tauk(yrn))*eik(yrn,ytm)*(1/pi(yrn,ytm)))+(2*bj(yrn)*mktj(yrn)*(1+tauk(yrn))*ejk(yrn,ytm)*(1/pj(yrn,ytm)))+(2*bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm)))))))); pricegd((yrn*mxfx)+ymx-mxfx,ytm)=(pgdn+((1-menu(yrn))*((mkti(yrn)*Ai)+(mktj(yrn)*Aj)+(mktk(yrn)*Ak)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(bj(yrn)*mktj(yrn)*(1/eij(yrn,ytm))*(1+tauij(yrn))*(1/pj(yrn,ytm)))+(bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm)))))))))/(pgdd+((1-menu(yrn))*(((2*bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(2*bj(yrn)*mktj(yrn)*(1+tauij(yrn))*(1/eij(yrn,ytm))*(1/pj(yrn,ytm)))+(2*bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm))))))); pricedindf((yrn*mxfx)+ymx-mxfx,ytm)=(pdindfn+(((1-menu(yrn))^2)*((mkti(yrn)*Ai)+(mktj(yrn)*Aj)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(bj(yrn)*mktj(yrn)*(1/eij(yrn,ytm))*(1+tauij(yrn))*(1/pj(yrn,ytm)))))))))/(pdindfd+(((1-menu(yrn))^2)*(((2*bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(2*bj(yrn)*mktj(yrn)*(1+tauij(yrn))*(1/eij(yrn,ytm))*(1/pj(yrn,ytm))))))); pricevinv((yrn*mxfx)+ymx-mxfx,ytm)=(pvinvn+(((1-menu(yrn))^2)*((mktk(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ak)+((wage(yrn,ytm)*((bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm)))))))))/(pvinvd+(((1-menu(yrn))^2)*(((1-tauk(yrn))*eik(yrn,ytm))*((2*bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm))))))); pricedind((yrn*mxfx)+ymx-mxfx,ytm)=(pdindn+(((1-menu(yrn))^2)*((mkti(yrn)*Ai)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))))))))/(pdindd+(((1-menu(yrn))^2)*(((2*bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm))))))); pricevinfv((yrn*mxfx)+ymx-mxfx,ytm)=(pvinfvn+(((1-menu(yrn))^2)*((mktj(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Aj)+(mktk(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ak)+((wage(yrn,ytm)*((bj(yrn)*mktj(yrn)*ejk(yrn,ytm)*(1+tauk(yrn))*(1/pj(yrn,ytm)))+(bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm)))))))))/(pvinfvd+(((1-menu(yrn))^2)*(((1-tauk(yrn))*eik(yrn,ytm))*((2*bj(yrn)*mktj(yrn)*(1+tauk(yrn))*ejk(yrn,ytm)*(1/pj(yrn,ytm)))+(2*bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm))))))); pricefindf((yrn*mxfx)+ymx-mxfx,ytm)=(pfindfn+(((1-menu(yrn))^2)*((mkti(yrn)*(1-tauij(yrn))*eij(yrn,ytm)*Ai)+(mktj(yrn)*(1-tauij(yrn))*eij(yrn,ytm)*Aj)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*eij(yrn,ytm)*(1+tauij(yrn))*(1/pi(yrn,ytm)))+(bj(yrn)*mktj(yrn)*(1/pj(yrn,ytm)))))))))/(pfindfd+(((1-menu(yrn))^2)*(((1-tauij(yrn))*eij(yrn,ytm))*((2*bi(yrn)*mkti(yrn)*(1+tauij(yrn))*eij(yrn,ytm)*(1/pi(yrn,ytm)))+(2*bj(yrn)*mktj(yrn)*(1/pj(yrn,ytm))))))); pricedindv((yrn*mxfx)+ymx-mxfx,ytm)=(pdindvn+(((1-menu(yrn))^2)*((mkti(yrn)*Ai)+(mktk(yrn)*Ak)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm)))))))))/(pdindvd+(((1-menu(yrn))^2)*(((2*bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(2*bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm))))))); pricefinf((yrn*mxfx)+ymx-mxfx,ytm)=(pfinfn+(((1-menu(yrn))^2)*((mktj(yrn)*(1-tauij(yrn))*eij(yrn,ytm)*Aj)+((wage(yrn,ytm)*((bj(yrn)*mktj(yrn)*(1/pj(yrn,ytm)))))))))/(pfinfd+(((1-menu(yrn))^2)*(((1-tauij(yrn))*eij(yrn,ytm))*((2*bi(yrn)*mkti(yrn)*(1+tauij(yrn))*eij(yrn,ytm)*(1/pi(yrn,ytm))))))); pricevindv((yrn*mxfx)+ymx-mxfx,ytm)=(pdindvn+(((1-menu(yrn))^2)*((mkti(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ai)+(mktk(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ak)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*eik(yrn,ytm)*(1+tauk(yrn))*(1/pi(yrn,ytm)))+(bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm)))))))))/(pdindvd+(((1-menu(yrn))^2)*(((2*bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(2*bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm))))))); pricevinf3p((yrn*mxfx)+ymx-mxfx,ytm)=(pvinfn+(((1-menu(yrn))^3)*((mktj(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Aj)+((wage(yrn,ytm)*((bj(yrn)*mktj(yrn)*ejk(yrn,ytm)*(1+tauk(yrn))*(1/pj(yrn,ytm)))))))))/(pvinfd+(((1-menu(yrn))^3)*(((1-tauk(yrn))*eik(yrn,ytm))*((2*bj(yrn)*mktj(yrn)*(1+tauk(yrn))*ejk(yrn,ytm)*(1/pj(yrn,ytm))))))); pricevind3p((yrn*mxfx)+ymx-mxfx,ytm)=(pvindn+(((1-menu(yrn))^3)*((mkti(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ai)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*eik(yrn,ytm)*(1+tauk(yrn))*(1/pi(yrn,ytm)))))))))/(pvindd+(((1-menu(yrn))^3)*(((1-tauk(yrn))*eik(yrn,ytm))*((2*bi(yrn)*mkti(yrn)*(1+tauk(yrn))*eik(yrn,ytm)*(1/pi(yrn,ytm))))))); pricedinf3p((yrn*mxfx)+ymx-mxfx,ytm)=(pdinfn+(((1-menu(yrn))^3)*((mktj(yrn)*Aj)+((wage(yrn,ytm)*((bj(yrn)*mktj(yrn)*(1/eij(yrn,ytm))*(1+tauij(yrn))*(1/pj(yrn,ytm)))))))))/(pdinfd+(((1-menu(yrn))^3)*(((2*bj(yrn)*mktj(yrn)*(1+tauij(yrn))*(1/eij(yrn,ytm))*(1/pj(yrn,ytm))))))); pricefind3p((yrn*mxfx)+ymx-mxfx,ytm)=(pfindn+(((1-menu(yrn))^3)*((mkti(yrn)*(1-tauij(yrn))*eij(yrn,ytm)*Ai)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*eij(yrn,ytm)*(1+tauij(yrn))*(1/pi(yrn,ytm)))))))))/(pfindd+(((1-menu(yrn))^3)*(((1-tauij(yrn))*eij(yrn,ytm))*((2*bi(yrn)*mkti(yrn)*(1+tauij(yrn))*eij(yrn,ytm)*(1/pi(yrn,ytm))))))); pricevinv3p((yrn*mxfx)+ymx-mxfx,ytm)=(pvinvn+(((1-menu(yrn))^3)*((mktk(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ak)+((wage(yrn,ytm)*((bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm)))))))))/(pvinvd+(((1-menu(yrn))^3)*(((1-tauk(yrn))*eik(yrn,ytm))*((2*bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm))))))); pricedind3p((yrn*mxfx)+ymx-mxfx,ytm)=(pdindn+(((1-menu(yrn))^3)*((mkti(yrn)*Ai)+((wage(yrn,ytm)*((bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))))))))/(pdindd+(((1-menu(yrn))^3)*(((2*bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm))))))); pricefinf3p((yrn*mxfx)+ymx-mxfx,ytm)=(pfinfn+(((1-menu(yrn))^3)*((mktj(yrn)*(1-tauij(yrn))*eij(yrn,ytm)*Aj)+((wage(yrn,ytm)*((bj(yrn)*mktj(yrn)*(1/pj(yrn,ytm)))))))))/(pfinfd+(((1-menu(yrn))^3)*(((1-tauij(yrn))*eij(yrn,ytm))*((2*bj(yrn)*mktj(yrn)*(1/pj(yrn,ytm))))))); pricedinv3p((yrn*mxfx)+ymx-mxfx,ytm)=(pdinvn+(((1-menu(yrn))^3)*((mktk(yrn)*Ak)+((wage(yrn,ytm)*((bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm)))))))))/(pdinvd+(((1-menu(yrn))^3)*(((2*bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm))))))); lastup(ymx)=ytm; end pgvn=pgvn+((mkti(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ai)+(mktj(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Aj)+(mktk(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ak)+((wgshdw(yrn,ytm)*((1+winfav(yrn))^ymx)*((bi(yrn)*mkti(yrn)*eik(yrn,ytm)*(1+tauk(yrn))*(1/pi(yrn,ytm)))+(bj(yrn)*mktj(yrn)*ejk(yrn,ytm)*(1+tauk(yrn))*(1/pj(yrn,ytm)))+(bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm))))))); pgvd=pgvd+(((1-tauk(yrn))*eik(yrn,ytm))*((2*bi(yrn)*mkti(yrn)*(1+tauk(yrn))*eik(yrn,ytm)*(1/pi(yrn,ytm)))+(2*bj(yrn)*mktj(yrn)*(1+tauk(yrn))*ejk(yrn,ytm)*(1/pj(yrn,ytm)))+(2*bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm))))); pgdn=pgdn+((mkti(yrn)*Ai)+(mktj(yrn)*Aj)+(mktk(yrn)*Ak)+((wgshdw(yrn,ytm)*((1+winfav(yrn))^ymx)*((bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(bj(yrn)*mktj(yrn)*(1/eij(yrn,ytm))*(1+tauij(yrn))*(1/pj(yrn,ytm)))+(bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm))))))); pgdd=pgdd+(((2*bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(2*bj(yrn)*mktj(yrn)*(1+tauij(yrn))*(1/eij(yrn,ytm))*(1/pj(yrn,ytm)))+(2*bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm))))); pdindfn=pdindfn+((mkti(yrn)*Ai)+(mktj(yrn)*Aj)+((wgshdw(yrn,ytm)*((1+winfav(yrn))^ymx)*((bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(bj(yrn)*mktj(yrn)*(1/eij(yrn,ytm))*(1+tauij(yrn))*(1/pj(yrn,ytm))))))); pdindfd=pdindfd+(((2*bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(2*bj(yrn)*mktj(yrn)*(1+tauij(yrn))*(1/eij(yrn,ytm))*(1/pj(yrn,ytm))))); pvinvn=pvinvn+((mktk(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ak)+((wgshdw(yrn,ytm)*((1+winfav(yrn))^ymx)*((bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm))))))); pvinvd=pvinvd+(((1-tauk(yrn))*eik(yrn,ytm))*((2*bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm))))); pdindn=pdindn+((mkti(yrn)*Ai)+((wgshdw(yrn,ytm)*((1+winfav(yrn))^ymx)*((bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm))))))); pdindd=pdindd+(((2*bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm))))); pvinfvn=pvinfvn+((mktj(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Aj)+(mktk(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ak)+((wgshdw(yrn,ytm)*((1+winfav(yrn))^ymx)*((bj(yrn)*mktj(yrn)*ejk(yrn,ytm)*(1+tauk(yrn))*(1/pj(yrn,ytm)))+(bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm))))))); pvinfvd=pvinfvd+(((1-tauk(yrn))*eik(yrn,ytm))*((2*bj(yrn)*mktj(yrn)*(1+tauk(yrn))*ejk(yrn,ytm)*(1/pj(yrn,ytm)))+(2*bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm))))); pdindvn=pdindvn+((mkti(yrn)*Ai)+(mktk(yrn)*Ak)+((wgshdw(yrn,ytm)*((1+winfav(yrn))^ymx)*((bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm))))))); pdindvd=pdindvd+(((2*bi(yrn)*mkti(yrn)*(1/pi(yrn,ytm)))+(2*bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm))))); pvindvn=pvindvn+((mkti(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ai)+(mktk(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ak)+((wgshdw(yrn,ytm)*((1+winfav(yrn))^ymx)*((bi(yrn)*mkti(yrn)*eik(yrn,ytm)*(1+tauk(yrn))*(1/pi(yrn,ytm)))+(bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm))))))); pvindvd=pvindvd+(((1-tauk(yrn))*eik(yrn,ytm))*((2*bi(yrn)*mkti(yrn)*(1+tauk(yrn))*eik(yrn,ytm)*(1/pi(yrn,ytm)))+(2*bk(yrn)*mktk(yrn)*(1/pk(yrn,ytm))))); pfinfn=pfinfn+((mktj(yrn)*(1-tauij(yrn))*eij(yrn,ytm)*Aj)+((wgshdw(yrn,ytm)*((1+winfav(yrn))^ymx)*((bj(yrn)*mktj(yrn)*(1/pj(yrn,ytm))))))); pfinfd=pfinfd+(((1-tauij(yrn))*eij(yrn,ytm))*((2*bj(yrn)*mktj(yrn)*(1/pj(yrn,ytm))))); pfindfn=pfindfn+((mkti(yrn)*(1-tauij(yrn))*eij(yrn,ytm)*Ai)+(mktj(yrn)*(1-tauij(yrn))*eij(yrn,ytm)*Aj)+((wgshdw(yrn,ytm)*((1+winfav(yrn))^ymx)*((bi(yrn)*mkti(yrn)*eij(yrn,ytm)*(1+tauij(yrn))*(1/pi(yrn,ytm)))+(bj(yrn)*mktj(yrn)*(1/pj(yrn,ytm))))))); pfindfd=pfindfd+(((1-tauij(yrn))*eij(yrn,ytm))*((2*bi(yrn)*mkti(yrn)*(1+tauij(yrn))*eij(yrn,ytm)*(1/pi(yrn,ytm)))+(2*bj(yrn)*mktj(yrn)*(1/pj(yrn,ytm))))); pvinfn=pvinfn+((mktj(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Aj)+((wgshdw(yrn,ytm)*((1+winfav(yrn))^ymx)*((bj(yrn)*mktj(yrn)*ejk(yrn,ytm)*(1+tauk(yrn))*(1/pj(yrn,ytm))))))); pvinfd=pvinfd+(((1-tauk(yrn))*eik(yrn,ytm))*((2*bj(yrn)*mktj(yrn)*(1+tauk(yrn))*ejk(yrn,ytm)*(1/pj(yrn,ytm))))); pvindn=pvindn+((mkti(yrn)*(1-tauk(yrn))*eik(yrn,ytm)*Ai)+((wgshdw(yrn,ytm)*((1+winfav(yrn))^ymx)*((bi(yrn)*mkti(yrn)*eik(yrn,ytm)*(1+tauk(yrn))*(1/pi(yrn,ytm))))))); pvindd=pvindd+(((1-tauk(yrn))*eik(yrn,ytm))*((2*bi(yrn)*mkti(yrn)*(1+tauk(yrn))*eik(yrn,ytm)*(1/pi(yrn,ytm))))); pdinfn=pdinfn+((mktj(yrn)*Aj)+((wgshdw(yrn,ytm)*((1+winfav(yrn))^ymx)*((bj(yrn)*mktj(yrn)*(1/eij(yrn,ytm))*(1+tauij(yrn))*(1/pj(yrn,ytm))))))); pdinfd=pdinfd+(((2*bj(yrn)*mktj(yrn)*(1+tauij(yrn))*(1/eij(yrn,ytm))*(1/pj(yrn,ytm))))); pdinvn=pdinvn+((mktk(yrn)*Ak)+((wgshdw(yrn,ytm)*((1+winfav(yrn))^ymx)*((bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm))))))); pdinvd=pdinvd+(((2*bk(yrn)*mktk(yrn)*(1+tauk(yrn))*(1/eik(yrn,ytm))*(1/pk(yrn,ytm))))); pfindn=pfindn+((mkti(yrn)*(1-tauij(yrn))*eij(yrn,ytm)*Ai)+((wgshdw(yrn,ytm)*((1+winfav(yrn))^ymx)*((bi(yrn)*mkti(yrn)*eij(yrn,ytm)*(1+tauij(yrn))*(1/pi(yrn,ytm))))))); pfindd=pfindd+(((1-tauij(yrn))*eij(yrn,ytm))*((2*bi(yrn)*mkti(yrn)*(1+tauij(yrn))*eij(yrn,ytm)*(1/pi(yrn,ytm))))); end end end for wrn=1:run for wmx=1:mxfx for wtm=48:time-1 if pricegv((wrn*mxfx)+wmx-mxfx,wtm+1)==pricegv((wrn*mxfx)+wmx-mxfx,wtm) profitgv((wrn*mxfx)+wmx-mxfx,wtm)=(((pricegv((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(1,wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ai-(bi(1,wrn)*(pricegv((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauk(1,wrn))*(eik(wrn,wtm))/(pi(wrn,wtm)))*mkti(1,wrn)))+((pricegv((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(1,wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Aj-(bj(1,wrn)*(pricegv((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauk(1,wrn))*(ejk(wrn,wtm))/(pj(wrn,wtm)))*mktj(1,wrn)))+((pricegv((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(1,wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ak-(bk(1,wrn)*(pricegv((wrn*mxfx)+wmx-mxfx,wtm)/(pk(wrn,wtm)))*mktk(1,wrn))))^risk(1,wrn); profitgd((wrn*mxfx)+wmx-mxfx,wtm)=(((pricegd((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricegd((wrn*mxfx)+wmx-mxfx,wtm)/(pi(wrn,wtm)))*mkti(wrn)))+((pricegd((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricegd((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauij(wrn))*((1)/(eij(wrn,wtm)))/(pj(wrn,wtm)))*mktj(wrn)))+((pricegd((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricegd((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauk(wrn))*((1)/(eik(wrn,wtm)))/(pk(wrn,wtm)))*mktk(wrn))))^risk(wrn); profitddv((wrn*mxfx)+wmx-mxfx,wtm)=(((pricedindf((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricedindf((wrn*mxfx)+wmx-mxfx,wtm)/(pi(wrn,wtm)))*mkti(wrn)))+((pricedindf((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricedindf((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauij(wrn))*((1)/(eij(wrn,wtm)))/(pj(wrn,wtm)))*mktj(wrn)))+((pricevinv((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricevinv((wrn*mxfx)+wmx-mxfx,wtm)/(pk(wrn,wtm)))*mktk(wrn))))^risk(wrn); profitdvv((wrn*mxfx)+wmx-mxfx,wtm)=(((pricedind((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricedind((wrn*mxfx)+wmx-mxfx,wtm)/(pi(wrn,wtm)))*mkti(wrn)))+((pricevinfv((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricevinfv((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauk(wrn))*((1)*(ejk(wrn,wtm)))/(pj(wrn,wtm)))*mktj(wrn)))+((pricevinfv((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricevinfv((wrn*mxfx)+wmx-mxfx,wtm)/(pk(wrn,wtm)))*mktk(wrn))))^risk(wrn); profitffv((wrn*mxfx)+wmx-mxfx,wtm)=(((pricefindf((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauij(wrn))*(eij(wrn,wtm-1)))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricefindf((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauij(wrn))*((1)*(eij(wrn,wtm)))/(pi(wrn,wtm)))*mkti(wrn)))+((pricefindf((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauij(wrn))*(eij(wrn,wtm-1)))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricefindf((wrn*mxfx)+wmx-mxfx,wtm)/(pj(wrn,wtm)))*mktj(wrn)))+((pricevinv((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricevinv((wrn*mxfx)+wmx-mxfx,wtm)/(pk(wrn,wtm)))*mktk(wrn))))^risk(wrn); profitdfd((wrn*mxfx)+wmx-mxfx,wtm)=(((pricedindv((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricedindv((wrn*mxfx)+wmx-mxfx,wtm)/(pi(wrn,wtm)))*mkti(wrn)))+((pricefinf((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauij(wrn))*(eij(wrn,wtm-1)))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricefinf((wrn*mxfx)+wmx-mxfx,wtm)/(pj(wrn,wtm)))*mktj(wrn)))+((pricedindv((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricedindv((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauk(wrn))*((1)/(eik(wrn,wtm)))/(pk(wrn,wtm)))*mktk(wrn))))^risk(wrn); profitvfv((wrn*mxfx)+wmx-mxfx,wtm)=(((pricevindv((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricevindv((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauk(wrn))*((1)*(eik(wrn,wtm)))/(pi(wrn,wtm)))*mkti(wrn)))+((pricefinf((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauij(wrn))*(eij(wrn,wtm-1)))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricefinf((wrn*mxfx)+wmx-mxfx,wtm)/(pj(wrn,wtm)))*mktj(wrn)))+((pricevindv((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricevindv((wrn*mxfx)+wmx-mxfx,wtm)/(pk(wrn,wtm)))*mktk(wrn))))^risk(wrn); profitdfv((wrn*mxfx)+wmx-mxfx,wtm)=(((pricedind((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricedind((wrn*mxfx)+wmx-mxfx,wtm)/(pi(wrn,wtm)))*mkti(wrn)))+((pricefinf((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauij(wrn))*(eij(wrn,wtm-1)))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricefinf((wrn*mxfx)+wmx-mxfx,wtm)/(pj(wrn,wtm)))*mktj(wrn)))+((pricevinv((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricevinv((wrn*mxfx)+wmx-mxfx,wtm)/(pk(wrn,wtm)))*mktk(wrn))))^risk(wrn); profitgv3p((wrn*mxfx)+wmx-mxfx,wtm)=(((pricevind3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricevind3p((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauk(wrn))*(eik(wrn,wtm))/(pi(wrn,wtm)))*mkti(wrn)))+((pricevinf3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricevinf3p((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauk(wrn))*(ejk(wrn,wtm))/(pj(wrn,wtm)))*mktj(wrn)))+((pricevinv3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricevinv3p((wrn*mxfx)+wmx-mxfx,wtm)/(pk(wrn,wtm)))*mktk(wrn))))^risk(wrn); profitgd3p((wrn*mxfx)+wmx-mxfx,wtm)=(((pricedind3p((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricedind3p((wrn*mxfx)+wmx-mxfx,wtm)/(pi(wrn,wtm)))*mkti(wrn)))+((pricedinf3p((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricedinf3p((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauij(wrn))*((1)/(eij(wrn,wtm)))/(pj(wrn,wtm)))*mktj(wrn)))+((pricedinv3p((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricedinv3p((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauk(wrn))*((1)/(eik(wrn,wtm)))/(pk(wrn,wtm)))*mktk(wrn))))^risk(wrn); profitddv3p((wrn*mxfx)+wmx-mxfx,wtm)=(((pricedind3p((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricedind3p((wrn*mxfx)+wmx-mxfx,wtm)/(pi(wrn,wtm)))*mkti(wrn)))+((pricedinf3p((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricedinf3p((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauij(wrn))*((1)/(eij(wrn,wtm)))/(pj(wrn,wtm)))*mktj(wrn)))+((pricevinv3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricevinv3p((wrn*mxfx)+wmx-mxfx,wtm)/(pk(wrn,wtm)))*mktk(wrn))))^risk(wrn); profitdvv3p((wrn*mxfx)+wmx-mxfx,wtm)=(((pricedind3p((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricedind3p((wrn*mxfx)+wmx-mxfx,wtm)/(pi(wrn,wtm)))*mkti(wrn)))+((pricevinf3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricevinf3p((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauk(wrn))*((1)*(ejk(wrn,wtm)))/(pj(wrn,wtm)))*mktj(wrn)))+((pricevinv3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricevinv3p((wrn*mxfx)+wmx-mxfx,wtm)/(pk(wrn,wtm)))*mktk(wrn))))^risk(wrn); profitffv3p((wrn*mxfx)+wmx-mxfx,wtm)=(((pricefind3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauij(wrn))*(eij(wrn,wtm-1)))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricefind3p((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauij(wrn))*((1)*(eij(wrn,wtm)))/(pi(wrn,wtm)))*mkti(wrn)))+((pricefinf3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauij(wrn))*(eij(wrn,wtm-1)))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricefinf3p((wrn*mxfx)+wmx-mxfx,wtm)/(pj(wrn,wtm)))*mktj(wrn)))+((pricevinv3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricevinv3p((wrn*mxfx)+wmx-mxfx,wtm)/(pk(wrn,wtm)))*mktk(wrn))))^risk(wrn); profitdfd3p((wrn*mxfx)+wmx-mxfx,wtm)=(((pricedind3p((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricedind3p((wrn*mxfx)+wmx-mxfx,wtm)/(pi(wrn,wtm)))*mkti(wrn)))+((pricefinf3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauij(wrn))*(eij(wrn,wtm-1)))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricefinf3p((wrn*mxfx)+wmx-mxfx,wtm)/(pj(wrn,wtm)))*mktj(wrn)))+((pricedinv3p((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricedinv3p((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauk(wrn))*((1)/(eik(wrn,wtm)))/(pk(wrn,wtm)))*mktk(wrn))))^risk(wrn); profitvfv3p((wrn*mxfx)+wmx-mxfx,wtm)=(((pricevind3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricevind3p((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauk(wrn))*((1)*(eik(wrn,wtm)))/(pi(wrn,wtm)))*mkti(wrn)))+((pricefinf3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauij(wrn))*(eij(wrn,wtm-1)))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricefinf3p((wrn*mxfx)+wmx-mxfx,wtm)/(pj(wrn,wtm)))*mktj(wrn)))+((pricevinv3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricevinv3p((wrn*mxfx)+wmx-mxfx,wtm)/(pk(wrn,wtm)))*mktk(wrn))))^risk(wrn); else profitgv((wrn*mxfx)+wmx-mxfx,wtm)=((1-menu(wrn))*((((pricegv((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricegv((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauk(wrn))*(eik(wrn,wtm))/(pi(wrn,wtm)))*mkti(wrn)))+((pricegv((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricegv((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauk(wrn))*(ejk(wrn,wtm))/(pj(wrn,wtm)))*mktj(wrn)))+((pricegv((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricegv((wrn*mxfx)+wmx-mxfx,wtm)/(pk(wrn,wtm)))*mktk(wrn))))))^risk(wrn); profitgd((wrn*mxfx)+wmx-mxfx,wtm)=((1-menu(wrn))*((((pricegd((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricegd((wrn*mxfx)+wmx-mxfx,wtm)/(pi(wrn,wtm)))*mkti(wrn)))+((pricegd((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricegd((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauij(wrn))*((1)/(eij(wrn,wtm)))/(pj(wrn,wtm)))*mktj(wrn)))+((pricegd((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricegd((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauk(wrn))*((1)/(eik(wrn,wtm)))/(pk(wrn,wtm)))*mktk(wrn))))))^risk(wrn); profitddv((wrn*mxfx)+wmx-mxfx,wtm)=(((1-menu(wrn))^2)*((((pricedindf((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricedindf((wrn*mxfx)+wmx-mxfx,wtm)/(pi(wrn,wtm)))*mkti(wrn)))+((pricedindf((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricedindf((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauij(wrn))*((1)/(eij(wrn,wtm)))/(pj(wrn,wtm)))*mktj(wrn)))+((pricevinv((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricevinv((wrn*mxfx)+wmx-mxfx,wtm)/(pk(wrn,wtm)))*mktk(wrn))))))^risk(wrn); profitdvv((wrn*mxfx)+wmx-mxfx,wtm)=(((1-menu(wrn))^2)*((((pricedind((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricedind((wrn*mxfx)+wmx-mxfx,wtm)/(pi(wrn,wtm)))*mkti(wrn)))+((pricevinfv((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricevinfv((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauk(wrn))*((1)*(ejk(wrn,wtm)))/(pj(wrn,wtm)))*mktj(wrn)))+((pricevinfv((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricevinfv((wrn*mxfx)+wmx-mxfx,wtm)/(pk(wrn,wtm)))*mktk(wrn))))))^risk(wrn); profitffv((wrn*mxfx)+wmx-mxfx,wtm)=(((1-menu(wrn))^2)*((((pricefindf((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauij(wrn))*(eij(wrn,wtm-1)))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricefindf((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauij(wrn))*((1)*(eij(wrn,wtm)))/(pi(wrn,wtm)))*mkti(wrn)))+((pricefindf((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauij(wrn))*(eij(wrn,wtm-1)))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricefindf((wrn*mxfx)+wmx-mxfx,wtm)/(pj(wrn,wtm)))*mktj(wrn)))+((pricevinv((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricevinv((wrn*mxfx)+wmx-mxfx,wtm)/(pk(wrn,wtm)))*mktk(wrn))))))^risk(wrn); profitdfd((wrn*mxfx)+wmx-mxfx,wtm)=(((1-menu(wrn))^2)*((((pricedindv((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricedindv((wrn*mxfx)+wmx-mxfx,wtm)/(pi(wrn,wtm)))*mkti(wrn)))+((pricefinf((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauij(wrn))*(eij(wrn,wtm-1)))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricefinf((wrn*mxfx)+wmx-mxfx,wtm)/(pj(wrn,wtm)))*mktj(wrn)))+((pricedindv((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricedindv((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauk(wrn))*((1)/(eik(wrn,wtm)))/(pk(wrn,wtm)))*mktk(wrn))))))^risk(wrn); profitvfv((wrn*mxfx)+wmx-mxfx,wtm)=(((1-menu(wrn))^2)*((((pricevindv((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricevindv((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauk(wrn))*((1)*(eik(wrn,wtm)))/(pi(wrn,wtm)))*mkti(wrn)))+((pricefinf((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauij(wrn))*(eij(wrn,wtm-1)))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricefinf((wrn*mxfx)+wmx-mxfx,wtm)/(pj(wrn,wtm)))*mktj(wrn)))+((pricevindv((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricevindv((wrn*mxfx)+wmx-mxfx,wtm)/(pk(wrn,wtm)))*mktk(wrn))))))^risk(wrn); profitdfv((wrn*mxfx)+wmx-mxfx,wtm)=(((1-menu(wrn))^3)*((((pricedind((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricedind((wrn*mxfx)+wmx-mxfx,wtm)/(pi(wrn,wtm)))*mkti(wrn)))+((pricefinf((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauij(wrn))*(eij(wrn,wtm-1)))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricefinf((wrn*mxfx)+wmx-mxfx,wtm)/(pj(wrn,wtm)))*mktj(wrn)))+((pricevinv((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricevinv((wrn*mxfx)+wmx-mxfx,wtm)/(pk(wrn,wtm)))*mktk(wrn))))))^risk(wrn); profitgv3p((wrn*mxfx)+wmx-mxfx,wtm)=(((1-menu(wrn))^3)*((((pricevind3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricevind3p((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauk(wrn))*(eik(wrn,wtm))/(pi(wrn,wtm)))*mkti(wrn)))+((pricevinf3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricevinf3p((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauk(wrn))*(ejk(wrn,wtm))/(pj(wrn,wtm)))*mktj(wrn)))+((pricevinv3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricevinv3p((wrn*mxfx)+wmx-mxfx,wtm)/(pk(wrn,wtm)))*mktk(wrn))))))^risk(wrn); profitgd3p((wrn*mxfx)+wmx-mxfx,wtm)=(((1-menu(wrn))^3)*((((pricedind3p((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricedind3p((wrn*mxfx)+wmx-mxfx,wtm)/(pi(wrn,wtm)))*mkti(wrn)))+((pricedinf3p((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricedinf3p((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauij(wrn))*((1)/(eij(wrn,wtm)))/(pj(wrn,wtm)))*mktj(wrn)))+((pricedinv3p((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricedinv3p((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauk(wrn))*((1)/(eik(wrn,wtm)))/(pk(wrn,wtm)))*mktk(wrn))))))^risk(wrn); profitddv3p((wrn*mxfx)+wmx-mxfx,wtm)=(((1-menu(wrn))^3)*((((pricedind3p((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricedind3p((wrn*mxfx)+wmx-mxfx,wtm)/(pi(wrn,wtm)))*mkti(wrn)))+((pricedinf3p((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricedinf3p((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauij(wrn))*((1)/(eij(wrn,wtm)))/(pj(wrn,wtm)))*mktj(wrn)))+((pricevinv3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricevinv3p((wrn*mxfx)+wmx-mxfx,wtm)/(pk(wrn,wtm)))*mktk(wrn))))))^risk(wrn); profitdvv3p((wrn*mxfx)+wmx-mxfx,wtm)=(((1-menu(wrn))^3)*((((pricedind3p((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricedind3p((wrn*mxfx)+wmx-mxfx,wtm)/(pi(wrn,wtm)))*mkti(wrn)))+((pricevinf3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricevinf3p((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauk(wrn))*((1)*(ejk(wrn,wtm)))/(pj(wrn,wtm)))*mktj(wrn)))+((pricevinv3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricevinv3p((wrn*mxfx)+wmx-mxfx,wtm)/(pk(wrn,wtm)))*mktk(wrn))))))^risk(wrn); profitffv3p((wrn*mxfx)+wmx-mxfx,wtm)=(((1-menu(wrn))^3)*((((pricefind3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauij(wrn))*(eij(wrn,wtm-1)))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricefind3p((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauij(wrn))*((1)*(eij(wrn,wtm)))/(pi(wrn,wtm)))*mkti(wrn)))+((pricefinf3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauij(wrn))*(eij(wrn,wtm-1)))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricefinf3p((wrn*mxfx)+wmx-mxfx,wtm)/(pj(wrn,wtm)))*mktj(wrn)))+((pricevinv3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricevinv3p((wrn*mxfx)+wmx-mxfx,wtm)/(pk(wrn,wtm)))*mktk(wrn))))))^risk(wrn); profitdfd3p((wrn*mxfx)+wmx-mxfx,wtm)=(((1-menu(wrn))^3)*((((pricedind3p((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricedind3p((wrn*mxfx)+wmx-mxfx,wtm)/(pi(wrn,wtm)))*mkti(wrn)))+((pricefinf3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauij(wrn))*(eij(wrn,wtm-1)))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricefinf3p((wrn*mxfx)+wmx-mxfx,wtm)/(pj(wrn,wtm)))*mktj(wrn)))+((pricedinv3p((wrn*mxfx)+wmx-mxfx,wtm))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricedinv3p((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauk(wrn))*((1)/(eik(wrn,wtm)))/(pk(wrn,wtm)))*mktk(wrn))))))^risk(wrn); profitvfv3p((wrn*mxfx)+wmx-mxfx,wtm)=(((1-menu(wrn))^3)*((((pricevind3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ai-(bi(wrn)*(pricevind3p((wrn*mxfx)+wmx-mxfx,wtm)*(1+tauk(wrn))*((1)*(eik(wrn,wtm)))/(pi(wrn,wtm)))*mkti(wrn)))+((pricefinf3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauij(wrn))*(eij(wrn,wtm-1)))-wage(wrn,wtm))*(Aj-(bj(wrn)*(pricefinf3p((wrn*mxfx)+wmx-mxfx,wtm)/(pj(wrn,wtm)))*mktj(wrn)))+((pricevinv3p((wrn*mxfx)+wmx-mxfx,wtm)*(1-tauk(wrn))*(eik(wrn,wtm-1)))-wage(wrn,wtm))*(Ak-(bk(wrn)*(pricevinv3p((wrn*mxfx)+wmx-mxfx,wtm)/(pk(wrn,wtm)))*mktk(wrn))))))^risk(wrn); end end end end %creating the total profit sumprofitgv=zeros(run,mxfx); sumprofitgd=zeros(run,mxfx); sumprofitddv=zeros(run,mxfx); sumprofitdvv=zeros(run,mxfx); sumprofitffv=zeros(run,mxfx); sumprofitdfd=zeros(run,mxfx); sumprofitvfv=zeros(run,mxfx); sumprofitdfv=zeros(run,mxfx); sumprofitgv3p=zeros(run,mxfx); sumprofitgd3p=zeros(run,mxfx); sumprofitddv3p=zeros(run,mxfx); sumprofitdvv3p=zeros(run,mxfx); sumprofitffv3p=zeros(run,mxfx); sumprofitdfd3p=zeros(run,mxfx); sumprofitvfv3p=zeros(run,mxfx); for vrn=1:run for vmx=1:mxfx for vtm=48:time sumprofitgv(vrn,vmx)=sumprofitgv(vrn,vmx)+profitgv((vrn*mxfx)+vmx-mxfx,vtm); sumprofitgd(vrn,vmx)=sumprofitgd(vrn,vmx)+profitgd((vrn*mxfx)+vmx-mxfx,vtm); sumprofitddv(vrn,vmx)=sumprofitddv(vrn,vmx)+profitddv((vrn*mxfx)+vmx-mxfx,vtm); sumprofitdvv(vrn,vmx)=sumprofitdvv(vrn,vmx)+profitdvv((vrn*mxfx)+vmx-mxfx,vtm); sumprofitffv(vrn,vmx)=sumprofitffv(vrn,vmx)+profitffv((vrn*mxfx)+vmx-mxfx,vtm); sumprofitdfd(vrn,vmx)=sumprofitdfd(vrn,vmx)+profitdfd((vrn*mxfx)+vmx-mxfx,vtm); sumprofitvfv(vrn,vmx)=sumprofitvfv(vrn,vmx)+profitvfv((vrn*mxfx)+vmx-mxfx,vtm); sumprofitdfv(vrn,vmx)=sumprofitdfv(vrn,vmx)+profitdfv((vrn*mxfx)+vmx-mxfx,vtm); sumprofitgv3p(vrn,vmx)=sumprofitgv3p(vrn,vmx)+profitgv3p((vrn*mxfx)+vmx-mxfx,vtm); sumprofitgd3p(vrn,vmx)=sumprofitgd3p(vrn,vmx)+profitgd3p((vrn*mxfx)+vmx-mxfx,vtm); sumprofitddv3p(vrn,vmx)=sumprofitddv3p(vrn,vmx)+profitddv3p((vrn*mxfx)+vmx-mxfx,vtm); sumprofitdvv3p(vrn,vmx)=sumprofitdvv3p(vrn,vmx)+profitdvv3p((vrn*mxfx)+vmx-mxfx,vtm); sumprofitffv3p(vrn,vmx)=sumprofitffv3p(vrn,vmx)+profitffv3p((vrn*mxfx)+vmx-mxfx,vtm); sumprofitdfd3p(vrn,vmx)=sumprofitdfd3p(vrn,vmx)+profitdfd3p((vrn*mxfx)+vmx-mxfx,vtm); sumprofitvfv3p(vrn,vmx)=sumprofitvfv3p(vrn,vmx)+profitvfv3p((vrn*mxfx)+vmx-mxfx,vtm); end end end %3 different results %1 best profit for all non-3ps %2 best profit for all 3ps + dfv %3 best profit for all prft functions bestpgv=[0,3,3,3,0,0]; bestpgd=[0,1,1,1,0,0]; bestpddv=[0,1,1,3,0,0]; bestpdvv=[0,1,3,3,0,0]; bestpffv=[0,2,2,3,0,0]; bestpdfd=[0,1,2,1,0,0]; bestpvfv=[0,3,2,3,0,0]; bestpdfv=[0,1,2,3,0,1]; bestpgv3p=[0,3,3,3,0,1]; bestpgd3p=[0,1,1,1,0,1]; bestpddv3p=[0,1,1,3,0,1]; bestpdvv3p=[0,1,3,3,0,1]; bestpffv3p=[0,2,2,3,0,1]; bestpdfd3p=[0,1,2,1,0,1]; bestpvfv3p=[0,3,2,3,0,1]; aaaY_1=zeros(run,6); aaaY_2=zeros(run,6); aaaY_3=zeros(run,6); for trn=1:run bestpgv=[0,3,3,3,0,0]; bestpgd=[0,1,1,1,0,0]; bestpddv=[0,1,1,3,0,0]; bestpdvv=[0,1,3,3,0,0]; bestpffv=[0,2,2,3,0,0]; bestpdfd=[0,1,2,1,0,0]; bestpvfv=[0,3,2,3,0,0]; bestpdfv=[0,1,2,3,0,1]; bestpgv3p=[0,3,3,3,0,1]; bestpgd3p=[0,1,1,1,0,1]; bestpddv3p=[0,1,1,3,0,1]; bestpdvv3p=[0,1,3,3,0,1]; bestpffv3p=[0,2,2,3,0,1]; bestpdfd3p=[0,1,2,1,0,1]; bestpvfv3p=[0,3,2,3,0,1]; for tmx=1:mxfx if sumprofitgv(trn,tmx)>bestpgv(1) bestpgv(1)=sumprofitgv(trn,tmx); bestpgv(5)=tmx; end if sumprofitgd(trn,tmx)>bestpgd(1) bestpgd(1)=sumprofitgd(trn,tmx); bestpgd(5)=tmx; end if sumprofitddv(trn,tmx)>bestpddv(1) bestpddv(1)=sumprofitddv(trn,tmx); bestpddv(5)=tmx; end if sumprofitdvv(trn,tmx)>bestpdvv(1) bestpdvv(1)=sumprofitdvv(trn,tmx); bestpdvv(5)=tmx; end if sumprofitffv(trn,tmx)>bestpffv(1) bestpffv(1)=sumprofitffv(trn,tmx); bestpffv(5)=tmx; end if sumprofitdfd(trn,tmx)>bestpdfd(1) bestpdfd(1)=sumprofitdfd(trn,tmx); bestpdfd(5)=tmx; end if sumprofitvfv(trn,tmx)>bestpvfv(1) bestpvfv(1)=sumprofitvfv(trn,tmx); bestpvfv(5)=tmx; end if sumprofitdfv(trn,tmx)>bestpdfv(1) bestpdfv(1)=sumprofitdfv(trn,tmx); bestpdfv(5)=tmx; end if sumprofitgv3p(trn,tmx)>bestpgv3p(1) bestpgv3p(1)=sumprofitgv3p(trn,tmx); bestpgv3p(5)=tmx; end if sumprofitgd3p(trn,tmx)>bestpgd3p(1) bestpgd3p(1)=sumprofitgd3p(trn,tmx); bestpgd3p(5)=tmx; end if sumprofitddv3p(trn,tmx)>bestpddv3p(1) bestpddv3p(1)=sumprofitddv3p(trn,tmx); bestpddv3p(5)=tmx; end if sumprofitdvv3p(trn,tmx)>bestpdvv3p(1) bestpdvv3p(1)=sumprofitdvv3p(trn,tmx); bestpdvv3p(5)=tmx; end if sumprofitffv3p(trn,tmx)>bestpffv3p(1) bestpffv3p(1)=sumprofitffv3p(trn,tmx); bestpffv3p(5)=tmx; end if sumprofitdfd3p(trn,tmx)>bestpdfd3p(1) bestpdfd3p(1)=sumprofitdfd3p(trn,tmx); bestpdfd3p(5)=tmx; end if sumprofitvfv3p(trn,tmx)>bestpvfv3p(1) bestpvfv3p(1)=sumprofitvfv3p(trn,tmx); bestpvfv3p(5)=tmx; end end if bestpgv(1)>bestpgd(1); aaaY_1(trn,1)=bestpgv(1); aaaY_1(trn,2)=bestpgv(2); aaaY_1(trn,3)=bestpgv(3); aaaY_1(trn,4)=bestpgv(4); aaaY_1(trn,5)=bestpgv(5); aaaY_1(trn,6)=bestpgv(6); else aaaY_1(trn,1)=bestpgd(1); aaaY_1(trn,2)=bestpgd(2); aaaY_1(trn,3)=bestpgd(3); aaaY_1(trn,4)=bestpgd(4); aaaY_1(trn,5)=bestpgd(5); aaaY_1(trn,6)=bestpgd(6); end if bestpddv(1)>aaaY_1(trn,1) aaaY_1(trn,1)=bestpddv(1); aaaY_1(trn,2)=bestpddv(2); aaaY_1(trn,3)=bestpddv(3); aaaY_1(trn,4)=bestpddv(4); aaaY_1(trn,5)=bestpddv(5); aaaY_1(trn,6)=bestpddv(6); end if bestpdvv(1)>aaaY_1(trn,1) aaaY_1(trn,1)=bestpdvv(1); aaaY_1(trn,2)=bestpdvv(2); aaaY_1(trn,3)=bestpdvv(3); aaaY_1(trn,4)=bestpdvv(4); aaaY_1(trn,5)=bestpdvv(5); aaaY_1(trn,6)=bestpdvv(6); end if bestpffv(1)>aaaY_1(trn,1) aaaY_1(trn,1)=bestpffv(1); aaaY_1(trn,2)=bestpffv(2); aaaY_1(trn,3)=bestpffv(3); aaaY_1(trn,4)=bestpffv(4); aaaY_1(trn,5)=bestpffv(5); aaaY_1(trn,6)=bestpffv(6); end if bestpdfd(1)>aaaY_1(trn,1) aaaY_1(trn,1)=bestpdfd(1); aaaY_1(trn,2)=bestpdfd(2); aaaY_1(trn,3)=bestpdfd(3); aaaY_1(trn,4)=bestpdfd(4); aaaY_1(trn,5)=bestpdfd(5); aaaY_1(trn,6)=bestpdfd(6); end if bestpvfv(1)>aaaY_1(trn,1) aaaY_1(trn,1)=bestpvfv(1); aaaY_1(trn,2)=bestpvfv(2); aaaY_1(trn,3)=bestpvfv(3); aaaY_1(trn,4)=bestpvfv(4); aaaY_1(trn,5)=bestpvfv(5); aaaY_1(trn,6)=bestpvfv(6); end if bestpdfv(1)>aaaY_1(trn,1) aaaY_1(trn,1)=bestpdfv(1); aaaY_1(trn,2)=bestpdfv(2); aaaY_1(trn,3)=bestpdfv(3); aaaY_1(trn,4)=bestpdfv(4); aaaY_1(trn,5)=bestpdfv(5); aaaY_1(trn,6)=bestpdfv(6); end if bestpgv3p(1)>bestpgd3p(1); aaaY_2(trn,1)=bestpgv3p(1); aaaY_2(trn,2)=bestpgv3p(2); aaaY_2(trn,3)=bestpgv3p(3); aaaY_2(trn,4)=bestpgv3p(4); aaaY_2(trn,5)=bestpgv3p(5); aaaY_2(trn,6)=bestpgv3p(6); else aaaY_2(trn,1)=bestpgd3p(1); aaaY_2(trn,2)=bestpgd3p(2); aaaY_2(trn,3)=bestpgd3p(3); aaaY_2(trn,4)=bestpgd3p(4); aaaY_2(trn,5)=bestpgd3p(5); aaaY_2(trn,6)=bestpgd3p(6); end if bestpddv3p(1)>aaaY_2(trn,1) aaaY_2(trn,1)=bestpddv3p(1); aaaY_2(trn,2)=bestpddv3p(2); aaaY_2(trn,3)=bestpddv3p(3); aaaY_2(trn,4)=bestpddv3p(4); aaaY_2(trn,5)=bestpddv3p(5); aaaY_2(trn,6)=bestpddv3p(6); end if bestpdvv3p(1)>aaaY_2(trn,1) aaaY_2(trn,1)=bestpdvv3p(1); aaaY_2(trn,2)=bestpdvv3p(2); aaaY_2(trn,3)=bestpdvv3p(3); aaaY_2(trn,4)=bestpdvv3p(4); aaaY_2(trn,5)=bestpdvv3p(5); aaaY_2(trn,6)=bestpdvv3p(6); end if bestpffv3p(1)>aaaY_2(trn,1) aaaY_2(trn,1)=bestpffv3p(1); aaaY_2(trn,2)=bestpffv3p(2); aaaY_2(trn,3)=bestpffv3p(3); aaaY_2(trn,4)=bestpffv3p(4); aaaY_2(trn,5)=bestpffv3p(5); aaaY_2(trn,6)=bestpffv3p(6); end if bestpdfd3p(1)>aaaY_2(trn,1) aaaY_2(trn,1)=bestpdfd3p(1); aaaY_2(trn,2)=bestpdfd3p(2); aaaY_2(trn,3)=bestpdfd3p(3); aaaY_2(trn,4)=bestpdfd3p(4); aaaY_2(trn,5)=bestpdfd3p(5); aaaY_2(trn,6)=bestpdfd3p(6); end if bestpvfv3p(1)>aaaY_2(trn,1) aaaY_2(trn,1)=bestpvfv3p(1); aaaY_2(trn,2)=bestpvfv3p(2); aaaY_2(trn,3)=bestpvfv3p(3); aaaY_2(trn,4)=bestpvfv3p(4); aaaY_2(trn,5)=bestpvfv3p(5); aaaY_2(trn,6)=bestpvfv3p(6); end if bestpdfv(1)>aaaY_2(trn,1) aaaY_2(trn,1)=bestpdfv(1); aaaY_2(trn,2)=bestpdfv(2); aaaY_2(trn,3)=bestpdfv(3); aaaY_2(trn,4)=bestpdfv(4); aaaY_2(trn,5)=bestpdfv(5); aaaY_2(trn,6)=bestpdfv(6); end if aaaY_2(trn,1)>=aaaY_1(trn,1) aaaY_3(trn,1)=aaaY_2(trn,1); aaaY_3(trn,2)=aaaY_2(trn,2); aaaY_3(trn,3)=aaaY_2(trn,3); aaaY_3(trn,4)=aaaY_2(trn,4); aaaY_3(trn,5)=aaaY_2(trn,5); aaaY_3(trn,6)=aaaY_2(trn,6); else aaaY_3(trn,1)=aaaY_1(trn,1); aaaY_3(trn,2)=aaaY_1(trn,2); aaaY_3(trn,3)=aaaY_1(trn,3); aaaY_3(trn,4)=aaaY_1(trn,4); aaaY_3(trn,5)=aaaY_1(trn,5); aaaY_3(trn,6)=aaaY_1(trn,6); end end %creating the vector for X in Y=XB+e aaaxreg=zeros(25,run); for urn=1:run aaaxreg(1,urn)=tauij(urn); aaaxreg(2,urn)=tauk(urn); aaaxreg(3,urn)=winfav(urn); aaaxreg(4,urn)=menu(urn); aaaxreg(5,urn)=mkti(urn); aaaxreg(6,urn)=mktj(urn); aaaxreg(7,urn)=mktk(urn); aaaxreg(8,urn)=bav(urn); aaaxreg(9,urn)=bvl(urn); aaaxreg(10,urn)=bi(urn); aaaxreg(11,urn)=bj(urn); aaaxreg(12,urn)=bk(urn); aaaxreg(13,urn)=pvlav(urn); aaaxreg(14,urn)=risk(urn); aaaxreg(15,urn)=varlpi(1,urn); aaaxreg(16,urn)=varlpj(1,urn); aaaxreg(17,urn)=varlpk(1,urn); aaaxreg(18,urn)=varleij(1,urn); aaaxreg(19,urn)=varleik(1,urn); aaaxreg(20,urn)=varlejk(1,urn); aaaxreg(21,urn)=wage(urn,1); aaaxreg(22,urn)=loop(1,urn); aaaxreg(23,urn)=corrpipj(urn,1); aaaxreg(24,urn)=corrpipk(urn,1); aaaxreg(25,urn)=corrpjpk(urn,1); end aaaxreg=aaaxreg'; toc beep