%openf93, randvärdesproblem u''+2u'+5u=0, u(0)=1,u(1)=0
x0=0;xslut=1;
n=19;h=(xslut-x0)/(n+1);
A=zeros(n,n);b=zeros(n,1);
for i=1:n-1
    A(i,i)=-2/h^2+5;
    A(i+1,i)=1/h^2-2/(2*h);
    A(i,i+1)=1/h^2+2/(2*h);
end
A(n,n)=-2/h^2+5;
b(1)=2/(2*h)-1/h^2;
u=A\b
up=[1;u;0];
xp=[0:h:1]';
plot(xp,up,'r')