E field in the Sphere along z-axis from incident plane calculation by mie absorption

Hello,

I am studying FDTD with Sullivan's book.

I have a problem on calculating bessel function expansion calculation along the main axis of a dielectric sphere in fig 4.7.
I used matlab-funtions in "Matlab functions for mie scattering and absorption" to get mie coefficients and calculated internal field by the formula in page 8.

However, I fail to get the same graph in Sullivan's book.

Could anyone let me know how to get the E field in the Sphere along z-axis from incident plane wave?

I attached my code and Matlab functions for mie.

Thank you.