function [w, lambda, TT,XX] = hiperbolikKDDFSoru7b(L,T,m,n,alpha)

% u_tt - alpha^2 u_xx = f(x,t), u(0,t)=fx0, du/dx(L,t)=fx1, 
% u(x,0)=ft0, du/dt(x,0)=ft1

f = @(x,t) ...; 
fx0 = @(x,t) ...; 
fx1 = @(x,t) ...;
ft0 = @(x,t) ...;
ft1 = @(x,t) ...;

h = L/m;
k = T/n;

x = 0:h:L;
t = 0:k:T;

lambda = (alpha*k)/h;
mu = 2*(1-lambda^2);

% BC'S

for j=1:n+1
    w(1,j) = fx0(x(1),t(j));
end

%IC'S

for i = 1:m+1
    w(i,1) = ft0(x(i),t(1));
end

for i = 2:m+1
    w(i,2) = w(i,1) + k*ft1(x(i),t(1));
end

% Sayısal çözümler

for j = 2:n
    for i = 2:m
        w(i,j+1) = (lambda^2)*w(i+1,j) + mu*w(i,j) + (lambda^2)*w(i-1,j) - w(i,j-1) + (k^2)*f(x(i),t(j));
    end
    w(m+1,j+1) = w(m,j+1) + h*fx1(x(m+1),t(j+1)); 
end

% Analitik çözüm

[TT,XX] = meshgrid(t,x);

yA = ...;
    
MH = abs(yA - w);