Olivier Ramare on Fri, 26 Feb 2010 15:09:22 +0100

[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

 Strange behaviour of intnum A12-140

```Hello,

(1)----------------------------------------------------------------------

Can anyone tell me what wrong with

intnum(t = [0,0], [1], (sin(t)/t)^2)

Gp answers

? intnum(t = [0,0], [1], (sin(t)/t)^2)
***   at top-level: intnum(t=[0,0],[1],(sin(t)/t)^2)
***                                     ^------------
*** sin: precision too low in mpsc1.

intnum(t = [0,0], 10000, (sin(t)/t)^2)
is ok, though the precision is of course fairly awful.

(2) ----------------------------------------------------------------------

intnum(t = [0,0], 1000, 9*((sin(t)-t*cos(t))/t^3)^2)

is -- 0 !!!
-------------------------------------------------------------------------
? intnum(t = [0,0], 1000, 9*((sin(t)-t*cos(t))/t^3)^2)
%19 = 0.E94
? intnum(t = 0.00001, 1000, 9*((sin(t)-t*cos(t))/t^3)^2)
%20 = 1.884938921515439810798085699
---------------------------------------------------------------------------
intnum(t = 0.00001, [1], 9*((sin(t)-t*cos(t))/t^3)^2)
***   at top-level: intnum(t=0.00001,[1],9*((sin(t)-t*cos(t))/t^3
***                                          ^--------------------
*** sin: precision too low in mpsc1.
---------------------------------------------------------------------------

Many thanks!
Amities,
Olivier

```