% jarkovskeho_jev_phi_th.m
% Jarkovskeho jev 1-D: uhel tepelneho zpozdeni.
% Miroslav Broz (miroslav.broz@email.cz), Jan 28th 2006

function y = phi_th(f,T_eq,K)

y = atan2(- Theta(f,T_eq,K), (1. + Theta(f,T_eq,K) ));

return

