Groups | Search | Server Info | Keyboard shortcuts | Login | Register [http] [https] [nntp] [nntps]

Groups > comp.lang.forth > #21669

RE: ANN: Port of Airy Function Code to Forth

Show all headers | 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