plot 1-exp(t^2/(t^2-1))*unitstep(1-t) from t=0 to 1.5