def METHOD_NAME(eu, om):
    sph = math.sin(eu[e_phi])
    cph = math.cos(eu[e_phi])
    tth = math.tan(eu[e_theta])
    cth = math.cos(eu[e_theta])
    return [om[r_p] + sph * tth * om[r_q] + cph * tth * om[r_r],
            cph * om[r_q] - sph * om[r_r],
            sph / cth * om[r_q] + cph / cth * om[r_r]]