%% Konstanten
% Pumpeturbine Wirkungsgrad
w=0.95;
% Umrechnungsfaktor m³-->MWh
c=998*(9.81)*(180.5)*(0.95)/(3600*1000000); % (kg/m³).(m/s²).m /(3600.1000000)
%% Systemgrenze Bounds
% [Turbine;Pumpe;Menge_max oberbecken;Menge_max_unterbecken;
Ub=[40*ones(N,1);40*ones(N,1);408000*ones(N,1);488000*ones(N,1);300000;400000];
Lb=[zeros(N,1);zeros(N,1);15885*ones(N,1);95885*ones(N,1);50000;100000];
%% Restriktion Aeq.x=Beq
t=ones(N,1);
z=zeros(N,1);
q=zeros(N,N);

Tr=-spdiags([t -t],[0 N],N,N);
Pp=w*spdiags([t -t],[0 N],N,N);
Fzob=spdiags([t -t],[0 -1],N,N);
Fzobx=q*Fzob;
Fzub=spdiags([t -t],[0 -1],N,N);
Fzubx=q*Fzub;
Fzob0=z;Fzob0(1,1)=1;
Fzub0=z;Fzub0(1,1)=1;
% 2 Gleihungen
Aeq=[Tr Pp Fzob Fzubx Fzob0 z;-Tr -Pp Fzobx Fzub z Fzub0];
Beq=[inFlow;z];