%% Löser ett linjärt randvärdesproblem med Dirichletvillkor
function [x,y] = rvp2(N)

a = 0;
b = 1;
h = (a+b)/N;    % intervallstorlek
x = a:h:b;      % same as x = linspace(a,b,N+1);

% bilda matrisen
d = [1; (-2-4*h^2)*ones(N-1,1); 1];	% diagonal
ud = [0; 0; ones(N-1,1)];
ld = [1; ones(N-2,1); 0; 0];
D = [ld d ud];
A = spdiags(D,[-1 0 1],N+1,N+1);

% högerled
b = [1; -50*h^2*cos(4*pi*x(2:end-1)'); 3];

% löser systemet
y = A\b;


