function y=interfun(x)
global c
y=c*x.*(1-x);
