function u = ukol2(k,l,n)
% u = ukol2(k,l,n)
%
% Medota siti pro danou 1D Dirichletovu ulohu 
%
% k   ... tuhost struny
% l   ... delka struny
% n+1 ... pocet pravidelnych uzlu na strune
%
% u   ... vysledny pruhyb struny 


% vzdalenost mezi sousednimi uzly
h = l / n;


% vytvoreni matice tuhosti K
K = diag(2*ones(n-1,1)) + diag(-ones(n-2,1),1) + diag(-ones(n-2,1),-1);
K = k/h^2 * K;


% vektor prave strany f
% doplnte kod
%
% doplnte kod 


% reseni soustavy u
u = K \ f;
u = [0; u; 0];


% vykresleni reseni
figure; hold on
x = 0 : h : l;
plot(x,u)

% konec

