function [y]=bv1(q,f,bound,x);
% Lös randvärdesproblemet
% y''-q(x) y = f,  y(x(0))=bound(1),  y(x(N+1))=bound(2)
% vektorn x(1:N)=x(0)+h:h:x(N+1)-h;
% q antas vara ickenegativ

N=length(x); h=x(2)-x(1);
qv=feval(q,x); fv=feval(f,x); fv=fv';

% Konstruera tridiagonal matris
A=-2*eye(N)+diag(ones(N-1,1),1)+diag(ones(N-1,1),-1);
A=A-h^2*diag(qv)

% Högerled
b=h^2*fv;
b(1)=b(1)-bound(1);
b(N)=b(N)-bound(2);

y=A\b;

