Matice
| Reaktanty | Produkty | |||
| ZnO | NaOH | H2O | Na2[Zn(OH)4] | |
| a | b | c | p | |
| Zn | 1 | 1 | ||
| O | 1 | 1 | 1 | 4 |
| Na | 1 | 2 | ||
| H | 1 | 2 | 4 | |
| náboj | ||||
Bilance prvků
|
|
+ 1·a | = | + 1·p |
|
|
+ 1·a + 1·b + 1·c | = | + 4·p |
|
|
+ 1·b | = | + 2·p |
|
|
+ 1·b + 2·c | = | + 4·p |
Bilance elektronů (náboje)
Zadání pro program Mathematica
eqns = {
+ 1*a== + 1*p,
+ 1*a + 1*b + 1*c== + 4*p,
+ 1*b== + 2*p,
+ 1*b + 2*c== + 4*p,
+0*a +0*b +0*c== +0*p}
Solve[eqns]
Neznámých koeficientů je: 4, počet nezávislých rovnic je: NaN. Počet stupňů volnosti je tedy: 4 - NaN = NaN. Jedno z možných řešení je:
Zadání (program Octave/Matlab) reaction_id-10-30.m
% % Jiri Jirat % Prague Institute of Chemical Technology % % % matice - 1. sloupec naboj, dalsi sloupce prvky % a = [ 0,0,0,1,1; 0,1,1,1,0; 0,2,0,1,0; 0,4,2,4,1] hodnost = rank(a) % hodnost matice = pocet nezavislych rovnic b = a' % transpozice matice c = null(b) % nalezeni baze nuloveho prostoru matice b reseni = rref(c') % upravy na "row reduced echelon form"
Řešení (program Octave/Matlab)
Zadání (program Mathematica)
m = {
{0,0,0,1,1},
{0,1,1,1,0},
{0,2,0,1,0},
{0,4,2,4,1}}
NullSpace[Transpose[m]]