function z=func(x,t,zeta,x0)
z=2*zeta^2*(sech(zeta*(x-4*zeta^2*t+x0))).^2;
