%%
%% Listing 7.14 w książce "Obliczenia inżynierskie i naukowe", PWN 2011
%% (C) Piotr Krzyżanowski, 2011
%%
%% Możesz używać tego kodu zgodnie z licencją Creative Commons BY-SA, 
%% zob. http://creativecommons.org/licenses/by-sa/3.0/pl/
%%

function L = lap1ddset(N, h)
% tworzy macierz rzadką, odpowiadajacą macierzy 
% (minus!) jednowymiarowego laplasjanu, z zerowym warunkiem Dirichleta
% N - liczba węzłów wewnętrznych, h - krok siatki

ii = [1:N, 2:N,   1:N-1];
jj = [1:N, 1:N-1, 2:N];
vv = [2*ones(1,N), -ones(1,2*(N-1))]/(h^2);
L = sparse(ii,jj,vv,N,N);
end
