| 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 |
|