Hello Pavel.

I already used this workaround and it is fine for me. Just wanted to let you know that there is a problem. Maybe you could mention it in the documentation.

Best regards,

Michael

Thank you for your answer.
I understand the problem, but it seems to me that it never had the problem in Matlab double-precision. Maybe the implemented some tolerance in the "==" operator. You can make the question public if you want. :)

Thank you Pavel. It will use bessel functions with arbitrary order and real or imaginary argument.

In general I only use quadruple precision. I will need the bessely function again, but it will have non-integer order. But it is a great improvement for me at this point already.

I used quadruple precision. It is probably even faster. Problem is, I don't have the previous version of the toolbox anymore to make more comparisons, but I think it is between 10-15 times faster. :)

Dear Pavel,

I got a speed improvement of factor 10-15. Thank you very much.