Fix summation limit in pipes._F_mk()

MassimoCimmino · Feb 12, 2018 · 1f59872 · 1f59872
1 parent 8f974e8
commit 1f59872
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/pygfunction/pipes.py b/pygfunction/pipes.py
@@ -2126,7 +2126,7 @@ def _F_mk(Q_p, P, n_p, J, r_b, r_p, z, pikg, sigma):
                                 *r_p[n]**(j+1)*(-r_p[m])**(k+1) \
                                 /(z[m] - z[n])**(j+k+2)
                     # Fourth term
-                    j_pend = np.min((k, j)) + 1
+                    j_pend = np.min((k, j)) + 2
                     for jp in range(j_pend):
                         fmk += sigma*np.conj(P[n,j])*binom(j+1, jp) \
                                 *binom(j+k-jp+1, j)*r_p[n]**(j+1) \