clear
clc
Nx = 5;
Ny = 6;
%Dieu kien sai so
tol = 1d-6;
%Sai so
err =1;
k = 0;
T = [ 0 300 300 300 0; 75 0 0 0 100; 75 0 0 0 100; 75 0 0 0 100; 75 0 0 0 100; 0 50 50 50 0 ];
Tkp1 = T;
fprintf( ' k|')
for j = 1 : Ny -2
    for i = 1 : Nx -2
        fprintf(' T(%li,%li) |',j,i)
    end
end
fprintf(' error\n')
while err > tol
    k = k + 1;
    fprintf('%4i |',k);
    % Tính T(i,j) theo phuong pháp sai phân trung tâm
    for j = 2: Ny -1
        for i = 2 : Nx -1
            Tkp1(j,i) = 0.25*(T(j+1,i) + T(j,i+1) + T(j-1,i) + T(j,i-1));
            fprintf( '%10.6f |', Tkp1(j,i));
        end
    end
    %Tính sai so
    err = sqrt(sum(sum(Tkp1 - T).^2));
    fprintf( '%10.6f \n', err);
    T =Tkp1;
    surf(T);
    getframe();
end