contour plot e(1-1/(1-(x-t)^2))*unitstep(1-(x-t)^2)