Zad 1

P=[2,3,1 ;
  4,1,2 ;
  3,5,1];

A=P;

[lw,lk]=size(A);

detA=A(1,1)*A(2,2)*A(3,3)+A(2,1)*A(3,2)*A(1,3)+A(3,1)*A(1,2)*A(2,3)-A(1,3)*A(2,2)*A(3,1)-A(
2,3)*A(3,2)*A(1,1)-A(3,3)*A(1,2)*A(2,1);
disp("Wyznacznik: ");
disp(detA);
X=det(P);
disp("Sprawdzenie wyznacznika: ");
disp(X);
if detA!=0
i=3;
j=3;
D=zeros(3,3);
M=P;

D=[M(2,2)*M(3,3)-M(2,3)*M(3,2),(-1*(M(2,1)*M(3,3)-M(2,3)*M(3,1))),M(2,1)*M(3,2)-M(2,2)*M(3,
1);

(-1*(M(1,2)*M(3,3)-M(1,3)*M(3,2))),M(1,1)*M(3,3)-M(1,3)*M(3,1),(-1*(M(1,1)*M(3,2)-M(1,2)*M(
3,1)));

M(1,2)*M(2,3)-M(1,3)*M(2,2),(-1*(M(1,1)*M(2,3)-M(1,3)*M(2,1))),M(1,1)*M(2,2)-M(1,2)*M(2,1)];
disp("macierz dopełnien: ");
disp(D);
DT=zeros(lw,lk);
for i=1:lw
for j=1:lk
DT(j,i)=D(i,j);
end
end
disp("Macierz dopełnień transponowana: ");
disp(DT);
P_odwrotne=(1/detA)*DT;
disp("Macierz odwrotna");
disp(P_odwrotne);
disp("Sprawdzenie macierzy odwrotnej: ");
disp(inv(P));
else
disp("Wyznacznik równy 0 -> brak macierzy odwrotnej");
end