function y = quad1(x,p) 
y = p.*x.*(1-x);