>
> In case I forgot the attachment...
>
>
--
You received this message because you are subscribed to the Google Groups
"sympy" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To view this discussion on the web visit
https://groups.google.com/d/msgid/sympy/31287eaf-9add-4500-b83a-6bf8749573f4%40googlegroups.com.
from sympy import symbols, Function, latex, Eq, Derivative, diff, sin, cos, Tuple, solve
a0, a1, a2, a3, a4, a5 = symbols('a_0 a_1 a_2 a_3 a_4 a_5')
am1, am2, am3, am4, am5 = symbols('a_{-1} a_{-2} a_{-3} a_{-4} a_{-5}')
E10 = symbols('E_{1\,0}')
E11 = symbols('E_{1\,1}')
E1m1 = symbols('E_{1\,-1}')
F1 = symbols('F_1')
F1 = symbols('F_1')
def tripLet(j, k):
strE = ('E'+str(k)+str(j)).replace('-', 'm')
stra1 = ('a'+str(j)).replace('-', 'm')
stra2 = ('a'+str(j-k)).replace('-', 'm')
strF = ('F'+str(k)).replace('-', 'm')
stra3 = ('a'+str(-j+k-1)).replace('-', 'm')
strG = ('G'+str(k)).replace('-', 'm')
stra4 = ('a'+str(-j-k-1)).replace('-', 'm')
return symbols(strE)*symbols(stra1)*symbols(stra2)+symbols(strF)*symbols(stra1)*symbols(stra3)+symbols(strG)*symbols(stra1)*symbols(stra4)
# For k = 1
#eq79 = Eq(0, Sum(Ekj*aj*ajmk + Fk*aj*amjpkm1 + Gk*aj*amjmkm1, (j,-1,1)))
#eq79 = eq79.subs(k, 1)
#print('eq79:')
#print(eq79)
#print(latex(eq79))
eq79 = Eq(0, tripLet(-1, 1)+tripLet(0, 1)+tripLet(1, 1))
print('eq79:')
print(eq79)
print(latex(eq79))
print('end')
