(*MATHEMATICA - grafy kapacitních plateb*)

(*Fixní náklady číselně [EUR/MWh]*)
(Subscript[FN, z])=269240;
(Subscript[FN, sh])=128498;
(Subscript[FN, sč])=146855;
(Subscript[FN, š])=91281;

(*Variabilní náklady číselně [EUR/MWh]*)
(Subscript[VN, z])=6.5;
(Subscript[VN, sh])=17.4;
(Subscript[VN, sč])=21.3;
(Subscript[VN, š])=37.31;

(*Roční doba využití maxima [h/rok]*)
(Subscript[h, z])=7446;
(Subscript[h, s])=4818;
(Subscript[h, š])=3504;

(*Celkové náklady [EUR/MW]*)
(Subscript[CN, z])=(Subscript[FN, z])+(Subscript[VN, z])*(Subscript[h, z]);
(Subscript[CN, sh])=(Subscript[FN, sh])+(Subscript[VN, sh])*(Subscript[h, s]);
(Subscript[CN, sč])=(Subscript[FN, sč])+(Subscript[VN, sč])*(Subscript[h, s]);
(Subscript[CN, š])=(Subscript[FN, š])+(Subscript[VN, š])*(Subscript[h, š]);

(*Příjmy [EUR/MW]*)
(Subscript[P, z])=(Subscript[TC, z])*(Subscript[h, z]);
(Subscript[P, s])=(Subscript[TC, s])*(Subscript[h, s]);
(Subscript[P, š])=(Subscript[TC, š])*(Subscript[h, š]);

(*Zisk/Ztráta [EUR/MW]*)
(Subscript[Z, z])=(Subscript[P, z])-(Subscript[CN, z]);
(Subscript[Z, sh])=(Subscript[P, s])-(Subscript[CN, sh]);
(Subscript[Z, sč])=(Subscript[P, s])-(Subscript[CN, sč]);
(Subscript[Z, š])=(Subscript[P, š])-(Subscript[CN, š]);

(*Fixní kapacitní platby [EUR/MW]*)
FKP =If[(Subscript[FN, š]) -((Subscript[TC, š])-(Subscript[VN, š]))*(Subscript[h, š])<0,0,(Subscript[FN, š]) -((Subscript[TC, š])-(Subscript[VN, š]))*(Subscript[h, š])];
(*Fixní kapacitní platby se stropem 50 000 EUR/MW [EUR/MW]*)
OFKP=If[(Subscript[FN, š]) -((Subscript[TC, š])-(Subscript[VN, š]))*(Subscript[h, š])<0,0,If[(Subscript[FN, š]) -((Subscript[TC, š])-(Subscript[VN, š]))*(Subscript[h, š])>50000,50000,(Subscript[FN, š]) -((Subscript[TC, š])-(Subscript[VN, š]))*(Subscript[h, š])]];

(*Dynamické kapacitní platby [EUR/MW]*)
(Subscript[DKP, z])=If[(Subscript[FN, z]) -((Subscript[TC, z])-(Subscript[VN, z]))*(Subscript[h, z])<0,0,(Subscript[FN, z]) -((Subscript[TC, z])-(Subscript[VN, z]))*(Subscript[h, z])];
(Subscript[DKP, sh])=If[(Subscript[FN, sh]) -((Subscript[TC, s])-(Subscript[VN, sh]))*(Subscript[h, s])<0,0,(Subscript[FN, sh]) -((Subscript[TC, s])-(Subscript[VN, sh]))*(Subscript[h, s])];
(Subscript[DKP, sč])=If[(Subscript[FN, sč]) -((Subscript[TC, s])-(Subscript[VN, sč]))*(Subscript[h, s])<0,0,(Subscript[FN, sč]) -((Subscript[TC, s])-(Subscript[VN, sč]))*(Subscript[h, s])];
(Subscript[DKP, š])=If[(Subscript[FN, š]) -((Subscript[TC, š])-(Subscript[VN, š]))*(Subscript[h, š])<0,0,(Subscript[FN, š]) -((Subscript[TC, š])-(Subscript[VN, š]))*(Subscript[h, š])];
(*Dynamické kapacitní platby se stropem 50 000 EUR/MW [EUR/MW]*)
(Subscript[ODKP, z])=If[(Subscript[FN, z]) -((Subscript[TC, z])-(Subscript[VN, z]))*(Subscript[h, z])<0,0,If[(Subscript[FN, z]) -((Subscript[TC, z])-(Subscript[VN, z]))*(Subscript[h, z])>50000,50000,(Subscript[FN, z]) -((Subscript[TC, z])-(Subscript[VN, z]))*(Subscript[h, z])]];
(Subscript[ODKP, sh])=If[(Subscript[FN, sh]) -((Subscript[TC, s])-(Subscript[VN, sh]))*(Subscript[h, s])<0,0,If[(Subscript[FN, sh]) -((Subscript[TC, s])-(Subscript[VN, sh]))*(Subscript[h, s])>50000,50000,(Subscript[FN, sh]) -((Subscript[TC, s])-(Subscript[VN, sh]))*(Subscript[h, s])]];
(Subscript[ODKP, sč])=If[(Subscript[FN, sč]) -((Subscript[TC, s])-(Subscript[VN, sč]))*(Subscript[h, s])<0,0,If[(Subscript[FN, sč]) -((Subscript[TC, s])-(Subscript[VN, sč]))*(Subscript[h, s])>50000,50000,(Subscript[FN, sč]) -((Subscript[TC, s])-(Subscript[VN, sč]))*(Subscript[h, s])]];
(Subscript[ODKP, š])=If[(Subscript[FN, š]) -((Subscript[TC, š])-(Subscript[VN, š]))*(Subscript[h, š])<0,0,If[(Subscript[FN, š]) -((Subscript[TC, š])-(Subscript[VN, š]))*(Subscript[h, š])>50000,50000,(Subscript[FN, š]) -((Subscript[TC, š])-(Subscript[VN, š]))*(Subscript[h, š])]];

