blob: 4730a2b6adc869f4e055f1e171d092e69eecc774 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
|
! driver for testing function * bessj0*
REAL bm5 bm4 bm3 bm2 bm1
b00 b01 b02 b03 b04
b05 b06 b07 b08 b09
b10 b11 b12 b13 b14
b15 :=
-0.1775968 -0.3971498 -0.2600520 0.2238908 0.7651976
1.0000000 0.7651977 0.2238908 -0.2600520 -0.3971498
-0.1775968 0.1506453 0.3000793 0.1716508 -0.0903336
-0.2459358 -0.1711903 0.0476893 0.2069261 0.1710735
-0.0142245 ;
REAL x := -5.0 ;
REAL y ;
PROCEDURE bessj0 ;
ECHO "Bessel Function J0" ;
WHILE x 16. < DO
IF x 0. <> THEN
bessj0 :: <<x>> >>y<< ;
ECHO "x=" x "bessj0(x)=" y "reference=" bm5 ;
ENDIF ;
EVALUATE x := x 1. + ;
EVALUATE bm5 bm4 bm3 bm2 bm1
b00 b01 b02 b03 b04
b05 b06 b07 b08 b09
b10 b11 b12 b13 b14
b15 :=
bm4 bm3 bm2 bm1
b00 b01 b02 b03 b04
b05 b06 b07 b08 b09
b10 b11 b12 b13 b14
b15 bm5 ;
ENDWHILE ;
QUIT " Program *xbessj0* XREF " .
|
