diff --git a/Scattering and Shit.lyx b/Scattering and Shit.lyx index 14ee495..1cdefc7 100644 --- a/Scattering and Shit.lyx +++ b/Scattering and Shit.lyx @@ -532,7 +532,7 @@ Definition [T](2.40); \begin_inset Formula \begin{eqnarray*} \widetilde{\vect N}_{mn}^{(j)} & = & \frac{n(n+1)}{kr}\sqrt{\frac{2n+1}{4\pi}\frac{\left(n-m\right)!}{\left(n+m\right)!}}P_{n}^{m}\left(\cos\theta\right)e^{im\phi}z_{n}^{j}\left(kr\right)\hat{\vect r}\\ - & & +\left[\tilde{\tau}_{mn}\left(\cos\theta\right)\hat{\vect{\theta}}+i\tilde{\pi}_{mn}\left(\cos\theta\right)\hat{\vect{\phi}}\right]e^{im\phi}z_{n}^{j}\left(kr\right)\\ + & & +\left[\tilde{\tau}_{mn}\left(\cos\theta\right)\hat{\vect{\theta}}+i\tilde{\pi}_{mn}\left(\cos\theta\right)\hat{\vect{\phi}}\right]e^{im\phi}\frac{1}{kr}\frac{\ud\left(kr\, z_{n}^{j}\left(kr\right)\right)}{\ud(kr)}\\ \widetilde{\vect M}_{mn}^{(j)} & = & \left[i\tilde{\pi}_{mn}\left(\cos\theta\right)\hat{\vect{\theta}}-\tilde{\tau}_{mn}\left(\cos\theta\right)\hat{\vect{\phi}}\right]e^{im\phi}z_{n}^{j}\left(kr\right) \end{eqnarray*} diff --git a/qpms/qpms_p.py b/qpms/qpms_p.py index e5e511b..a7b8bc3 100644 --- a/qpms/qpms_p.py +++ b/qpms/qpms_p.py @@ -423,7 +423,7 @@ def Ã(m,n,μ,ν,kdlj,θlj,φlj,r_ge_d,J): +math.lgamma(n+ν+1) - math.lgamma(2*(n+ν)+1)) presum = math.exp(exponent) presum = presum * np.exp(1j*(μ-m)*φlj) * (-1)**m * 1j**(ν+n) / (4*n) - qmax = floor(q_max(-m,n,μ,ν)) #nemá tu být +m? + qmax = math.floor(q_max(-m,n,μ,ν)) #nemá tu být +m? q = np.arange(qmax+1, dtype=int) # N.B. -m !!!!!! a1q = a_q(-m,n,μ,ν) # there is redundant calc. of qmax @@ -455,7 +455,7 @@ def B̃(m,n,μ,ν,kdlj,θlj,φlj,r_ge_d,J): presum = math.exp(exponent) presum = presum * np.exp(1j*(μ-m)*φlj) * (-1)**m * 1j**(ν+n+1) / ( (4*n)*(n+1)*(n+m+1)) - Qmax = floor(q_max(-m,n+1,μ,ν)) + Qmax = math.floor(q_max(-m,n+1,μ,ν)) q = np.arange(Qmax+1, dtype=int) if (μ == ν): # it would disappear in the sum because of the factor (ν-μ) anyway ã2q = 0