function z=twosol(x,t)
p=2;
q=1.5;
theta=p*x-4*p^3*(t-0.5);   
chi=q*x-4*q^3*(t-0.5);     
z=2*(p^2-q^2)*((p^2+q^2*(sech(chi)).^2.*(sinh(theta)).^2)./ ... 
    ((p*cosh(theta)-q*tanh(chi).*sinh(theta)).^2));