f2c04ccdad 2021-03-03 trnsz@pobox.c: %%%%%%%%%%%%%%%%%%%%%
f2c04ccdad 2021-03-03 trnsz@pobox.c: % LINEAR ALGEBRA
f2c04ccdad 2021-03-03 trnsz@pobox.c: %%%%%%%%%%%%%%%%%%%%%
f2c04ccdad 2021-03-03 trnsz@pobox.c:
f2c04ccdad 2021-03-03 trnsz@pobox.c: % Enable fancy output for easier viewing of
f2c04ccdad 2021-03-03 trnsz@pobox.c: % matrix output
f2c04ccdad 2021-03-03 trnsz@pobox.c: fancy_output$
f2c04ccdad 2021-03-03 trnsz@pobox.c:
f2c04ccdad 2021-03-03 trnsz@pobox.c: % Load linear algebra package
f2c04ccdad 2021-03-03 trnsz@pobox.c: load_package linalg$
f2c04ccdad 2021-03-03 trnsz@pobox.c:
f2c04ccdad 2021-03-03 trnsz@pobox.c: % Define a complex 3x3 matrix
f2c04ccdad 2021-03-03 trnsz@pobox.c: m1 := mat((1+i*3, 2-i*5, 7-i), (4-i*2, 6+i*9,-8+i*4), (-3-i*7, 3+i*2, -1+i*6));
f2c04ccdad 2021-03-03 trnsz@pobox.c:
f2c04ccdad 2021-03-03 trnsz@pobox.c: % Determinant of matrix
f2c04ccdad 2021-03-03 trnsz@pobox.c: write "|m1| = ", det(m1)$
f2c04ccdad 2021-03-03 trnsz@pobox.c:
f2c04ccdad 2021-03-03 trnsz@pobox.c: % Trace of matrix
f2c04ccdad 2021-03-03 trnsz@pobox.c: write "trace(m1) = ", trace(m1)$
f2c04ccdad 2021-03-03 trnsz@pobox.c:
f2c04ccdad 2021-03-03 trnsz@pobox.c: % Characteristic polynomial
f2c04ccdad 2021-03-03 trnsz@pobox.c: write "characteristic polynomial of m1:";
f2c04ccdad 2021-03-03 trnsz@pobox.c: char_poly(m1,eta);
f2c04ccdad 2021-03-03 trnsz@pobox.c:
f2c04ccdad 2021-03-03 trnsz@pobox.c: % Enable real arithmetic
f2c04ccdad 2021-03-03 trnsz@pobox.c: on rounded$
f2c04ccdad 2021-03-03 trnsz@pobox.c:
f2c04ccdad 2021-03-03 trnsz@pobox.c: % Singular value decomposition of a matrix.
f2c04ccdad 2021-03-03 trnsz@pobox.c: a := mat((1,3),(-4,3));
f2c04ccdad 2021-03-03 trnsz@pobox.c: write "Singular Value Decomposition of a:"$
f2c04ccdad 2021-03-03 trnsz@pobox.c: svd(a);