% symm_ss_set.m - Setzt die Werte des symmetrischen SS nachdem man mit
% symm_ss den SS-Zins Rss gefunden hat.


function ys=symm_ss_set(Rss)
global M_                  % D. h., dass M_ im folgenden als globale Variable aufzufassen ist.
                           % M_ wird schon in weg_steadystate als global deklariert. 
                           
for i=1:size(M_.params,1);
    setval(M_.param_names(i,:),M_.params(i,1));
end

lneab=0;
lnebc=0;
lnRa=ln(Rss);
lnRb=ln(Rss);
lnRc=ln(Rss);
rss=1/ba1-(1-daa);
wpss=(1-aaa)*(exp(lnAa)/(1+lapa)*(aaa/rss)^aaa)^(1/(1-aaa));
klss=wpss*aaa/(1-aaa)/rss;
ylss=exp(lnAa)*klss^aaa;
css=(ba1*nuc1a*(1/Rss-ba1)/tab1/eta1a)^(1/sa1);
pss=mass/3/css;
lass=tab1*eta1a/pss/(1/Rss-ba1);
lss=(wpss*pss*lass/(1+lala)/tan1)^(1/xi1);
yss=ylss*lss;
lnwa=ln(wpss*pss);           % Nominallöhne
lnwb=ln(wpss*pss); 
lnwc=ln(wpss*pss);
lnla=ln(lss);           % Arbeitsnachfrage von Firmen in a
lnlb=ln(lss);
lnlc=ln(lss);
lnn1=ln(lss);           % Arbeitsangebot Haushalt 1
lnn2=ln(lss);
lnn3=ln(lss);
ln_ra=ln(rss);          % Bruttorealzins
ln_rb=ln(rss);
ln_rc=ln(rss);
lnka=ln(klss*lss);           % Kapital
lnkb=ln(klss*lss); 
lnkc=ln(klss*lss);
lnmca=ln(1/(1+lapa));          % Marginal cost
lnmcb=ln(1/(1+lapb));
lnmcc=ln(1/(1+lapc));
lndpa=0;          % price dispersion
lndpb=0;
lndpc=0;
lninfa=0;         % Inflation
lninfb=0;
lninfc=0;
lnpa=ln(pss);                                                                                    
lnpb=ln(pss); 
lnpc=ln(pss);
lnpsa=ln(pss);          % Optimaler Preis in A geg. die Preisrigidität                                                                          
lnpsb=ln(pss); 
lnpsc=ln(pss);
lnla1=ln(lass);                         % Lagrange-Multiplikatoren Budgetrestriktion
lnla2=ln(lass); 
lnla3=ln(lass); 
lnmu1=ln((1-ba1)/ba1)+lnla1;                       % Lagrange-Multiplikatoren Cash-in-advance constraint
lnmu2=ln((1-ba1)/ba1)+lnla1;  
lnmu3=ln((1-ba1)/ba1)+lnla1;  
lnq1=lnpa+lnla1;                        % Lagrange-Multiplikatoren Kapitalakkumulationsgleichung
lnq2=lnpa+lnla1;   
lnq3=lnpa+lnla1;  
lnifs1a=((1+lapa)/lapa*lnpa+ln(yss))-ln(1-ba1*xipa);
lnifs1b=lnifs1a;
lnifs1c=lnifs1a;
lnifs2a=(ln(pss)/lapa+ln(yss))-ln(1-ba1*xipa); 
lnifs2b=lnifs2a; 
lnifs2c=lnifs2a;
lnifs3a=ln(tan1*(1+lala)/lala)+(1+xi1)*((1+lala)/lala*lnwa+ln(lss))-ln(1-ba1*xiwa);
lnifs3b=lnifs3a;
lnifs3c=lnifs3a;
lnifs4a=ln(lass/lala*exp(lnwa)^((1+lala)/lala)*lss)-ln(1-ba1*xiwa);
lnifs4b=lnifs4a;
lnifs4c=lnifs4a;
lnwsa=lnwa;                       % Optimaler Lohn in A geg. die Lohnrigidität
lnwsb=lnwa; 
lnwsc=lnwa;
lnc1=ln(css);                        % Aggregierter Konsum
lnc2=ln(css);  
lnc3=ln(css);                         % Die länderspezifischen Konsumkomponenten stehen als Variablen von fsolve ganz oben 
lnc1a=lnc1;
lnc1b=lnc1a;
lnc1c=lnc1a;
lnc2a=lnc1a;
lnc2b=lnc1a;
lnc2c=lnc1a;
lnc3a=lnc1a;
lnc3b=lnc1a;
lnc3c=lnc1a;
lnya=ln(yss);                        % Produktion in Land A
lnyb=ln(yss); 
lnyc=ln(yss); 
lnia=ln(daa)+lnka;                        % Investitionsnachfrage in Land A 
lnib=ln(daa)+lnka;
lnic=ln(daa)+lnka;
lnm1=ln(mass);                        % Nominale Geldnachfragen der Haushalte
lnm2=ln(mass);  
lnm3=ln(mass);  
lnpia=ln(exp(lnpa+lnya)-exp(ln_ra+lnpa+lnka)/aaa);                       % Nominale Gewinne
lnpib=ln(exp(lnpa+lnya)-exp(ln_ra+lnpa+lnka)/aaa);
lnpic=ln(exp(lnpa+lnya)-exp(ln_ra+lnpa+lnka)/aaa);
lnta=ln(tass);                        % Reale Steuern
lntb=ln(tass);  
lntc=ln(tass);  
lnga=ln(gass);                          % reale Staatsausgaben
lngb=ln(gass);   
lngc=ln(gass);  
lneac=0;
eab=1;          % Ab hier die Niveaus
eac=1;
ebc=1;
pa=exp(lnpa);
pb=exp(lnpa);
pc=exp(lnpa);
dpa=1;                         % Price dispersion, nur in Approximation zweiter Ordnung wesentlich
dpb=1;
dpc=1;
ra=exp(ln_ra);                          % Bruttorealzins
rb=exp(ln_ra); 
rc=exp(ln_ra); 
infa=1;                        % Inflationsraten
infb=1;
infc=1;
ifs1a=exp(lnifs1a);                       % infinite sums for Calvo pricing 
ifs1b=exp(lnifs1a);
ifs1c=exp(lnifs1a);
ifs2a=exp(lnifs2a);
ifs2b=exp(lnifs2a);
ifs2c=exp(lnifs2a);
ifs3a=exp(lnifs3a);                       % infinite sums for Calvo wage setting
ifs3b=exp(lnifs3a);                      
ifs3c=exp(lnifs3a);
ifs4a=exp(lnifs4a);
ifs4b=exp(lnifs4a);
ifs4c=exp(lnifs4a);
m1=mass;                          % Geldnachfragen
m2=mass;
m3=mass;
la1=exp(lnla1);
la2=exp(lnla1);
la3=exp(lnla1);
mu1=exp(lnmu1);
mu2=exp(lnmu1);
mu3=exp(lnmu1);
q1=exp(lnq1);
q2=exp(lnq1);
q3=exp(lnq1);
pia=exp(lnpia);
pib=exp(lnpia);
pic=exp(lnpia);
% bca=pss*(tass-gass)/3/(1-1/Rss);                         % Bondnachfrage der Zentralbank von A
% bcb=pss*(tass-gass)/3/(1-1/Rss);
% bcc=pss*(tass-gass)/3/(1-1/Rss);
bca=0;
bcb=0;
bcc=0;
ta=tass;
tb=tass;
tc=tass;
ga=gass;
gb=gass;
gc=gass;
ma=mass;                          % Geldangebot in A
mb=mass;
mc=mass;
s1a=1/3;                         % Profit shares
s1b=1/3;  
s1c=1/3;  
s2a=1/3;  
s2b=1/3;  
s2c=1/3;    
s3a=1/3;  
s3b=1/3;  
s3c=1/3;  
x1b=(tax1b/lass)^(1/sax1b);                         % Nominal Remittances
x1c=(tax1b/lass)^(1/sax1b);
x2a=(tax1b/lass)^(1/sax1b);
x2c=(tax1b/lass)^(1/sax1b);
x3a=(tax1b/lass)^(1/sax1b);
x3b=(tax1b/lass)^(1/sax1b);
spa=pia/(1-ba1)*(ba1+ka1a/la1*tas1/pa);                         % share price for profits of firms in country a
spb=pia/(1-ba1)*(ba1+ka1a/la1*tas1/pa); 
spc=pia/(1-ba1)*(ba1+ka1a/la1*tas1/pa); 
b1a=pa*(ta-ga)/3/(1-1/exp(lnRa));                         % Bondnachfragen nach Bonds des eigenen Landes
b2b=b1a;
b3c=b1a;
b2a=b1a;                         % Bondnachfragen nach Bonds anderer Länder
b2c=b1a;
b1b=b1a;
b1c=b1a;
b3a=b1a;
b3b=b1a;
lnlncut1sa=ln(ln(3/2));    % cut1 und cut2 sind Zahlen zwischen 0 und 1 in der Form cut=exp(-exp(lnlncut)). 
lnlncut2sa=ln(ln(2));    % Daraus definieren wir uns drei zu Eins addierbare shares als s1=cut2*cut1, s1+s2=cut1, s3=1-s1-s2
lnlncut1sb=ln(ln(3/2));    % Damit sind die shares immer positiv und immer Eins in der Summe.
lnlncut2sb=ln(ln(2));    % Die explizite Formel der shares lautet: s1=cut1*cut2; s2=cut1*(1-cut2); s3=1-cut1; 
lnlncut1sc=ln(ln(3/2));    % sa bezeichnet die shares für die Gewinne von a, etc.  
lnlncut2sc=ln(ln(2)); 
lnx1b=ln(x1b);          % Transfers von Haushalt 1 in Land B
lnx1c=ln(x1b);
lnx2a=ln(x1b);
lnx2c=ln(x1b);
lnx3a=ln(x1b);
lnx3b=ln(x1b);
lnspa=ln(spa);         % share price for profits of firms in country a
lnspb=ln(spb);
lnspc=ln(spc);     

