function main()
%Zeithorizont
T=??;
%Wieviele Zeitintervalle fuer die Steuerungsparametrisierung 
diss=??;
%Startwerte
x0=[??;??;??];
u=zeros(diss,1);
%Untere und obere Schranken fuer die Steuerungen
lb=??;
ub=??;
options = optimset('Display','iter','TolX',1e-3,'MaxFunEvals',1e6);
%Wonach wird optimiert?
[x,fval,exitflag,output] = fminsearchbnd(@(??)obj_fun(x0,u,T,diss),u,lb,ub,options);
t=[0:T/diss:T];
figure
stairs(t',[x;x(length(x))])
title('Control')
xlabel('t')
ylabel('Fishing')


%Berechnung der Zielfunktion. 
%Loesungsalgorithmus:
%1)Loese die ODE im ersten Intervall 
%2)Veraendere die Steuerung und die Startwerte und loese die ODE auf dem folgenden Intervall
%3) usw.
function L=obj_fun(x0,u,T,diss)
hilf_mir_t=[];
hilf_mir_x=[];
%Berechne Loesung auf dem ersten Intervall
[t,y] = ode45(@(t,y)lot_vol(t,y,u(??)),[0,??],x0);
%Abspeichern
hilf_mir_t=[hilf_mir_t t'];
clear t;
hilf_mir_x=[hilf_mir_x y'];
clear y;
k=2;
%Schleife fuer die restlichen Intervalle.
for i=??:??:??
    [t,y] = ode45(@(t,y)lot_vol(t,y,u(??)),[i,i+??],[hilf_mir_x(:,??)]);
    hilf_mir_t=[hilf_mir_t t'];
    clear t;
    hilf_mir_x=[hilf_mir_x y'];
    clear y;
    k=k+1;
end
%Berechnung des Zielfunktionals
L=??
%Plotting;
plot(hilf_mir_t',hilf_mir_x(1:2,:)')
title('Species')
xlabel('t')
ylabel('Species')
pause(0.005)
clear hilf_mir_x;
clear hilf_mir_t;

%ODE Achtung: yp(1)=erste Variable, yp(2)=zweite Variable, yp(3)=Zielfunktional
function yp = lot_vol(t,y,u)
c1=1.;c2=1.;
yp = [??;??;??];