(*MATHEMATICA - grafy celkových nákladů*)

(*Předpokládaná doba provozu elektrárny [rok]*)
(Subscript[T, j])=60;
(Subscript[T, u])=40;
(Subscript[T, p])=30;

(*Produkce emisí CO2 [t/MWh]*)
(Subscript[Em, j])=0;
(Subscript[Em, uh])=1.2;
(Subscript[Em, uč])=0.9;
(Subscript[Em, p])=0.33;

(*Roční doba využití maxima [%]*)
(Subscript[k, j])=85;
(Subscript[k, u])=55;
(Subscript[k, p])=40;

(*Přepočet oční doby využití maxima na hodiny/rok*)
(Subscript[h, j])=(((Subscript[k, j])*8760)/100);
(Subscript[h, u])=(((Subscript[k, u])*8760)/100);
(Subscript[h, p])=(((Subscript[k, p])*8760)/100);

(*Výpočet poměrné anuity [-]*)
(Subscript[a, j])=(((r/100)*(((1+(r/100)))^(Subscript[T, j])))/((((1+(r/100)))^(Subscript[T, j]))-1));
(Subscript[a, u])=(((r/100)*(((1+(r/100)))^(Subscript[T, u])))/((((1+(r/100)))^(Subscript[T, u]))-1));
(Subscript[a, p])=(((r/100)*(((1+(r/100)))^(Subscript[T, p])))/((((1+(r/100)))^(Subscript[T, p]))-1));

(*Počáteční investice [EUR/MW]*)
(Subscript[IN, j])=2800000;
(Subscript[IN, uh])=1400000;
(Subscript[IN, uč])=1600000;
(Subscript[IN, p])=1000000;

(*Náklady na provoz a údržbu [EUR/MW]*)
(Subscript[Nprovoz, j])=0.025*(Subscript[IN, j]);
(Subscript[Nprovoz, uh])=0.022*(Subscript[IN, uh]);
(Subscript[Nprovoz, uč])=0.022*(Subscript[IN, uč]);
(Subscript[Nprovoz, p])=0.015*(Subscript[IN, p]);

(*Náklady na vyřazení z provozu [EUR/MW]*)
(Subscript[Nvyřazení, j])=0.15*(Subscript[IN, j]);
(Subscript[Nvyřazení, uh])=0.05*(Subscript[IN, uh]);
(Subscript[Nvyřazení, uč])=0.05*(Subscript[IN, uč]);
(Subscript[Nvyřazení, p])=0.05*(Subscript[IN, p]);

(*Náklady na palivo [EUR/MWh]*)
(Subscript[Npalivo, j])=((100+kp)/100)*6.5;
(Subscript[Npalivo, uh])=((100+kp)/100)*9;
(Subscript[Npalivo, uč])=((100+kp)/100)*15;
(Subscript[Npalivo, p])=((100+kp)/100)*35;

(*Náklady na emisní povolenky [EUR/MWh]*)
(Subscript[Nemise, j])=povolenka*(Subscript[Em, j]);
(Subscript[Nemise, uh])=povolenka*(Subscript[Em, uh]);
(Subscript[Nemise, uč])=povolenka*(Subscript[Em, uč]);
(Subscript[Nemise, p])=povolenka*(Subscript[Em, p]); 

(*Fixní náklady [EUR/MW]*)
(Subscript[FN, j])=((Subscript[IN, j])+(Subscript[Nvyřazení, j]))*(Subscript[a, j])+(Subscript[Nprovoz, j]);
(Subscript[FN, uh])=((Subscript[IN, uh])+(Subscript[Nvyřazení, uh]))*(Subscript[a, u])+(Subscript[Nprovoz, uh]);
(Subscript[FN, uč])=((Subscript[IN, uč])+(Subscript[Nvyřazení, uč]))*(Subscript[a, u])+(Subscript[Nprovoz, uč]);
(Subscript[FN, p])=((Subscript[IN, p])+(Subscript[Nvyřazení, p]))*(Subscript[a, p])+(Subscript[Nprovoz, p]);

(*Variabilní náklady [EUR/MWh]*)
(Subscript[VN, j])=(Subscript[Nemise, j])+(Subscript[Npalivo, j]);
(Subscript[VN, uh])=(Subscript[Nemise, uh])+(Subscript[Npalivo, uh]);
(Subscript[VN, uč])=(Subscript[Nemise, uč])+(Subscript[Npalivo, uč]);
(Subscript[VN, p])=(Subscript[Nemise, p])+(Subscript[Npalivo, p]);

(*Celkové náklady [EUR/MWh]*)
(Subscript[CN, j])=(Subscript[FN, j])/(Subscript[h, j])+(Subscript[VN, j]);
(Subscript[CN, uh])=(Subscript[FN, uh])/(Subscript[h, u])+(Subscript[VN, uh]);
(Subscript[CN, uč])=(Subscript[FN, uč])/(Subscript[h, u])+(Subscript[VN, uč]);
(Subscript[CN, p])=(Subscript[FN, p])/(Subscript[h, p])+(Subscript[VN, p]);

