function[u2o] = pois1d_2o(N)

h = 1/N;

%
% Randbedingungen
%
g0 = cos(0);
g1 = cos(1);

%
% Aufsetzen der Matrix und der rechten Seite des LGS
%
A = sparse(zeros(N-1,N-1));
b = sparse(zeros(N-1,1));

for i = 1:N-1
    A(i,i) = 2;
    f = cos(i*h);
    b(i) = h*h*f;
end

for i = 1:N-2
    A(i,i+1) = -1;
    A(i+1,i) = -1;
end

b(1) = b(1) + g0;
b(N-1) = b(N-1) + g1;

%
% Lösung des LGS 
%
u2o = A\b;