diff --git a/src/DistQuads.jl b/src/DistQuads.jl
index 818a166..e7c28f8 100644
--- a/src/DistQuads.jl
+++ b/src/DistQuads.jl
@@ -17,7 +17,7 @@ distribution(dq::DistQuad) = dq.d
 
 DistQuad(d; N = 32) = DistQuad(d, N)
 function DistQuad(d::Distributions.Beta, N)
-    gj = gaussjacobi(N, d.α-1, d.β-1)
+    gj = FastGaussQuadrature.JacobiRec(N, d.α-1, d.β-1)
     G = gamma(d.α)*gamma(d.β)/gamma(d.α+d.β)
     w = gj[2]/((2.0^(d.α+d.β-1.0))*G)
     n = (1.-gj[1])./2