/* CUSTOM FUNCTIONS */
+/* Didn't find any euclid norm in Maxima, so here it is. */
+norm(x) := sqrt(transpose(x) . x);
+
/* Maxima's norm() can't calculate the 2-norm. Only works for non-complex
* matrices. */
norm_2(x) := sqrt(maxi(eigenvalues(transpose(x) . x)[1]));