%% Initialwerte definieren
N=length(inFlow);
%% Restriktionen definieren
constraints;
%% zu maximierende Funktion definieren
P=c.*[price;-price;zeros(N,1);zeros(N,1);0;0];
%% optimierung
[x1,fval1,exitflag,output,lambda]=linprog(P,[],[],Aeq,Beq,Lb,Ub);
%% visualisierung
plotResults(x1,price,N);
%plotResults(x1,price,N);