Remove["Global`*"]

(*Simulera population av fisk för ett visst antal fiskare och ett viss kvot i varje tidssteg*)
x0=0.5;
TMax=100;

Manipulate[Plot[Evaluate[x[t]/.NDSolve[{x'[t]==x[t]*(1-x[t])-n*k*x[t],x[0]==x0},x,{t,0,TMax}]],{t,0,TMax},PlotRange->{{0,TMax},{0,1}}],{k,0,1},{n,1,100}]

(*2 fiskare och kvot av varje tidssteg (de vet hur mkt fisk det är i varje tidssteg)*)
n=2;
TMax=100;
x0=0.8;

Manipulate[Plot[Evaluate[x[t]/.NDSolve[{x'[t]==x[t]*(1-x[t])-k1*x[t]-k2*x[t],x[0]==x0},x,{t,0,TMax}]],{t,0,TMax},PlotRange->{{0,TMax},{0,1}}],{k1,0,1},{k2,0,1}]

(*Optimal andel är (1-x0)/n*, då är populationen stabil*)

(*3 fiskare som tar ut en viss del av x0 i varje tidssteg*)
Needs["PlotLegends`"]
TMax=30;
x0=1;
Manipulate[Plot[Evaluate[{x[t],y1[t],y2[t],y3[t],y1[t]+y2[t]+y3[t]}/.NDSolve[{x'[t]==x[t]*(1-x[t])-(k1+k2+k3)*x0,y1'[t]==k1*x[t],y2'[t]==k2*x[t],y3'[t]==k3*x[t],y1[0]==0,y2[0]==0,y3[0]==0,x[0]==x0},{x,y1,y2,y3},{t,0,TMax}]],{t,0,TMax},PlotRange->{{0,TMax},{0,10}}],{k1,0,1},{k2,0,1},{k3,0,1}]