a=8.898; (*součinitel a,viz BP str.40*)

b15=((1.5/a)*((((x/a)))^(1.5-1)))*(E^((-((((x/a)))^1.5)))); (*součinitel b, viz BP str.41*)
b225=((2.25/a)*((((x/a)))^(2.25-1)))*(E^((-((((x/a)))^2.25))));
b3=((3/a)*((((x/a)))^(3-1)))*(E^((-((((x/a)))^3))));
graph=Plot[{b15,b225,b3},{x,0.0,25},PlotLabel->"Weibullovo rozdělení",LabelStyle->{Directive[Bold,22]},PlotStyle->{{Thick,Red},{Thick,Green},{Thick,Dashed,Blue}},AxesLabel->{"rychlost větru (m/s)","poměrná četnost rychlostí větru (-)"},AxesStyle->{Directive[Black,18],Black},TargetUnits->Automatic,PlotLegends->LineLegend[{"1,50","2,25","3,00"},LegendFunction->Frame]]

Integrate[b15,{x,0,25}](*integrace funkci grafu*)

Integrate[b225,{x,0,25}]

Integrate[b3,{x,0,25}]

rocniprumteplota =10; (*průměrná teplota Moldavy v Krušných horách*)
baromtlakmore=1013*100; (*hPa*)
r=287; (*plynová konstanta [J/kg*K]*)
vyskaVtE=850; (*nadmořská výška uvažované VtE v m.n.m.*)
As=(71^2); (*plocha protékaného vzduchu skrz VtE v m2*)
T=273.15+rocniprumteplota; (*přepočet na termodynamickou teplou (K)*)
prepocetbaromtlaku=baromtlakmore*(1-0.118358*(10^(-3))*vyskaVtE+0.5291*(10^(-8))*(((vyskaVtE))^2)); (*přepočet barometrického tlaku na danou nadmořskou výšku a teplotu*)
hustotafinal=(prepocetbaromtlaku/(r*T))(*přepočtená hustota vzduchu*)

energie={5390.33,5169.56}; (*vypočtená suma energie v MWh*)

vykup=1930*energie (*penezni zisk z vykupni casti*)

bonus=1340*energie (*penezni zisk z bonusu*)

0.02*17626379 (*vlastni spotreba elektrarny*)

S=30*(10^6)*((((1.1^10)*0.1)/((1.1^10)-1)))

10*% (*složene uroceni, rocni splatna castka za uver*)

(0.04+0.9*0.1)*((9100000/39100000))+(0.04+0.07)*(1-0.88)*((29000000/39100000))(*výpočet WAAC, nakonec navýšen na 9 %, tato hodnota se týká spíše státních zakázek, ne soukromého sektoru*)

-39100000+((∑_(t=1)^10)(((-5230000+17626400))/(((1+0.09))^t)))+((∑_(t=11)^20)(((-350000+17626400))/(((1+0.09))^t))) (*výpočet NPV*)

Solve[(((∑_(t=1)^x)(17.6264*(10^6)-5.23*(10^6))*(((1+0.09))^(-t)))-39.1*(10^6))==0,x](*výpočet Reálné doby návratnosti*)