Fortran output



>    My problem involves some mundane but huge expansions of products of
>    polynomials. I don't care whether the Fortran is human-legible or
>    not. Is it worth optimizing the expressions with optimize()
(turning
>    the block statements into a Fortran sequence by hand) or should I
not
>    bother, and just let the compiler do the optimization?

A good Fortran compiler at full optimization level should do better than
Optimize (which is really very simple-minded).

>    Is there a way to get "optimize" to handle multiple expressions,
>    i.e. am computing a sequence of array elements, and there are
>    apparently lots of common subexpressions among the individual array
>    elements.

optimize([expr1,expr2,expr3]) or even optimize(matrix([a,b],[c,d]).

> 4. On the subject of optimization, is there a way to get the
> Newton form of a polynomial?

I don't think there is one built-in, or in the Share libraries.  If you
write it as Maxima code, I'm sure it would be a good addition to the
Share libraries.

       -s