RE: ANN: Port of Airy Function Code to Forth

Path	csiph.com!v102.xanadu-bbs.net!xanadu-bbs.net!feeder.erje.net!eu.feeder.erje.net!xlned.com!feeder5.xlned.com!news2.euro.net!postnews2.euro.net!news.wanadoo.nl!not-for-mail
From	mhx@iae.nl (Marcel Hendrix)
Subject	RE: ANN: Port of Airy Function Code to Forth
Newsgroups	comp.lang.forth
Message-ID	<68571488998434@frunobulax.edu> (permalink)
Date	Mon, 15 Apr 2013 21:46:35 +0200
References	<994580ac-8143-4114-b96c-427397ca7eb4@googlegroups.com>
X-Newsreader	iForth 2.0 console (October 21, 2006)
Lines	27
Organization	Wanadoo
NNTP-Posting-Date	15 Apr 2013 19:46:31 GMT
NNTP-Posting-Host	s529d937f.adsl.online.nl
X-Trace	1366055191 dr4.euro.net 9205 82.157.147.127:60974
X-Complaints-To	abuse@wanadoo.nl
Xref	csiph.com comp.lang.forth:21669

Show key headers only | View raw

krishna.myneni@ccreweb.org writes RE: ANN: Port of Airy Function Code to Forth

> A first version of my port of the Fortran routine, JAIRY, for computing the
> Airy function, and its derivative, for a real argument, is available at

> ftp://ccreweb.org/software/fsl/extras/jairy.fs
[..]
> TESTING JAIRY
> ok

I get one error, for -10.0e jairy .
According to Wolfram Alpha, the result should be:

  t{ -10.0e jairy -> 4.0241238486443190689430e-2  9.962650441327900e-1   ?}t
\                                  ^^^^^^^^^^
However, the current file says
  t{ -10.0e jairy -> 4.024123848644296e-2         9.962650441327906e-1  ?}t
\                                  ^^^

iForth computes:
FORTH> -10e jairy cr fswap +e. space +e.
 4.0241238486443971872e-0002  9.9626504413279069178e-0001 ok

This tests OK with the Alpha values, but not with the original test result.

-marcel

Back to comp.lang.forth | Previous | Next — Previous in thread | Next in thread | Find similar

Thread

ANN: Port of Airy Function Code to Forth krishna.myneni@ccreweb.org - 2013-04-14 19:36 -0700
  RE: ANN: Port of Airy Function Code to Forth mhx@iae.nl (Marcel Hendrix) - 2013-04-15 21:46 +0200
    Re: ANN: Port of Airy Function Code to Forth krishna.myneni@ccreweb.org - 2013-04-15 16:32 -0700
      Re: ANN: Port of Airy Function Code to Forth krishna.myneni@ccreweb.org - 2013-04-15 17:54 -0700

csiph-web