def METHOD_NAME() -> tf.Tensor:
    mu: AnyNDArray = np.array([1.0, 1.3])
    var: AnyNDArray = np.array([3.0, 3.5])
    num_gauss_hermite_points = 25
    return quadrature.ndiagquad(
        [lambda *X: tf.exp(X[0])], num_gauss_hermite_points, [mu], [var]
    )