function y = quad4(x,a) 
y = x.*(1+a.*(1-x));