function xdot=vx(x,t)

if x>0
xdot=-x*log(x);
else
    xdot=0
end