If I’m not mistaken, GitHub - SciML/Integrals.jl: A common interface for quadrature and numerical integration for the SciML scientific machine learning organization already provides such rules on top of QuadGK functionality. See also Array issue with optimzing a function with quadgk for discussion on some other differentiable alternatives (N.B. that Quadrature.jl → Integrals.jl).