maxima/maxima-init.mac: Add. Define norm_2().

diff --git a/maxima/maxima-init.lisp b/maxima/maxima-init.lisp
index ec1dcf06cdd9c6da54d064bdd6532592670e7a2d..82b90d4496872321f6d4b97fa21b9fb0064d88ec 100644 (file)
--- a/maxima/maxima-init.lisp
+++ b/maxima/maxima-init.lisp
@@ -1,4 +1,6 @@
 ; Maxima (computer algebra system) configuration file, lisp part.
+;
+; Also have a look at maxima/maxima-init.mac.
 
 
 ; Set directory for temporary files to ~/.maxima (used for plotting for

diff --git a/maxima/maxima-init.mac b/maxima/maxima-init.mac
new file mode 100644 (file)
index 0000000..438989f
--- /dev/null
+++ b/maxima/maxima-init.mac
@@ -0,0 +1,20 @@
+/*
+ * Maxima (computer algebra system) configuration file, maxima part.
+ *
+ * Also have a look at maxima/maxima-init.mac.
+ */
+
+
+/* PACKAGES */
+
+/* Provides maxi() and mini() (and more). */
+load (descriptive);
+
+
+/* CUSTOM FUNCTIONS */
+
+/* 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]));
+
+/* vim: ft=maxima */