function F=optimierung(parameter,S)

mue_max=parameter(1);

K_m=parameter(2);

K_i=parameter(3);

F=1/K_m*(S./(mue_max./K_m + S+(S.^2/K_i)));