Bill,
Thank you for following up on this.
This is a nice illustration of how to extend pan-Axiom to deal with a
specific case. One of the advantages of algorithmic differentiation (AD)
is that we can rapidly differentiate across a range of operators and
expressions. One could re-write the previou
The problem is just the expansion of (1+x)^1.
(1) -> )set message time on
(1) -> g(x) == (1+x)^1
Type: Void
Time: 0 sec
(2) -> eval(D(g(x),x), x=0.1)
Com
Oops, I forgot to add one line of output
(3) -> distribute(%)::Float
(3) 0.7681727501_0912254214 E 418
Type: Float
Time: 0 sec
'distribute' is the operator that