This post is my first in this Mathematica community. Thanks to all in advance who contribute to this site's success and to its users' learning experiences.

My problem is such that plotting a function is taking a painfully long time (and I have multiple plots for similar functions). I am hoping someone can offer suggestions to speed things up. Currently the plot takes more than 1 hour to finish. The image below is clipped from Mathematica.

Fk = Re[HankelH2[1, k]/(HankelH2[1, k] + I HankelH2[0, k])]

Phi = 2/Pi Integral from 0 to Infinity [Fk/k Sin[k s]] dk

I was thinking about using a Taylor Series expansion for Sin[ks], which would sacrifice exactness for speed. However my gut tells me that it's the Hankel funtions that are slowing down the plot calculations. What can I do?

Many thanks!

This post is my first in this Mathematica community. Thanks to all in advance who contribute to this site's success and to its users' learning experiences.

My problem is such that plotting a function is taking a painfully long time (and I have multiple plots for similar functions). I am hoping someone can offer suggestions to speed things up. Currently the plot takes more than 1 hour to finish. The image below is clipped from Mathematica.

Fk = Re[HankelH2[1, k]/(HankelH2[1, k] + I HankelH2[0, k])]

I was thinking about using a Taylor Series expansion for Sin[k s], which would sacrifice exactness for speed. However my gut tells me that it's the Hankel funtions that are slowing down the plot calculations. What can I do?