X-Git-Url: https://ruderich.org/simon/gitweb/?a=blobdiff_plain;f=maxima%2Fmaxima-init.mac;h=1dedac1ba79932fc63ab1222536bf25db42be6a6;hb=576296b2eac6ae7620c8bacdce64b6841df3f05c;hp=f270c10f3f5bef5153c78a901547447d4fc16c57;hpb=f4b44743d9546dbe42c109b2a25f25a1a51d8bbb;p=config%2Fdotfiles.git diff --git a/maxima/maxima-init.mac b/maxima/maxima-init.mac index f270c10..1dedac1 100644 --- a/maxima/maxima-init.mac +++ b/maxima/maxima-init.mac @@ -1,7 +1,32 @@ /* - * Maxima (computer algebra system) configuration file, maxima part. + * Maxima (computer algebra system) configuration file, Maxima part. * - * Also have a look at maxima/maxima-init.mac. + * Also have a look at maxima/maxima-init.lisp. + * + * Notes: + * + * Solve Ax = b, thanks to Florian. + * + * A : matrix([1,2,0],[0,2,3],[0,0,4])$ + * b : matrix([1],[3],[5])$ + * linsolve_by_lu(A,b); + */ + +/* + * Copyright (C) 2011-2012 Simon Ruderich + * + * This file is free software: you can redistribute it and/or modify + * it under the terms of the GNU General Public License as published by + * the Free Software Foundation, either version 3 of the License, or + * (at your option) any later version. + * + * This file is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the + * GNU General Public License for more details. + * + * You should have received a copy of the GNU General Public License + * along with this file. If not, see . */ @@ -13,7 +38,7 @@ norm(x) := sqrt(transpose(x) . x); /* Maxima's mat_norm() can't calculate the 2-norm. Thanks to Wolfgang Lindner * (http://www.ma.utexas.edu/pipermail/maxima/2007/006300.html) for an * improved version. */ -mat_norm2(x) := sqrt(lmax(eigenvalues(transpose(conjugate(x)).x)[1])); +mat_norm2(x) := sqrt(lmax(eigenvalues(transpose(conjugate(x)) . x)[1])); /* Shortcut to perform lu factorization. */ lu(x) := block(x : lu_factor(x), get_lu_factors(x));