/* 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")
*/

/*
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:
*/