--- /dev/null
+/* Test of jn() function.
+ Copyright (C) 2010 Free Software Foundation, Inc.
+
+ This program is free software: you can redistribute it and/or modify
+ it under the terms of the GNU General Public License as published by
+ the Free Software Foundation; either version 3 of the License, or
+ (at your option) any later version.
+
+ This program is distributed in the hope that it will be useful,
+ but WITHOUT ANY WARRANTY; without even the implied warranty of
+ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+ GNU General Public License for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with this program. If not, see <http://www.gnu.org/licenses/>. */
+
+/* Written by Bruno Haible <bruno@clisp.org>, 2010. */
+
+#include <config.h>
+
+#include <math.h>
+
+#include "signature.h"
+SIGNATURE_CHECK (jn, double, (int, double));
+
+#include "macros.h"
+
+volatile double x;
+double y;
+
+int
+main ()
+{
+ /* A particular value. */
+ x = 3.8;
+ y = jn (0, x);
+ ASSERT (y >= -0.4025564102 && y <= -0.4025564101);
+
+ /* A particular value. */
+ x = 3.8;
+ y = jn (1, x);
+ ASSERT (y >= 0.0128210029 && y <= 0.0128210030);
+
+ /* A particular value. */
+ x = 3.8;
+ y = jn (2, x);
+ ASSERT (y >= 0.4093043064 && y <= 0.4093043065);
+
+ /* A particular value. */
+ x = 3.8;
+ y = jn (3, x);
+ ASSERT (y >= 0.4180256354 && y <= 0.4180256355);
+
+ /* A particular value. */
+ x = 3.8;
+ y = jn (4, x);
+ ASSERT (y >= 0.2507361705 && y <= 0.2507361706);
+
+ /* A particular value. */
+ x = 3.8;
+ y = jn (5, x);
+ ASSERT (y >= 0.1098399867 && y <= 0.1098399868);
+
+ /* A particular value. */
+ x = 3.8;
+ y = jn (6, x);
+ ASSERT (y >= 0.0383164262 && y <= 0.0383164263);
+
+ /* A particular value. */
+ x = 3.8;
+ y = jn (7, x);
+ ASSERT (y >= 0.0111592540 && y <= 0.0111592541);
+
+ /* A particular value. */
+ x = 3.8;
+ y = jn (8, x);
+ ASSERT (y >= 0.0027966149 && y <= 0.0027966150);
+
+ /* A particular value. */
+ x = 3.8;
+ y = jn (9, x);
+ ASSERT (y >= 0.0006159669 && y <= 0.0006159670);
+
+ /* A particular value. */
+ x = 3.8;
+ y = jn (10, x);
+ ASSERT (y >= 0.0001211233 && y <= 0.0001211234);
+
+ return 0;
+}