 
function  f=xdoth(x)

global a h

f=a*x.*(1-x)-h;
%f=sin(pi*x);
