function g=fun_g(v,u)
%uwaga ! zamieniiona kolejno zmiennych
global b

g=b*v.*(1-v./u);