Added vectorized RNG, allowed euler_sde to work with normal rhs functions