function eq = defHCEQ(filename)
% DEFHCEQ Funkce, kde jsou definována data přímé úlohy vedení tepla
% Vstupy: filename - cesta k soubotu, kde je uložena triangulace z PDE
% Modelleru
% Výstupy: struktura
% Koeficienty rovnice
eq.c = @(x,y,t,T) 0;
eq.p = @(x,y,t) 0;
eq.q = @(x,y,t,T) 0.0;
eq.g = @(x,y,t,T) 0.0;
% Načtení triangulace oblasti
if ~exist('filename','var')
    filename = 'regions/wall';
end
load(strcat(filename,'.mat')); % Soubor s exportovanou sítí z PDE Modelleru
eq.domain.p = p; % Points: Uzly
eq.domain.e = e; % Edges: Hrany
eq.domain.t = t; % Triangles: Prvky
% Časový interval
eq.time = []; % 0, až konečný čas
eq.num_tsteps = 0; % Počet časových bodů
% Okrajová podmínka
eq.f = {@(x,y,t) 0,@(x,y,t) 0,@(x,y,t) 0,@(x,y,t) 0};
eq.alpha = {@(x,y,t,T) 0,@(x,y,t,T) 0,@(x,y,t,T) 0,@(x,y,t,T) 0};
% Testovací tepelné toky
eq.beta = {@(x,y,t,T) 0,@(x,y,t,T) 0,@(x,y,t,T) 0,@(x,y,t,T) 0};
eq.dir_edge = []; % Indexy hran v eq.domain.e, kde je předepsána Dirichletova OP, na zbytku je Neumann
eq.neum_edge = []; % Indexy hran v edges, kde je předepsána Neumannova OP
% Počáteční podmínka
eq.phi = @(x,y) 0;
% Typ úlohy
eq.stationary = true;
eq.linear = true;
end

