/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
/* [ Created with wxMaxima version 12.04.0 ] */

/* [wxMaxima: comment start ]
Spusťte Ctrl + R;
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
kill(all);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
A: matrix(
 [1,2,0,0], 
 [2,1,1,3], 
 [0,1,0,1], 
 [0,4,7,1]
);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
Nachystání matice s nulovými prvky. Sem budeme vkládat poměry determinantů submatice a matice.
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
inverzniA: matrix(
 [0,0,0,0], 
 [0,0,0,0], 
 [0,0,0,0], 
 [0,0,0,0]
);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
Zjištění determinantu pomocí příkazu determinant(%).
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
detA:determinant(A);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
Zjištění submatice A(i,j) vynecháním i-tého řádku a j-tého sloupce pomocí příkazu submatrix (i,A,j).
Výpočet determinantu submatice.
Výpočet prvku (j,i) inverzní matice inverzniA jako součin (-1)^(i+j) s podílem determinantu submatice A(i,j) a determinantu matice A
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
for i:1 thru 4 do{
for j:1 thru 4 do{
detAs:determinant(submatrix(i,A,j)),
inverzniA[j,i] : (-1)^(i+j) * detAs / detA,
display(i,j,submatrix(i,A,j),detAs,inverzniA[j,i])
}};
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
Vypsání vypočtené inverzní matice.
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
inverzniA;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
Zjištění inverzní matice pomocí příkazu invert(%).
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
invert(A);
/* [wxMaxima: input   end   ] */

/* Maxima can't load/batch files which end with a comment! */
"Created with wxMaxima"$