n = 6;
h = 0.1;
k = h^2/n;

//считаем T
t(1) = 0;
for j = 2:n+1
    t(j) = t(j-1)+k;
end

//объявляем матрицу
u = zeros(n+1, n+1);

index = 1;
for x = 0:0.1:0.6
    u(1, index) = 1.2 + log(x + 0.4);
    index = index + 1;
end

index = 1;
for i = t(1):k:t(n+1)
    u(index, 1) = 0.8 + i;
    index = index + 1;
end

for i = 1:n+1
    u(i, n+1) = 1.2;
end

for i = 2:n+1
    for j = 2:n
        u(i, j) = 1/6 * (u(i-1, j-1) + 4*u(i-1, j) + u(i-1, j+1));
    end
end

mprintf("| j  |    i   |");
for i = 0:n
    mprintf("    %d   |", i);
end
mprintf("\n------------------------------------------------------------------------------\n");
mprintf("|    | tj\\xi  |");
for x = 0:0.1:0.6
    mprintf("   %1.1f  |", x);
end
mprintf("\n------------------------------------------------------------------------------\n");

for i = 1:n+1
    mprintf("| %d ", i-1);
    mprintf(" | %2.4f | ", t(i));
    for j = 1:n+1
        mprintf("%2.4f | ", u(i, j));
    end
mprintf("\n------------------------------------------------------------------------------\n");
end