How to compare approximations and near exact solutions of SDE with Euler scheme?

Question

The code I used to apply the SDEs and the Black-Scholes scheme in python:

# dw Wiener process
def dW(delta_t):
    """" Random sample