first_lab

x=-10:10;
y=ones(1,21);
X=x'*y;
Y=y'*x;

R=max(abs(X),abs(Y));
% surf(X, Y, -R);
% title('Исходная пирамида');

% det(R);    % проверяем, что определитель равен 0

R = R + eye(21);
% surf(X, Y, -R);
% title('Прибавление к исходной матрице единичной');

% det(R);
  % складываем с единичной матрицей того же размера
% %det(R1);           % проверяем, что определитель теперь не равен 0
% 
% surf(X, Y, inv(R));

       % график обращения матрицы командой inv 
% surf(X, Y, ones(21)./R) 
% title('Матрица поэлементного деления матрицы ones(n,n) на R');% график поэлементного деления 
%                             % матрицы из единиц на R1
% 
% A = [1 2 3;
%      4 5 6;
%      7 8 9];
% 
% disp(expm(A));

% RSqrt = sqrtm(R);
% surf(X, Y, RSqrt);
% title('Результат операции sqrtm(R)');

% Rsqrt = sqrt(R1);
% % surf(X, Y, Rsqrt);
% RSqrt = sqrtm(R1);
% % surf(X, Y, RSqrt);
% Rlog = log(R1);
% % surf(X, Y, Rlog);
% RLog = logm(R1);
% % surf(X, Y, RLog);
% Rexp = exp(R1);
% % surf(X, Y, Rexp);
% RExp = expm(R1);
% % surf(X, Y, RExp);
% 
% R(11:21,1:11)=10;
% surf(X, Y, -R);
% title('"Отрежем" угол у пирамиды');

R1 = sqrt(X.^2 + Y.^2);
% surf(X, Y, -R1);
% title('Конус');

A = [1, 1, 1; 1, 2, 1; 1, 1, 1];
A = kron(A, R1);
A(22:42, 22:42) = A(22:42, 22:42) -15;
% surf(-A);
% title('Увеличенный в 2 раза центральный конус');

imagesc([22,30], [22,30], A);