Warning: this is an htmlized version!
The original is here, and the conversion rules are here. |
/* This file: * http://angg.twu.net/MAXIMA/matrixify.mac.html * http://angg.twu.net/MAXIMA/matrixify.mac * (find-angg "MAXIMA/matrixify.mac") * Author: Eduardo Ochs <eduardoochs@gmail.com> * * (defun m () (interactive) (find-angg "MAXIMA/matrixify.mac")) * See: (find-es "maxima" "matrixify") * Superseded by: * (defun m2 () (interactive) (find-angg "MAXIMA/matrixify2.mac")) */ /* matrixifyf : lambda([xmin,ymin, xmax,ymax, f], apply(matrix, makelist(makelist(f(x,y), x, xmin, xmax), y, ymax, ymin, -1))); matrixify : lambda([xmin,ymin, xmax,ymax, expr], apply(matrix, makelist(makelist(ev(expr), x, xmin, xmax), y, ymax, ymin, -1))); */ matrixify(xmin,ymin, xmax,ymax, expr) := apply(matrix, makelist(makelist(ev(expr), x, xmin, xmax), y, ymax, ymin, -1)); /* * (eepitch-maxima) * (eepitch-kill) * (eepitch-maxima) ** (find-angg "MAXIMA/matrixify.mac") load ("~/MAXIMA/matrixify.mac"); [Dx,Dy] : [x-x0,y-y0]; [x0,y0] : [5,4]; matrixify (x0-3,y0-3, x0+3,y0+3, [x,y]); matrixify (x0-3,y0-3, x0+3,y0+3, [Dx,Dy]); matrixify (x0-3,y0-3, x0+3,y0+3, Dx-Dy); matrixify (x0-3,y0-3, x0+3,y0+3, Dy); matrixify (x0-3,y0-3, x0+3,y0+3, Dy*(Dx-Dy)); z : Dy*(Dx-Dy); matrixify (x0-3,y0-3, x0+3,y0+3, z); [x0,y0] : [0,0]; [d,e,f] : [1,1,2]; [d,e,f] : [2,-1,1]; matrixify (-2,-2, 2,2, d*x^2 + e*x*y + f*y^2); [d,e,f] : [1,-1,-2]; matrixify (-2,-2, 2,2, d*x^2 + e*x*y + f*y^2); [d,e,f] : [2,-1,-1]; matrixify (-2,-2, 2,2, d*x^2 + e*x*y + f*y^2); F : d*x^2 + e*x*y + f*y^2; Fx : diff(F,x); Fy : diff(F,y); matrixify (-2,-2, 2,2, [Fx, Fy]); matrixify (-3,-3, 3,3, x*y - y^2); */ /* * Local Variables: * mode: c * modes: (c-mode maxima-mode) * End: */