exact_function_heat_equation.m

function res = exact_function_heat_equation(glob, params)

t = params.t;

%disp(['time: ', num2str(t)]);

k0 = params.diff_k0;

epsilon = pi;

X = glob(:,1);
Y = glob(:,2);

res = exp(-k0.* epsilon .^2 .* t) ...
       .* sin(1/sqrt(2).*epsilon .* X) ...
       .* cos(1./sqrt(2) .* epsilon .* Y);