1 | test_besselj = ->
|
2 | run_test [
|
3 |
|
4 | "besselj(x,n)",
|
5 | "besselj(x,n)",
|
6 |
|
7 | "besselj(0,0)",
|
8 | "1",
|
9 |
|
10 | "besselj(0,1)",
|
11 | "0",
|
12 |
|
13 | "besselj(0,-1)",
|
14 | "0",
|
15 |
|
16 | "besselj(x,1/2)-sqrt(2/pi/x)*sin(x)",
|
17 | "0",
|
18 |
|
19 | "besselj(x,-1/2)-sqrt(2/pi/x)*cos(x)",
|
20 | "0",
|
21 |
|
22 | "besselj(x,3/2)-sqrt(2/pi/x)*(sin(x)/x-cos(x))",
|
23 | "0",
|
24 |
|
25 | "besselj(x,-3/2)-sqrt(2/pi/x)*(-cos(x)/x-sin(x))",
|
26 | "0",
|
27 |
|
28 | "besselj(x,5/2)-sqrt(2/pi/x)*((3/x^2-1)*sin(x)-3/x*cos(x))",
|
29 | "0",
|
30 |
|
31 | "besselj(x,-5/2)-sqrt(2/pi/x)*((3/x^2-1)*cos(x)+3/x*sin(x))",
|
32 | "0",
|
33 |
|
34 |
|
35 |
|
36 | "besselj(x,3/2)-(1/x)*besselj(x,1/2)+besselj(x,-1/2)",
|
37 | "0",
|
38 |
|
39 | "besselj(x,-3/2)+(1/x)*besselj(x,-1/2)+besselj(x,1/2)",
|
40 | "0",
|
41 |
|
42 |
|
43 |
|
44 | "y=besselj(x,5/2)",
|
45 | "",
|
46 |
|
47 | "x^2*d(y,x,x)+x*d(y,x)+(x^2-(5/2)^2)*y",
|
48 | "0",
|
49 |
|
50 | "y=quote(y)",
|
51 | "",
|
52 | ]
|
53 |
|
54 |
|