Monte Carlo simulation for SABR


I want to construct the SABR Implied Volatility across differerent strikes (hence the vol smile) using the Euler Scheme on spot and volatility as described in equation 3.1 here:
or page 12 here:

Let's say I have the SABR parameters (for instance T=1, alpha=25%, Beta=0.5, rho = 0.5, spot = 100) given. Can anyone please explain and hopefuly provide a "step-by-step" guide on how to construct the SABR volatility smile/skew using this Monte Carlo method.

I have implemented the algorithm itself (it's pretty simple) but don't know how to use it.