1 %Dane
2 P=[-3, -2, 3; 1, -1, 5; -8, -4, 4];
3 A=P;
4
5 [lw, lk]=size(A);
6
7 detA=A(1,1)*A(2,2)*A(3,3) + A(1,2)*A(2,3)*A(3,1) + A(1,3)*A(2,1)*A(3,2) -
8 A(1,3)*A(2,2)*A(3,1) - A(1,1)*A(2,3)*A(3,2) - A(1,2)*A(2,1)*A(3,3);
9
10 disp("Wyznacznik: ");
11 disp(detA);
12 X=det(P);
13 disp("Sprawdzenie wyznacznika: ");
14 disp(X);
15 if detA!=0
16 i=3;
17 j=3;
18 D=zeros(3,3);
19 M=P;
20 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)
21 *M(3,2)-M(2,2)*M(3,1);
22 (-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)
23 *M(3,2)-M(1,2)*M(3,1)));
24 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)
25 -M(1,2)*M(2,1)];
26 disp("macierz dopełnien: ");
27 disp(D);
28 DT=zeros(lw,lk);
29 for i=1:lw
30 for j=1:lk
31 DT(j,i)=D(i,j);
32 end
33 end
34 disp("Macierz dopełnień transponowana: ");
35 disp(DT);
36 P_odwrotne=(1/detA)*DT;
37 disp("Macierz odwrotna");
38 disp(P_odwrotne);
39 disp("Sprawdzenie macierzy odwrotnej: ");
40 disp(inv(P));
41 else
42 disp("Wyznacznik równy 0 -> brak macierzy odwrotnej");
43 end