qpms/tests/s1_constfacs.sage

   1 def s1_constfacs(m, n):
   2     if (m+n) % 2 == 1:
   3         return []
   4     s1consts = list()
   5     for j in range((n-abs(m))+1):
   6         s1consts.append(
   7             -I**(n+1)/2 * (-1)**((n+m)/2) * sqrt((2*n+1)*factorial(n-m)*factorial(n+m))
   8             * (-1)**j / 2**(n-2*j)
   9             / (factorial(j) * factorial((n-m)/2-j) * factorial((n+m)/2-j)) / 2**(2*j-1)
  10             )
  11     return s1consts
  12
  13 for l in range(1, 11):
  14     for m in range (-l, l+1):
  15         print(l, m, [N(x, prec=66) for x in s1_constfacs(m,l)])
  16