Hello,

Just wanted to ask if anyone got this issue solved. I am trying to create something like the below example and would like to know how to do it in bokeh.

"

r = np.arange(0, 2*a ,0.01)

theta= np.linspace(0, 2 * np.pi ,100)

R, Theta = np.meshgrid(r, theta)

A,B = R*np.cos(Theta),R*np.sin(Theta)

def I(l,r,Theta):

Z=np.zeros(R.shape,np.float64)

Z[np.where(R<=a)] = jv(l, (x[m] * R[np.where(R<=a)]/ a)) * np.cos(l_Theta[np.where(R<=a)])

Z[np.where(R>a)] = jv(l,x[m]) / kv(l, t[m]) * kv(l,(t[m] * R[np.where(R>a)]/ a)) * np.cos(l_Theta[np.where(R>a)])

return Z

plt.pcolormesh(A, B, I(l,r,Theta))

plt.show()

"

this is the result:

￼Cheers