using numerical techniques plot the intensity

Question

Using numerical techniques, plot the intensity distribution for LP₀₁ and LP₁₁ modes. Vary the value of V (reduce it to 1.5) and plot several of the intensity plots (of LP₀₁) to demonstrate the effects of V on the evanescent field. Show that at V=2.405, the core contains about 83% of the power in the fiber for the LP₀₁ mode. You must use a 3-D plot to demonstrate the effect of increase in cladding.

I will need the matlab codes for the plot.