xref: /freebsd/contrib/bc/tests/bc/scripts/bessel.bc (revision e17f5b1d)
1#! /usr/bin/bc -q
2
3t[0] = 0
4t[1] = 0.5
5t[2] = 1
6t[3] = 1.5
7t[4] = 1.74
8t[5] = 2
9t[6] = 3.2345
10t[7] = 100
11t[8] = -0.5
12t[9] = -1
13t[10] = -1.5
14t[11] = -1.74
15t[12] = -2
16t[13] = -3.2345
17t[14] = -100
18
19l = 15
20
21a[0] = t[0]
22
23for (i = 1; i < l; ++i) {
24	a[i * 2 - 1] = t[i]
25	a[i * 2] = -t[i]
26}
27
28l *= 2
29l -= 1
30
31for (i = 0; i < l; ++i) {
32	for (j = 0; j < l; ++j) {
33		print "j(", a[i], ", ", a[j], ")\n"
34	}
35}
36
37# These are specific tests that bc could not pass at one time.
38print "j(3, 0.75)\n"
39print "scale = 0; j(40, 0.75)\n"
40
41print "halt\n"
42
43halt
44