Well, it's sort of better now. It does SOMETHING for contour(T, Y, f(T, Y), 30, 'k'), but then errors out saying :
Warning: Divide by zero.
> In @(t, y) (-t)./(exp(-y).*cos(y)-exp(-y))
??? Error using ==> contourc
Contour levels must be finite.
Warning: Error occurred while evaluating listener...