slab waveguide matlab code