function s2=ff(x, t)

      ss0=8.0*(x^2-1);
      ss1=ss0*(x+t);
      ss2=8+16.*x*(x+t)+ss0+2.*(x+t)*(8.*x+ss0*(x+t));
      s2=(ss1-ss2)*exp((x+t)^2);