%lab3%
clear all
clc
close all
%--------------------------------------------------------------------------
a = 1;
h = 0.1;
l = 0.1;
n=10;
m=10;
i = 0:n;
j = 0:m;
x = i*h;
y = j*l;
u = zeros(n+1,m+1);
%-------------------
u( 0+1, j +1) = y(j+1); %mu1
u( n+1, j +1) = y(:) - 1 + y(:)*exp(y(:));%mu2
u( i+1, 0 +1) = x(:) - x(:).^3;
u( i+1, m +1) = x(:) + 0.11;

um = sum( u( 0+1, j +1) );
um = um + sum( u( n+1, j +1) );
um = um + sum( u( i+1, 0 +1) );
um = um + sum( u( i+1, m +1) );
um=um/40;
%
for ind2 = 2:(m)
    for ind1 = 2:n
        u( ind1, ind2 ) = um;
    end
end

for ind2 = 2:(m)
    for ind1 = 2:n
        u( ind1, ind2 ) =  (  u( ind1+1  , ind2 ) ...
                            + u( ind1-1  , ind2 ) ...
                            + u( ind1,     ind2+1 ) ...
                            + u( ind1,     ind2-1 ))/4;

    end
end
%вывод результатов: - - - - - - - - - - - - - - - - - - - - - - - - - - - -
figure('name','fig1');surf(y,x,u);xlabel('x');ylabel('y');
disp('Значения u:');
disp( u );
%--------------------------------------------------------------------------