Ilya Zakharevich on Sat, 10 May 2003 15:29:23 -0700


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

CVS: seriesprecision too large


? default(seriesprecision,9)
   seriesprecision = 9 significant terms
? sin(x+y)
      /     1 3    1  5    1   7     1    9     10  \
%2 =  | y - -y  + ---y  - ----y  + ------y + O(y  ) | +
      \     6     120     5040     362880           /
/     1 2   1  4    1  6     1   8      1    10     11  \
| 1 - -y  + --y  - ---y  + -----y  - -------y  + O(y  ) | x +
\     2     24     720     40320     3628800            /
/ -1    1  3    1  5     1   7     1    9     10  \  2
| --y + --y  - ---y  + -----y  - ------y + O(y  ) | x  +
\ 2     12     240     10080     725760           /
/ -1   1  2    1  4    1   6     1    8      1     10     11  \  3
| -- + --y  - ---y  + ----y  - ------y  + --------y  + O(y  ) | x  +
\ 6    12     144     4320     241920     21772800            /
/ 1      1  3    1   5     1    7      1    9     10  \  4
| --y - ---y  + ----y  - ------y  + -------y + O(y  ) | x  +
\ 24    144     2880     120960     8709120           /
/  1     1  2    1   4     1   6      1    8       1     10     11  \  5
| --- - ---y  + ----y  - -----y  + -------y  - ---------y  + O(y  ) | x  +
\ 120   240     2880     86400     4838400     435456000            /
/ -1      1   3     1   5      1    7       1     9     10  \  6   /  -1
| ---y + ----y  - -----y  + -------y  - ---------y + O(y  ) | x  + | ---- +
\ 720    4320     86400     3628800     261273600           /      \ 5040
---8<---8<---
   1   2     1    4      1    6       1     8        1      10     11  \  7
 -----y  - ------y  + -------y  - ---------y  + -----------y  + O(y  ) | x
 10080     120960     3628800     203212800     18289152000            /
  /   1        1    3      1    5       1     7        1      9     10  \
+ | -----y - ------y  + -------y  - ---------y  + -----------y + O(y  ) |
  \ 40320    241920     4838400     203212800     14631321600           /
 8      9
x  + O(x )

See the terms with y^10.  It is fine with 2.2.0.

Yours,
Ilya