function [w, yA, MH, lambda] = hiperbolikKDD(L,T,m,n,alpha)

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

fx0 = @(x,t) 0; 
fxL = @(x,t) 0;
ft0 = @(x,t) sin(pi*x);
ft1 = @(x,t) 0;

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));
    w(m+1,j)=fxL(x(m+1),t(j));
end

%IC'S

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

for i = 2:m
    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);
    end
end

% Analitik çözüm

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

yA = sin(pi*XX).*cos(2*pi*TT);
    
MH = abs(yA - w);