[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
## [Axiom-math] indexed operators

**From**: |
Martin Rubey |

**Subject**: |
[Axiom-math] indexed operators |

**Date**: |
Thu, 9 Dec 2004 18:15:23 +0100 |

I worked out how to program the indexed operators stuff, see below. Also, I'd
like to supplement something I called a bug before:
Given
qs:=[operator(subscript('q, [i::OutputForm])::Symbol) for i in 1..3]
I said, that
(qs.1(2*t)
does not work.
(qs.1)(2*t::EXPR INT)
does, however.
Here is the promised code for indexed operators and setting their derivatives:
)abbrev package TEST Test
Test(): Exports == Implementation where
EI == Expression Integer
Exports == with
tst: () -> EI
Implementation == add
tst () ==
qs:=[operator(subscript('q, [i::OutputForm])::Symbol) for i in 1..3]
ps:=[operator(subscript('p, [i::OutputForm])::Symbol) for i in 1..3]
for i in 1..3 repeat
derivative((qs.i)::BasicOperator,(ps.i)(#1))$BasicOperatorFunctions1(EI)
var:=t::Symbol::EI
D((qs.1)(2*var)+(qs.2)(var),t::Symbol)
Martin

[Prev in Thread] |
**Current Thread** |
[Next in Thread] |

**[Axiom-math] indexed operators**,
*Martin Rubey* **<=**