def METHOD_NAME(theta, theta0, ktheta):
    """Integrand for harmonic angle restraints. Domain is on [0,pi].
    Args:
        theta (float): Angle to be integrated, in radians
        theta0 (float): Equilibrium angle, in radians
        ktheta (float): Force constant, in kcal mol-1 rad-2
    Returns:
        float: Value of integrand
    """
    return sin(theta)*np.exp(-(ktheta*(theta-theta0)**2)/(2*R*T))