[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 <=