clear;

N=100; % Rozmer matice T_N
h=1/(N+1);

F = @(x) x.^3; % Prava strana Poissonovej rovnice

% Diskretizacia pravej strany Poissonovej rovnice
Delenie = linspace(0,1,N+2);
f = F(Delenie(2:N+1))';

% Matica T_N
T_N = 2*eye(N) + diag(-ones(N-1,1),-1)+diag(-ones(N-1,1),1);

% Implementujte Jacobiho metodu pre system T_n v=h^2f
% Vysledne riesenie ulozte do premennej X ako stlpcovy vektor

% ----- Vas program -------


% ------Zobrazenie grafov----
% Zobrazenie presneho riesenia
PresneRiesenie = @(x) -x.^5/20+1/20*x;
figure('Name','Presne riesenie');
plot(Delenie(2:N+1),PresneRiesenie(Delenie(2:N+1)));

% Zobrazenie Jacobiho riesenia
figure('Name','Riesenie Jacobiho metodou');
plot(Delenie(2:N+1),X');

