share/numeric/submac.dem

   1 /* Demo of macsyma substitution macros. */
   2
   3 if get('submac,'version)=false
   4 then load('submac)$
   5
   6 /* define a macro to compute a numerical integration.
   7    We want it to take an expression, a variable, and upper and
   8    lower limits. The expression will be evaluated for values
   9    of the variable between the upper and lower limits. */
  10
  11 rect_rule('exp,'x,a,b,dx)=>
  12  (modedeclare([a,b,dx],float),
  13   block([%_sum:0.0],modedeclare(%_sum,float),
  14         for x:a thru b step dx do %_sum:%_sum+exp,
  15         dx*%_sum))$
  16
  17 grind(rect_rule);
  18
  19 define(g(a),macroexpand(rect_rule(x^2,x,0.0,a,0.1)))$
  20 grind(g);