function dx = valentinefunction(t,x)
global V K  

dx = zeros(size(x));

dx(1) = -(V*x(1)/(x(1)+K));
dx(2) = -(dx(1)) ;

end