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