(*Průměrná tržní cena*)
(Subscript[TC, z]) = 1*(Subscript[TC, z]);
(Subscript[TC, s]) = 1.12*(Subscript[TC, z]);
(Subscript[TC, š]) = 1.2*(Subscript[TC, z]);

(*Graf - Zisk/Ztráta na 1 MWh vyrobené energie při změně ceny energie na burze*)
Needs["PlotLegends`"]
Plot[{(Subscript[Z, z])/(Subscript[h, z]),(Subscript[Z, sh])/(Subscript[h, s]),(Subscript[Z, sč])/(Subscript[h, s]),(Subscript[Z, š])/(Subscript[h, š])},{(Subscript[TC, z]),25,65},PlotStyle->{Green,Cyan,Blue,Red},AxesLabel->{"Baseload cena (EUR/MWh)","Zisk/Ztráta (EUR/MWh)"},PlotLegend->{"Jaderná el.","Uhelná el.- hnědé uhlí","Uhelná el.- černé uhlí","Paroplynová el."}]

(*Graf - Zisk/Ztráta na 1 MW instalovaného výkonu při změně ceny energie na burze*)
Plot[{(Subscript[Z, z]),(Subscript[Z, sh]),(Subscript[Z, sč]),(Subscript[Z, š])},{(Subscript[TC, z]),25,65},PlotStyle->{Green,Cyan,Blue,Red},AxesLabel->{"Baseload cena (EUR/MWh)","Zisk/Ztráta (EUR/MW)"},PlotLegend->{"Jaderná el.","Uhelná el.- hnědé uhlí","Uhelná el.- černé uhlí","Paroplynová el."}]

(*Graf - Vliv fixních kapacitních plateb na zisk/ztrátu elektráren*)
Plot[{(Subscript[Z, z])+FKP,(Subscript[Z, z]),(Subscript[Z, sh])+FKP,(Subscript[Z, sh]),(Subscript[Z, sč])+FKP,(Subscript[Z, sč]),(Subscript[Z, š])+FKP,(Subscript[Z, š]),FKP},{(Subscript[TC, z]),25,65},PlotStyle->{Green,{Green,Dashed},Cyan,{Cyan,Dashed},Blue,{Blue,Dashed},Red,{Red,Dashed},Black},AxesLabel->{"Baseload cena (EUR/MWh)","Zisk/Ztráta (EUR/MW)"},PlotLegend->{"Jaderná el. s FKP","Jaderná el. bez FKP","Uhelná el.- hnědé uhlí s FKP","Uhelná el.- hnědé uhlí bez FKP","Uhelná el.- černé uhlí s FKP","Uhelná el.- černé uhlí bez FKP","Paroplynová el. s FKP","Paroplynová el. bez FKP","Fixní kapacitní platby"}]

(*Graf - Vliv fixních kapacitních plateb se stropem 50 000 EUR/MW na zisk/ztrátu elektráren*)
Plot[{(Subscript[Z, z])+OFKP,(Subscript[Z, sh])+OFKP,(Subscript[Z, sč])+OFKP,(Subscript[Z, š])+OFKP,OFKP},{(Subscript[TC, z]),25,65},PlotStyle->{Green,Cyan,Blue,Red,Black},AxesLabel->{"Baseload cena (EUR/MWh)","Zisk/Ztráta (EUR/MW)"},PlotLegend->{"Jaderná el. s FKP","Uhelná el.- hnědé uhlí s FKP","Uhelná el.- černé uhlí s FKP","Paroplynová el. s FKP","Fixní kapacitní platby"}]

