import numpy as np

def rand_quat():
    q = np.random.normal(0.0, 1.0, 4)
    q = q / np.linalg.norm(q)
    return q