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

% u_t - alpha^2 u_xx = 0, u(0,t)=fx0, u(L,t)=fxL, u(x,0)=ft0

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

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

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

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

%Sınır koşulları
for j = 1:n+1
    
    w(1,j) = fx0(x(1),t(j));
    w(m+1,j) = fxL(x(m+1),t(j));
    
end

%Başlangıç koşulu
for i = 1:m+1
    
    w(i,1) = ft0(x(i),t(1));
    
end

% Bilinmeyenlerin elde edilmesi

for j = 1:n
    for i = 2:m
        
        w(i,j+1) = mu*w(i,j) + lambda*(w(i+1,j) + w(i-1,j));
        
    end
end

% Analitik çözüm

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

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