(*Graf - Vliv dynamických kapacitních plateb na zisk/ztrátu elektráren*)
Plot[{(Subscript[Z, z])+(Subscript[DKP, z]),(Subscript[Z, z]),(Subscript[Z, sh])+(Subscript[DKP, sh]),(Subscript[Z, sh]),(Subscript[Z, sč])+(Subscript[DKP, sč]),(Subscript[Z, sč]),(Subscript[Z, š])+(Subscript[DKP, š]),(Subscript[Z, š]),(Subscript[DKP, z]),(Subscript[DKP, sh]),(Subscript[DKP, sč]),(Subscript[DKP, š])},{(Subscript[TC, z]),30,65},PlotStyle->{Green,{Green,Dashed},Cyan,{Cyan,Dashed},Blue,{Blue,Dashed},Red,{Red,Dashed},{Green,DotDashed},{Cyan,DotDashed},{Blue,DotDashed},{Red,DotDashed}},
AxesLabel->{"Baseload cena (EUR/MWh)","Zisk/Ztráta (EUR/MW)"},PlotLegend->{"Jaderná el. s FKP","Jaderná el. bez FKP","Uhelná el.- hnědé uhlí s FKP","Uhelná el.- hnědé uhlí bez FKP","Uhelná el.- černé uhlí s FKP","Uhelná el.- černé uhlí bez FKP","Paroplynová el. s FKP","Paroplynová el. bez FKP","Dynamické kapacitní platby - jaderná el.","Dynamické kapacitní platby - hnědouhelná el.","Dynamické kapacitní platby - černouhelná el.","Dynamické kapacitní platby - paroplynová el."}]

(*Graf - Vliv dynamických kapacitních plateb na zisk/ztrátu elektráren na 1 MWh*)
Plot[{((Subscript[Z, z])+(Subscript[DKP, z]))/(Subscript[h, z]),(Subscript[Z, z])/(Subscript[h, z]),((Subscript[Z, sh])+(Subscript[DKP, sh]))/(Subscript[h, s]),(Subscript[Z, sh])/(Subscript[h, s]),((Subscript[Z, sč])+(Subscript[DKP, sč]))/(Subscript[h, s]),(Subscript[Z, sč])/(Subscript[h, s]),((Subscript[Z, š])+(Subscript[DKP, š]))/(Subscript[h, š]),(Subscript[Z, š])/(Subscript[h, š]),(Subscript[DKP, z])/(Subscript[h, z]),(Subscript[DKP, sh])/(Subscript[h, s]),(Subscript[DKP, sč])/(Subscript[h, s]),(Subscript[DKP, š])/(Subscript[h, š])},{(Subscript[TC, z]),30,65},PlotStyle->{Green,{Green,Dashed},Cyan,{Cyan,Dashed},Blue,{Blue,Dashed},Red,{Red,Dashed},{Green,DotDashed},{Cyan,DotDashed},{Blue,DotDashed},{Red,DotDashed}},
AxesLabel->{"Baseload cena (EUR/MWh)","Zisk/Ztráta (EUR/MWh)"},PlotLegend->{"Jaderná el. s FKP","Jaderná el. bez FKP","Uhelná el.- hnědé uhlí s FKP","Uhelná el.- hnědé uhlí bez FKP","Uhelná el.- černé uhlí s FKP","Uhelná el.- černé uhlí bez FKP","Paroplynová el. s FKP","Paroplynová el. bez FKP","Dynamické kapacitní platby - jaderná el.","Dynamické kapacitní platby - hnědouhelná el.","Dynamické kapacitní platby - černouhelná el.","Dynamické kapacitní platby - paroplynová el."}]

(*Graf - Vliv dynamických kapacitních plateb se stropem 50 000 EUR/MW na zisk/ztrátu elektráren*)
Plot[{(Subscript[Z, z])+(Subscript[ODKP, z]),(Subscript[Z, sh])+(Subscript[ODKP, sh]),(Subscript[Z, sč])+(Subscript[ODKP, sč]),(Subscript[Z, š])+(Subscript[ODKP, š]),(Subscript[ODKP, z]),(Subscript[ODKP, sh]),(Subscript[ODKP, sč]),(Subscript[ODKP, š])},{(Subscript[TC, z]),20,65},PlotStyle->{Green,Cyan,Blue,Red,{Green,DotDashed},{Cyan,DotDashed},{Blue,DotDashed},{Red,DotDashed}},
AxesLabel->{"Baseload cena (EUR/MWh)","Zisk/Ztráta (EUR/MW)"},PlotLegend->{"Jaderná el. s FKP","Uhelná el.- hnědé uhlí s FKP","Uhelná el.- černé uhlí s FKP","Paroplynová el. s FKP","Dynamické kapacitní platby - jaderná el.","Dynamické kapacitní platby - hnědouhelná el.","Dynamické kapacitní platby - černouhelná el.","Dynamické kapacitní platby - paroplynová el."}]