% moje testovaci uloha cislo 1
%L_x = 0; R_x = 4; T = 5; % vymezeni oblasti vypoctu
%a_advec = 0.3; % parametr advekce
%RHS = @(x,t) 0 * x + 0 * t; % prava strana rovnice
%c = 25; IC = @(x) exp(-c * (x - (L_x + R_x) / 2) .^ 2); % pocatecni podminka
%BC = @(t) 0 * t; % okrajova podminka
%
%[P,E] = MeshGeneration1D(L_x, R_x, 16); % moje vlastni generovani site pro 1D ulohy
%
%axis_limits = [L_x, R_x, -0.1, 1.1]; % pro mou testovaci ulohu


% moje testovaci uloha cislo 2
L_x = -5; R_x = 15; T = 10; % vymezeni oblasti vypoctu
a_advec = 1; % parametr advekce
RHS = @(x,t) 0 * x + 0 * t; % prava strana rovnice
b = 1; IC = @(x) exp((-log(2) / (b ^ 2)) * (x .^ 2)); % pocatecni podminka
BC = @(x,t) IC(x - a_advec * t); % okrajova podminka

[P,E] = MeshGeneration1D(L_x, R_x, 20); % moje vlastni bgenerovani site pro 1D ulohy

axis_limits = [L_x, R_x, -0.1, 1.1]; % pro mou testovaci ulohu


% testovaci uloha od Bajka
%L_x = -50; R_x = 150; T = 100; % vymezeni oblasti vypoctu
%a_advec = 1; % parametr advekce
%RHS = @(x,t) 0 * x + 0 * t; % prava strana rovnice
%IC = @(x) (2 + cos((pi / 2) * x)) .* exp(-log(2) * ((x / 10) .^ 2)); % pocatecni podminka
%BC = @(x,t) IC(x - a_advec * t); % okrajova podminka
%
%[P,E] = MeshGeneration1D(L_x,R_x,800); % moje vlastni generovani site pro 1D ulohy
%
%axis_limits = [L_x, R_x, -0.25, 3.25]; % pro mou testovaci ulohu


% uloha z knihy Hesthaven
% L_x = 0; R_x = 2 * pi; T = 10; % vymezeni oblasti vypoctu
% a_advec = 2 * pi; % parametr advekce
% RHS = @(x,t) 0 * x + 0 * t; % prava strana rovnice
% IC = @(x) sin(x); % pocatecni podminka
% BC = @(x,t) sin(x - a_advec * t); % okrajova podminka
% 
% [P,E] = MeshGeneration1D(L_x, R_x, 32); % moje vlastni generovani site pro 1D ulohy
% 
% axis_limits = [L_x, R_x, -1.1, 1.1]; % pro ulohu z knihy