calculating of electric field distribution in 1D photonic crystal

I am trying to study the localization of light in crystals periodic and quasi-periodic photonic.

i have to simulate the distribution of the electric field E as a function of the thickness z of the periodic structure with the matrix transfer method.

Is there any matlab code for the Electric field intensity distribution of 1D photonic crystal?

and thank you in advance.