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

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

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

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

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

lambda = ((alpha^2)*k)/(h^2);
mu1 = 1+lambda;
mu2 = 1-lambda;

% 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

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

% Katsayılar matrisi

for i = 1:m-1
    K(i,i) = mu1;
end

for i = 1:m-2
    K(i,i+1) = -lambda/2;
    K(i+1,i) = -lambda/2;
end

% Denklem sistemi çözümleri

for j = 1:n
    
    F(1,1) = (lambda/2)*w(1,j) + mu2*w(2,j) + (lambda/2)*w(3,j) + (lambda/2)*w(1,j+1) + (k/2)*(f(x(2),t(j)) + f(x(2),t(j+1)));
    F(m-1,1) = (lambda/2)*w(m-1,j) + mu2*w(m,j) + (lambda/2)*w(m+1,j) + (lambda/2)*w(m+1,j+1) + (k/2)*(f(x(m),t(j)) + f(x(m),t(j+1)));
    
    for i = 2:m-2
        F(i,1) = (lambda/2)*w(i,j) + mu2*w(i+1,j) + (lambda/2)*w(i+2,j) + (k/2)*(f(x(i+1),t(j)) + f(x(i+1),t(j+1)));
    end
    
    ww = K\F;
    
    for i = 2:m
        w(i,j+1) = ww(i-1,1);
    end
    
end

% Analitik çözüm

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

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