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);