Optical fiber Matlab code TE modes

Hello,

I want to solve the this expression from weakly guiding modes in matlab. (optics, D. Marcuse). Can anyone help me with this. For the region r < a:

Combined two solutions with v+ 1 and í— 1 in order to achieve cancellation of the Ex and Hy components.

Thanks in advance.