function pinco=f(x,t)
%pinco=10*exp(-t)*(x.^2.*(1-x).^2-2*(1-6*x+6*x.^2));
pinco=zeros(size(x));
