function [w, K, F, yA, MH] = LineerSonluFarklar(a,b,alpha,beta,N)

% y'' = p(x)y' + q(x)y + r(x), y(a) = alpha, y(b) = beta

p = @(x) 2;
q = @(x) -1;
r = @(x) x*exp(x) - x;

h = (b-a)/N;

x = a:h:b;

% Sınır koşulları

w(1,1) = alpha;

w(N+1,1) = beta; 

% Katsayılar matrisi (K)

for j = 1:N-1
    
    K(j,j) = 4 + 2*(h^2)*q(x(j+1));
    
end

for j = 1:N-2
    
    K(j,j+1) = -2 + h*p(x(j+1));
    K(j+1,j) = -2 - h*p(x(j+2));
    
end

% Sağ taraf vektörü (F)

F(1,1) = -2*(h^2)*r(x(2)) - (-2-h*p(x(2)))*alpha;
F(N-1,1) = -2*(h^2)*r(x(N-1)) - (-2+h*p(x(N)))*beta;

for j = 2:N-2
    
    F(j,1) = -2*(h^2)*r(x(j+1));
    
end


% Bilinmeyen vektörü

w(2:N) = K\F;    % inv(K)*F

% Analitik çözüm

yA = (1/6)*(x.^3).*exp(x) - (5/3)*x.*exp(x) + 2*exp(x) - x - 2;
    
MH = abs(yA' - w);