% Theta.m
% Jarkovskeho jev 1-D: tepelny parametr Theta.
% Miroslav Broz (miroslav.broz@email.cz), Jan 28th 2006

function y = Theta(f,T_eq,K)

global C rho eps sigma;

y = sqrt(pi .* f*K*C*rho) ./ (4.*pi*eps*sigma * T_eq .^ 3.);

return

