% T_x_t.m
% Jarkovskeho jev 1-D: zmeny teploty v hloubce a v case.
% Miroslav Broz (miroslav.broz@email.cz), Jan 28th 2006

function y = T(x,t)

global T_eq T_1_0_ f phi_th_ chi;

y = T_eq + T_1_0_ .* cos(2.*pi.*f .* t + phi_th_ - sqrt(pi.*f./chi) .* x) .* exp(-sqrt(pi.*f./chi) .* x);

return

