function [A, f] = LaplaceMatVek1D(J)
%LAPLACEMATVEK1D Baut Matrix und rechte Seite einer Laplace-Gleichung auf
% Mittels 1D-Hutfunktionen auf einem dyadischen Gitter wird die Laplace-Gleichung
%   -Laplace y = 1  in (0, 1)
%            y = 0  auf Rand
% geloest. Dies fuehrt auf ein lineares Gleichungssystem A*y = f, von dem
% die Matrix A und die rechte Seite f zurueckgeliefert werden.
%
% Eingabe:
%   - J: Level der Verfeinerung (2^J - 1 Elemente)
%
% Ausgabe:
%   - A: Matrix des linearen Gleichungssystems A*y = f, Gram-Matrix der
%        Bilinearform a(y, v) = \int_\Omega (\nabla y) (\nabla v) d\Omega
%   - f: Rechte Seite b des linearen Gleichungssystems
%
% Copyright (C) 2017 Samuel Leweke (Universitaet zu Koeln)

	N = 2^J - 1; %Anz. Elemente

	A = 2^(J) * spdiags([-ones(N, 1), 2 * ones(N, 1), -ones(N, 1)], -1:1, N, N);
	f = 2^(-J) .* ones(size(A, 1), 1);
end