Download An Integrated System Dynamics Model for

z(10) = ((FC1+tau1)*exp(short_El)./exp(a1p) - exp(TCe1)) - ((FC2+tau2)*exp(short_El)./exp(a2p) - exp(TCe2));
z(11) = ((FC2+tau2)*exp(short_El)./exp(a2p) - exp(TCe2)) - ((FC3+tau3)*exp(short_El)./exp(a3p) - exp(TCe3));
% Heat energy production
z(12) = exp(Eh) - Ah*(b1*exp(Fh1).^eta_h + b2*exp(Fh2).^eta_h + b3*exp(Fh3).^eta_h + b4*exp(Fh4).^eta_h).^(1./eta_h);
z(13) = exp(Ph) - (exp(TCh1) + exp(TCh2) + exp(TCh3) + exp(TCh4))./exp(Eh);
z(14) = ((FC1+tau1)*exp(short_H)./(b1*exp(Fh1).^(eta_h-1)) - exp(TCh1)) - ((FC2+tau2)*exp(short_H)./(b2*exp(Fh2).^(eta_h-1))
- exp(TCh2));
z(15) = ((FC2+tau2)*exp(short_H)./(b2*exp(Fh2).^(eta_h-1)) - exp(TCh2)) - ((FC3+tau3)*exp(short_H)./(b3*exp(Fh3).^(eta_h-1))
- exp(TCh3));
z(16) = ((FC3+tau3)*exp(short_H)./(b3*exp(Fh3).^(eta_h-1)) - exp(TCh3)) ((1+mu2)*exp(FC4)*exp(short_H)./(b4*exp(Fh4).^(eta_h-1)) - exp(TCh4));
% Nuclear and hydro electricity production
z(17) = exp(e4) - a4;
z(18) = exp(e5) - a5;
%Fuel prices
z(19) = fstar1 - (exp(Fe1) + exp(Fh1) + Exp1);
z(20) = fstar2 - (exp(Fe2) + exp(Fh2) + Exp2);
z(21) = fstar3 - (exp(Fe3) + exp(Fh3) + Exp3);
% Total cost functions for electricity
z(22) = exp(TCe1) - exp(r)*k1 - exp(Fe1)*(FC1+tau1);
z(23) = exp(TCe2) - exp(r)*k2 - exp(Fe2)*(FC2+tau2);
z(24) = exp(TCe3) - exp(r)*k3 - exp(Fe3)*(FC3+tau3);
% Total cost functions for heat energy
z(25) = exp(TCh1) - exp(Fh1)*(FC1+tau1);
108