if detA!=0 i=3;
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