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