|Bill Allombert on Tue, 18 Jun 2002 16:26:24 +0200|
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
On Mon, Jun 17, 2002 at 07:31:31PM +0200, Karim BELABAS wrote: > Received the following bug report: > > (19:10) gp > intnum(x=0,Pi,sin(16*x)^2) > %1 = 1.855125350366532182 E-54 \\ should be Pi/2 > > Cause: integral computed using Romberg method, with equally spaced points > (+ interpolation of successive values for a few runs with increasing number > of points). For highly oscillating functions (esp. integrating over 2^k * > period), we may be unlucky. > > Possible solution: pick random point in each subinterval ? ? intnum(x=0,Pi,sin(16*x)^2,3) %1 = 1.570796326794896619231321691 So what trick does this flag use ? Cheer, Bill.