save lnya;
save lnc1a;
save lnga;
save lnia;




ys=[
lneab         % Wechselkurs Anzahl Einheiten A pro Einheit B
lnebc
lnRa          % Nominalzinsen
lnRb           
lnRc
lnpa 
lnpb
lnpc
lnc1a         % Konsum von Haushalt 1 des in A produzierten Gutes
lnc2b
lnc3c
b1a                         % Bondnachfragen nach Bonds des eigenen Landes
b2b
b3c
b2a                         % Bondnachfragen nach Bonds anderer Länder
b2c
b1b
b1c
b3a
b3b
lnlncut1sa    % cut1 und cut2 sind Zahlen zwischen 0 und 1 in der Form cut=exp(-exp(lnlncut)). 
lnlncut2sa    % Daraus definieren wir uns drei zu Eins addierbare shares als s1=cut2*cut1, s1+s2=cut1, s3=1-s1-s2
lnlncut1sb    % Damit sind die shares immer positiv und immer Eins in der Summe.
lnlncut2sb    % Die explizite Formel der shares lautet: s1=cut1*cut2; s2=cut1*(1-cut2); s3=1-cut1; 
lnlncut1sc    % sa bezeichnet die shares für die Gewinne von a, etc.  
lnlncut2sc    % Die Zahlen 1 und 2 haben nichts mit den Haushalten zu tun sondern erlauben die Berechnung von zwei Schnittpunkten im Intervall (0,1) 
lnx1b          % Transfers von Haushalt 1 in Land B
lnx1c
lnx2a
lnx2c
lnx3a
lnx3b
lnspa          % share price for profits of firms in country a
lnspb
lnspc         % Diese ersten Variablen (bis hier) werden in fsolve gefunden.
lnc1b         % Konsum der Güter anderer Länder
lnc1c
lnc2a
lnc2c
lnc3a
lnc3b
lnwa           % Nominallöhne
lnwb 
lnwc
lnla           % Arbeitsnachfrage von Firmen in a
lnlb
lnlc
lnn1           % Arbeitsangebot Haushalt 1
lnn2
lnn3
ln_ra          % Bruttorealzins
ln_rb
ln_rc
lnka           % Kapital
lnkb 
lnkc
lnmca          % Marginal cost
lnmcb
lnmcc
lndpa          % price dispersion
lndpb
lndpc
lninfa         % Inflation
lninfb
lninfc
lnpsa          % Optimaler Preis in A geg. die Preisrigidität                                                                          
lnpsb 
lnpsc
lnifs1a 
lnifs2a 
lnifs3a 
lnifs4a 
lnifs1b 
lnifs2b 
lnifs3b 
lnifs4b 
lnifs1c
lnifs2c
lnifs3c
lnifs4c
lnwsa                       % Optimaler Lohn in A geg. die Lohnrigidität
lnwsb 
lnwsc
lnc1                        % Aggregierter Konsum
lnc2 
lnc3                        % Die länderspezifischen Konsumkomponenten stehen als Variablen von fsolve ganz oben 
lnya                        % Produktion in Land A
lnyb
lnyc
lnia                        % Investitionsnachfrage in Land A 
lnib 
lnic
lnm1                        % Nominale Geldnachfragen der Haushalte
lnm2
lnm3
lnla1                       % Lagrange-Multiplikatoren Budgetrestriktion
lnla2
lnla3
lnmu1                       % Lagrange-Multiplikatoren Cash-in-advance constraint
lnmu2
lnmu3
lnq1                        % Lagrange-Multiplikatoren Kapitalakkumulationsgleichung
lnq2 
lnq3
lnpia                       % Nominale Gewinne
lnpib
lnpic
lnta                        % Reale Steuern
lntb
lntc
lnga                        % reale Staatsausgaben
lngb 
lngc
lneac
eab          % Ab hier die Niveaus
eac
ebc
pa
pb
pc
dpa                         % Price dispersion, nur in Approximation zweiter Ordnung wesentlich
dpb
dpc
ra                          % Bruttorealzins
rb
rc
infa                        % Inflationsraten
infb
infc
ifs1a                       % infinite sums for Calvo pricing 
ifs1b
ifs1c
ifs2a
ifs2b
ifs2c
ifs3a                       % infinite sums for Calvo wage setting
ifs3b                      
ifs3c
ifs4a
ifs4b
ifs4c
m1                          % Geldnachfragen
m2
m3
la1
la2
la3
mu1
mu2
mu3
q1
q2
q3
pia
pib
pic
bca                         % Bondnachfrage der Zentralbank von A
bcb
bcc
ta
tb
tc
ga
gb
gc
ma                          % Geldangebot in A
mb
mc
s1a                         % Profit shares
s1b
s1c
s2a
s2b
s2c  
s3a
s3b
s3c
x1b                         % Nominal Remittances
x1c
x2a
x2c
x3a
x3b
spa                         % share price for profits of firms in country a
spb
spc
];

save ys;
end


