function [x,u1,u2,v1,v2,w,yA,MutlakHata] = lineeratis(a,b,alpha,beta,N)

% y'' = p(x)y' + q(x)y + r(x)

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

% Adım uzunluğu

h = (b-a)/N; 

x = a:h:b;

u1(1) = alpha;
u2(1) = 0;
v1(1) = 0;
v2(1) = 1; 

for i = 1:N
    
    k1 = h*u2(i);
    m1 = h*(p(x(i))*u2(i) + q(x(i))*u1(i) + r(x(i)));
    k2 = h*(u2(i) + m1/2);
    m2 = h*(p(x(i)+h/2)*(u2(i)+m1/2) + q(x(i)+h/2)*(u1(i)+k1/2) + r(x(i)+h/2));
    k3 = h*(u2(i) + m2/2);
    m3 = h*(p(x(i)+h/2)*(u2(i)+m2/2) + q(x(i)+h/2)*(u1(i)+k2/2) + r(x(i)+h/2));
    k4 = h*(u2(i)+m3);
    m4 = h*(p(x(i)+h)*(u2(i)+m3) + q(x(i)+h)*(u1(i)+k3) + r(x(i)+h));
    
    kk1 = h*v2(i);
    mm1 = h*(p(x(i))*v2(i) + q(x(i))*v1(i));
    kk2 = h*(v2(i)+mm1/2);
    mm2 = h*(p(x(i)+h/2)*(v2(i)+mm1/2) + q(x(i)+h/2)*(v1(i)+kk1/2));
    kk3 = h*(v2(i)+mm2/2);
    mm3 = h*(p(x(i)+h/2)*(v2(i)+mm2/2) + q(x(i)+h/2)*(v1(i)+kk2/2));
    kk4 = h*(v2(i) + mm3);
    mm4 = h*(p(x(i)+h)*(v2(i)+mm3) + q(x(i)+h)*(v1(i)+kk3));
    
    u1(i+1) = u1(i) + (1/6)*(k1 + 2*k2 + 2*k3 + k4);
    u2(i+1) = u2(i) + (1/6)*(m1 + 2*m2 + 2*m3 + m4);
    v1(i+1) = v1(i) + (1/6)*(kk1 + 2*kk2 + 2*kk3 + kk4);
    v2(i+1) = v2(i) + (1/6)*(mm1 + 2*mm2 + 2*mm3 + mm4);
    
end

for i = 1:N+1
    
    w(i) = u1(i) + ((beta - u1(N+1))/v1(N+1))*v1(i);
    
end

%Analitik çözüm

yA = ...; 

% Hata 

MutlakHata = abs(yA - w);