On Sat, 2011-05-21 at 22:27 +0300, Jussi Lahtinen wrote:
> Gambas!
> http://www.wolframalpha.com/input/?i=besselj0%285%29
>
> Jussi
>
Here is the result of the C version.
gcc -Wall intro.c -lgslcblas -lgsl -o intro
J0(5) = -1.775967713143382642471124199801124632358551025391e-01
For grins I DEFINED the printf() and used it instead of the SB FORMAT()
function.
DEFINE "printf NONE DOUBLE 3 STRING DOUBLE DOUBLE"
DLL("printf \"J0(%g) = %.48e\n\" 5.0 " & DLL("gsl_sf_bessel_J0 " & 5.0))
J0(5) = -1.775970000000000048601123125990852713584899902344e-01
<further testing ...>
v = 1.775967713143382642471124199801124632358551025391e-01
PRINT FORMAT("J0(%g) = %.48e\n", 5.0, v)
J0(5) = 1.775967713143382642471124199801124632358551025391e-01
I'm convinced I have a GTK-Server issue with DOUBLE returns.
------------------------------------------------------------------------------
What Every C/C++ and Fortran developer Should Know!
Read this article and learn how Intel has extended the reach of its
next-generation tools to help Windows* and Linux* C/C++ and Fortran
developers boost performance applications - including clusters.
http://p.sf.net/sfu/intel-dev2devmay
_______________________________________________
Gambas-user mailing list
Gambas-user@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/gambas-user