(*Graf - Závíslost celkových nákladů na zvoleném diskontu*)
(Subscript[CN, jn])=Evaluate[(Subscript[CN, j])/.{kp->1,povolenka->7}];
(Subscript[CN, uhn])=Evaluate[(Subscript[CN, uh])/.{kp->1,povolenka->7}];
(Subscript[CN, učn])=Evaluate[(Subscript[CN, uč])/.{kp->1,povolenka->7}];
(Subscript[CN, pn])=Evaluate[(Subscript[CN, p])/.{kp->1,povolenka->7}];
Needs["PlotLegends`"]
Plot[{(Subscript[CN, jn]),(Subscript[CN, uhn]),(Subscript[CN, učn]),(Subscript[CN, pn])},{r,2,10},
PlotStyle->{Green,Cyan,Blue,Red},
PlotRange->{{1,11},{25,75}},
GridLines->{None,Automatic},
AxesLabel->{"Diskont (%)","CN (EUR/MWh)"},PlotLegend->{"Jaderná el.","Uhelná el.- hnědé uhlí","Uhelná el.- černé uhlí","Paroplynová el."}]


(*Graf - Závíslost celkových nákladů na ceně emisních povolenek*)
(Subscript[CN, jn])=Evaluate[(Subscript[CN, j])/.{kp->1,r->6}];
(Subscript[CN, uhn])=Evaluate[(Subscript[CN, uh])/.{kp->1,r->6}];
(Subscript[CN, učn])=Evaluate[(Subscript[CN, uč])/.{kp->1,r->6}];
(Subscript[CN, pn])=Evaluate[(Subscript[CN, p])/.{kp->1,r->6}];
Needs["PlotLegends`"]
Plot[{(Subscript[CN, jn]),(Subscript[CN, uhn]),(Subscript[CN, učn]),(Subscript[CN, pn])},{povolenka,0,20},
PlotStyle->{Green,Cyan,Blue,Red},
PlotRange->{{0,22},{25,75}},
GridLines->{None,Automatic},
AxesLabel->{"Emisní povolenka (EUR/t CO2)","CN (EUR/MWh)"},PlotLegend->{"Jaderná el.","Uhelná el.- hnědé uhlí","Uhelná el.- černé uhlí","Paroplynová el."}]

(*Graf - Závíslost celkových nákladů na relativní změně ceny paliva*)
(Subscript[CN, jn])=Evaluate[(Subscript[CN, j])/.{r->6,povolenka->7}];
(Subscript[CN, uhn])=Evaluate[(Subscript[CN, uh])/.{r->6,povolenka->7}];
(Subscript[CN, učn])=Evaluate[(Subscript[CN, uč])/.{r->6,povolenka->7}];
(Subscript[CN, pn])=Evaluate[(Subscript[CN, p])/.{r->6,povolenka->7}];
Needs["PlotLegends`"]
Plot[{(Subscript[CN, jn]),(Subscript[CN, uhn]),(Subscript[CN, učn]),(Subscript[CN, pn])},{kp,-50,50},
PlotStyle->{Green,Cyan,Blue,Red},
PlotRange->{{-60,60},{25,90}},
GridLines->{None,Automatic},
AxesLabel->{"Relativní změna ceny paliva (%)","CN (EUR/MWh)"},PlotLegend->{"Jaderná el.","Uhelná el.- hnědé uhlí","Uhelná el.- černé uhlí","Paroplynová el."}]

(*Fixní náklady číselně [EUR/MW]*)
(Subscript[FN, j])=269240;
(Subscript[FN, uh])=128498;
(Subscript[FN, uč])=146855;
(Subscript[FN, p])=91281;

(*Variabilní náklady číselně [EUR/MWh]*)
(Subscript[VN, j])=6.5;
(Subscript[VN, uh])=17.4;
(Subscript[VN, uč])=21.3;
(Subscript[VN, p])=37.31;

(*Graf - Závislost celkových nákladů na roční době využití maxima elektrárny*)
Needs["PlotLegends`"]
Plot[{(Subscript[FN, j])/h+(Subscript[VN, j]),(Subscript[FN, uh])/h+(Subscript[VN, uh]),(Subscript[FN, uč])/h+(Subscript[VN, uč]),(Subscript[FN, p])/h+(Subscript[VN, p])},{h,500,8760},
PlotRange->{0,200},
PlotStyle->{Green,Cyan,Blue,Red},
AxesLabel->{"Roční doba využití maxima (h/rok)","CN (EUR/MWh)"},
GridLines->{None,Automatic},
PlotLegend->{"Jaderná el.","Uhelná el.- hnědé uhlí","Uhelná el.- černé uhlí","Paroplynová el."}]

(*Graf - Závislost celkových nákladů na provoz 1 MW instalovaného výkonu na roční době využití maxima elektrárny*)
Plot[{(Subscript[FN, j])+(Subscript[VN, j])*h,(Subscript[FN, uh])+(Subscript[VN, uh])*h,(Subscript[FN, uč])+(Subscript[VN, uč])*h,(Subscript[FN, p])+(Subscript[VN, p])*h},{h,500,8760},
PlotStyle->{Green,Cyan,Blue,Red},
GridLines->{None,Automatic},
AxesLabel->{"Roční doba využití maxima (h/rok)","CN na provoz instalovaného výkonu 1 MW (EUR/MW)"},PlotLegend->{"Jaderná el.","Uhelná el.- hnědé uhlí","Uhelná el.- černé uhlí","Paroplynová el."}]