Give the Euler-Milstein approximation scheme for the following SDE:

\[d S_{t}=\mu_{S} S_{t} d t+\sigma S_{t}^{\beta} d B_{t}\]

where \(\beta \in(0,1]\). Generate five paths and plot them for the following parameter values:

\[\mu_{S}=0.1, \quad \sigma=0.7, \quad \beta=\frac{1}{3}, \quad S_{0}=100\]