[sage-support] Re: Approximating integral with infinite bounds

2020-03-01 Thread Emmanuel Charpentier
Le samedi 29 février 2020 23:06:36 UTC+1, Simon King a écrit : > > Hi Emmanuel, > > On 2020-02-29, Emmanuel Charpentier > > wrote: > > This question would have been more properly posed on ask.sagemath.org... > > Why? > As I understand our current documentation : sage-support is oriented tow

[sage-support] Re: Approximating integral with infinite bounds

2020-02-29 Thread Simon King
Hi Emmanuel, On 2020-02-29, Emmanuel Charpentier wrote: > This question would have been more properly posed on ask.sagemath.org... Why? I for one totally dislike ask.sagemath.org and would never post a question or an answer there. It is of course a matter of taste. But it is certainly not appro

Re: [sage-support] Re: Approximating integral with infinite bounds

2020-02-29 Thread Vincent Delecroix
Le 29/02/2020 à 21:41, Emmanuel Charpentier a écrit : You may try to separate the real and imaginary part of your function, and sample from these part. Look at Stan and the companion R package bridgesampling

[sage-support] Re: Approximating integral with infinite bounds

2020-02-29 Thread Emmanuel Charpentier
This question would have been more properly posed on ask.sagemath.org... The firs argument of numerical integral must be a function of one (real) variable returning a real. It is indeed possible to nest such integrals, but the second one *nests* the first one, hence long computation times. Exam

Re: [sage-support] Re: Approximating integral with infinite bounds

2020-02-27 Thread Montgomery-Smith, Stephen
Split it into it's real and imaginary parts maybe? From: sage-support@googlegroups.com on behalf of saad khalid Sent: Thursday, February 27, 2020 2:51:07 PM To: sage-support Subject: Re: [sage-support] Re: Approximating integral with infinite bounds Hi,

Re: [sage-support] Re: Approximating integral with infinite bounds

2020-02-27 Thread saad khalid
o:* sage-support > *Subject:* [sage-support] Re: Approximating integral with infinite bounds > > Also, I tried running the integral with finite bounds, and I get a giac > error. Here is the code: > var('t1,t2,u,w,k') > T = 1 > m = 100 > E = 1 > v = 0 > y=1 > O = 1

Re: [sage-support] Re: Approximating integral with infinite bounds

2020-02-27 Thread Montgomery-Smith, Stephen
, 2020 9:41:19 PM To: sage-support Subject: [sage-support] Re: Approximating integral with infinite bounds Also, I tried running the integral with finite bounds, and I get a giac error. Here is the code: var('t1,t2,u,w,k') T = 1 m = 100 E = 1 v = 0 y=1 O = 1 integral(integral(integral(

[sage-support] Re: Approximating integral with infinite bounds

2020-02-26 Thread saad khalid
Also, I tried running the integral with finite bounds, and I get a giac error. Here is the code: var('t1,t2,u,w,k') T = 1 m = 100 E = 1 v = 0 y=1 O = 1 integral(integral(integral( integral(integral( e^(-t1^2/T^2)*e^(-t2^2/T^2)*e^(I*O*t1)* e^(-I*O*t2)*e^(-I*E*y^2*(1 - v)*t1^2/2)*