function u=solexa(x,t)

L = 1000; t0=.1; x0=.5;
u = exp(-L*(t+t0)*(x-x0).^2);