clear all, close all, clc
%Losning av randvardesproblem

N=10;   %Antal delintervall
n=N-1;  %Antal obekanta
h=1/N;
x=h*(1:n)'; %x-koordinaterna dar vi soker losningen

A=diag(1/h^2*ones(n-1,1),-1)...
    +diag(-2/h^2-100*x.^2)+diag(1/h^2*ones(n-1,1),1);
b=zeros(n,1);
b(1)=-1/h^2;
b(end)=-2/h^2;
y=A\b;
plot([0;x;1],[1;y;2])
xlabel('x')
ylabel('y')
title(['Losning av randvardesproblem med ', num2str(N),' delintervall.'])