function [x,y,K,F,MutlakHata,yA] = lineersonlufarklar(a,b,alpha,beta,N)

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

p = @(x) ...;
q = @(x) ...;
r = @(x) ...; 

h = (b-a)/N;

x = a:h:b; 

y(1) = alpha;
y(N+1) = beta; 

% Katsayılar matrisi

K = zeros(N-1);

%Köşegen elemanlar
for i = 1:N-1
    
    K(i,i) = -4 -2*(h^2)*q(x(i+1));
    
end

% alt ve üst köşegen elemanlar

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

% Sağtaraf vektörü F

F = zeros(N-1,1);

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

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

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

% Sayısal çözüm

w = K\F; % w = inv(K)*F

for i = 2:N
    
    y(i) = w(i-1);
    
end

% Analitik çözüm

yA = ...; 
    
MutlakHata = abs(yA - y);