Dear both.
Thanks for your responses. Daniel, it really helps to speed this up, thanks for that. Do you think if there is any way to speed up the code even further, maybe by manipulating the module? I recently found out that Mathematica seems to be really slow if arbitary precision numbers are combined with Bessel functions, which I need to use for another waveguide problem. Is there any solution to this?
Thanks again, Markus