 
function  f=xdot(x)

global a b

f=a*x.*(b-x);
