Maybe the following will help.
Also note that I've replaced GENSYM()$Lisp by new(). Look up that
function in HyperDoc.
Ralf
mysum := operator 'mysum
mi -- represents minusInfinity
pi -- represents plusInfinity
B(N, l, h, a, t) == mysum(l, h, N, besselJ(N, a)*exp(sqrt(-1)*N*t))
bess(a,t) == (N:Symbol:=new('p);B(N,mi,pi,a,t))
test_eqn := exp(sqrt(-1)*a*sin(omg*t))*exp(sqrt(-1)*b*sin(2*omg*t))
exp_to_bessel_rule := rule exp(sqrt(-1)*z*sin(t)) == bess(z,t)
r := exp_to_bessel_rule( test_eqn )
ks := [first kernels x for x in isTimes r]
summationVars := [argument(k).3 for k in ks]
eval(r,[summationVars.1=v1,summationVars.2=v2])
--
You received this message because you are subscribed to the Google Groups
"FriCAS - computer algebra system" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at http://groups.google.com/group/fricas-devel.
For more options, visit https://groups.google.com/d/optout.
