[Axiom-math] indexed operators

[Martin Rubey]

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