To comment on the following update, log in, then open the issue:
http://www.openoffice.org/issues/show_bug.cgi?id=19991





------- Additional comments from [EMAIL PROTECTED] Thu Oct  6 19:08:02 -0700 
2005 -------
Compile and run the test program as follows:

[EMAIL PROTECTED] Erf]$ g++ -Os -o erf erf.cxx
[EMAIL PROTECTED] Erf]$ ./erf 5
<maclaurin>
erf:  1.000000000072773
erfc: -7.2773443412881e-11
<asymptotic>
erf:  0.9999999999984626
erfc: 1.537436844500917e-12

The maclaurin series expansion gives Calc's current answers for erf and erfc,
while the asymptotic series expansion (proposed new algorithm) give more
accurate answers.

You can also run this program with no arguments.  It will compute erf and erfc
using both the Maclaurin expansion (current algorithm) and the Asymptotic
expansion (new algorithm) for x between -10 and 10 incremented by 0.1.

[EMAIL PROTECTED] Erf]$ ./erf
x       erf(x) <maclaurin>      erfc(x) <maclaurin>     erf(x) <asymptotic>    
erfc(x) <asymptotic>
-10     7.331381385975313e+21   -7.331381385975313e+21  -1      2
-9.9    -3.548111434363548e+19  3.548111434363548e+19   -1      2
-9.8    -1.744418020013381e+20  1.744418020013381e+20   -1      2
-9.7    -1.560844956904942e+19  1.560844956904942e+19   -1      2
<snip>
9.7     -7.659310690090643e+18  7.659310690090643e+18   1       0
9.8     3.709705234685252e+19   -3.709705234685252e+19  1       0
9.9     4.120589789696997e+20   -4.120589789696997e+20  1       0
10      9.970406838904948e+21   -9.970406838904948e+21  1       0



---------------------------------------------------------------------
Please do not reply to this automatically generated notification from
Issue Tracker. Please log onto the website and enter your comments.
http://qa.openoffice.org/issue_handling/project_issues.html#notification

---------------------------------------------------------------------
To unsubscribe, e-mail: [EMAIL PROTECTED]
For additional commands, e-mail: [EMAIL PROTECTED]


---------------------------------------------------------------------
To unsubscribe, e-mail: [EMAIL PROTECTED]
For additional commands, e-mail: [EMAIL PROTECTED]

Reply via email to