Skip to main content
20 events
when toggle format what by license comment
Jun 8, 2019 at 23:52 vote accept J Bausch
Jun 5, 2019 at 20:33 comment added J Bausch I'll check and accept asap. Thanks so much!!
Jun 1, 2019 at 3:45 history edited Max Alekseyev CC BY-SA 4.0
SageMath code added
May 31, 2019 at 18:28 history edited Max Alekseyev CC BY-SA 4.0
edited body
May 31, 2019 at 18:15 history edited Max Alekseyev CC BY-SA 4.0
Formulae are corrected.
May 31, 2019 at 18:12 comment added Max Alekseyev @JBausch: Good catch! There was an error in my formulae, now corrected. Essentially, we should use $\frac{b}{a}$ instead of $1$ in $I(1,\cdot)$ and $c$ instead of $\infty$ in $I(\cdot,\infty)$. This also addresses the issue with $\infty$ you pointed out earlier.
May 31, 2019 at 18:09 history edited Max Alekseyev CC BY-SA 4.0
Formulae are corrected.
May 30, 2019 at 9:21 comment added J Bausch @MaxAlekseyev did you have a chance to have a look? I've compared this to a numerical integration and yours matches up when $n\rightarrow\infty$, but e.g. in the case I wrote just above there's some discrepancy. No rush!! I really appreciate your help with this - I've worked through your argument but can't really pinpoint a mistake, which is why I ask.
May 27, 2019 at 22:30 comment added J Bausch Another edge case that doesn't work: take $k=3/2, l=0, n=2, a=10, b=1, c=2$, then $I_2^{3/2,0}(1,2)=-2(0-0-I^{3/2,0}_1(1,\infty)+I^{1,0}_1(1/5,\infty)/\sqrt 2)$; both of the latter integrals can be calculated exactly; the total epression turns out to be $\approx-0.52$, which cannot be since the integrand is strictly positive.
May 27, 2019 at 21:43 comment added J Bausch You're right; treating $\infty$ as a symbolic large constant and defining its product with zero to be zero does the trick. Will check the rest and report back whether it matches a numeric integration!
May 27, 2019 at 19:51 comment added J Bausch @MaxAlekseyev: I've just tried implementing it in MM, but there seems to be something amiss: for e.g. $a=20, b=1, c=40, n=3, l=0$ we have: $I_3^{k,0}(1,40)$ expands the rightmost term $I_{2}^{1,0}(2,\infty)$; the latter then includes a term that has $\log(\infty)I^{1,0}_1(\infty,\infty)$ in it. Am I missing a break condition?
May 24, 2019 at 16:01 comment added J Bausch Very useful answer, much appreciated. I'll try to implement this in MM and cross-check with the original expression.
May 23, 2019 at 15:36 comment added Max Alekseyev @MattF.: My estimate is that $I^{1,l}_n(b,c)$ has $O(n)$ terms and $I^{k,0}_n(b,c)$ has $O(n^2)$ terms.
May 22, 2019 at 21:15 history edited Max Alekseyev CC BY-SA 4.0
added a simpler formula; deleted 6 characters in body
May 22, 2019 at 20:24 comment added user44143 Do you have an estimate for the number of terms in this calculation? If the calculation of $I_n$ uses at least $4$ terms of $I_{n-1}$, and so it uses at least $4^n$ basic terms, then it'd be nicer to get something with polynomial growth.
May 22, 2019 at 18:46 history edited Max Alekseyev CC BY-SA 4.0
added 43 characters in body
May 22, 2019 at 18:41 history edited Max Alekseyev CC BY-SA 4.0
added 43 characters in body
May 22, 2019 at 18:28 history edited Max Alekseyev CC BY-SA 4.0
deleted 6 characters in body
May 22, 2019 at 18:21 history edited Max Alekseyev CC BY-SA 4.0
added 926 characters in body
May 22, 2019 at 16:50 history answered Max Alekseyev CC BY-SA 4.0