function f=fun_f(y,x)
%uwaga ! zamieniiona kolejnosc zmiennyv
global a d
f=x.*(1-x)-a*x.*y./(x+d);
