load bbear12
sigma1 = 1;
sigma2 = 1;
n = 20;
Z=twobumpy(x,y,xdata,ydata,tdata,sigma1,sigma2,n);
who
mesh(x,y,Z);
