Part B from problem 1
It can be seen from figures 5-7 that our analytical values for \(\beta\) agree very well with the simulation.
Part C from problem 1
These figures also plot radial dependence with r on the bottom axis. This looks very close to the ones I found in mathematica.
Part D from problem 1
I dont know where in the code we can find the spot size let alone to graph it with respect to wavelength. For the sake of time I had to move on.