[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Axiom-developer] series for sums

From: Martin Rubey
Subject: Re: [Axiom-developer] series for sums
Date: Fri, 7 Oct 2005 17:57:22 +0200

Dear all,

I'm just trying to do the same for products, however, I came across a little
mathematical difficulty. It seems to obtain the coefficient of x^k I need to 

* generate all compositions of k

* for each composition (i1,i2,..,il), generate a multiple sum of the form


Of course this is doable (although I will have to use ugly newsymbols for the
j's), but I wonder whether somebody can think of a better alternative.


Martin Rubey writes:
 > Well, here is a stupid alternative, which works better on the surface:
 >       if FE has CombinatorialOpsCategory then
 >         is?(ker,"%defsum" :: Symbol) =>
 >           (summand := exprToUPS(first args,posCheck?,atanFlag)) case 
 > %problem => 
 >             summand
 > --          op := operator("%defsum"::Symbol)$CommonOperators
 > --          mp:FE->FE:=op cons(#1, rest args)
 >           eq := equation(retract(third args)@Symbol, 
 > args.4..args.5)$SegmentBinding(FE)
 >           mp:FE->FE:=summation(eval(#1, retract(second args)@K, third args), 
 > eq)
 >           [map(mp, summand.%series)$UPS]
 > now Axiom displays beautiful sums, but they won't be evaluated because of a 
 > bug
 > equivalent to issue #212.
 > The summation design is seriously broken, I'd say.
 > Martin

reply via email to

[Prev in Thread] Current Thread [Next in Thread]