Hello Ray, thanks for your reply.
Hmmm, maybe my question was not clear enough.
Anyway, unfortunately I think you're misunderstanding the
situation. Let me give you another couple of examples to
see what's wrong.
(C7) j0(x)-bessel_j(0, x);
(D7) 0
This means that, on algebraic level, the first parameter in
bessel_j(v, z) is interpreted as the order and the second
is the argument (which is consistent with the documentation).
Then, when evaluating to the floating point, the same thing
is interpreted as the argument.
(C10) float(j0(1)-bessel_j(0,1));
(D10) .7651976865579666
(C11) j0(1)-j1(0);
(D11) .7651976865579666
I hope this helps you to see the problem.
Best regards,
Keita