archive/share/trash/numer.dmo

   1
   2 /* a demo of defining floating point functions,
   3    and floating point complex functions,
   4    macsyma macros, and the macsyma compiler
   5    for more details, see the numer directory.*/
   6 /* 5:43pm  thursday, 20 march 1980
   7    if this demo screws up it is because something has changed and
   8    it is out of date. please tell gjc. */
   9
  10 /* showtime all should be done to keep track of the efficiency of
  11    the code. */
  12
  13 (showtime:'all,loadfile(numer,autoload,dsk,numer))$
  14
  15 /* make an expression */
  16
  17 exp1:taylor(exp(-x^2),x,0,5);
  18
  19 /* make a numerical function. floatdefunk takes a math expression. */
  20
  21 floatdefunk('f,['x],exp1);
  22
  23 /* lets see what translation and compilation can do for us
  24    in this case. use a function testfn which will call romberg
  25    of f, n times. */
  26
  27 testfn(n):=(modedeclare(n,fixnum),for j:2 thru n do romberg(f,