diff --git a/Project.toml b/Project.toml
index 64f6867..5a7e4eb 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "DiffEqProblemLibrary"
 uuid = "a077e3f3-b75c-5d7f-a0c6-6bc4c8ec64a9"
 authors = ["Chris Rackauckas <accounts@chrisrackauckas.com>"]
-version = "4.8.0"
+version = "4.8.1"
 
 [deps]
 DiffEqBase = "2b5f629d-d688-5b77-993f-72d75c75574e"
diff --git a/src/ode/ode_simple_nonlinear_prob.jl b/src/ode/ode_simple_nonlinear_prob.jl
index 740b2a3..2028f3b 100644
--- a/src/ode/ode_simple_nonlinear_prob.jl
+++ b/src/ode/ode_simple_nonlinear_prob.jl
@@ -55,7 +55,7 @@ prob_ode_fitzhughnagumo = ODEProblem(fitz,[1.0;1.0],(0.0,1.0),(0.7,0.8,1/12.5,0.
 eqs = [D(x) ~ μ*((1-x^2)*y - x),
        D(y) ~ y]
 de = ODESystem(eqs)
-van = ODEFunction(de, [x,y], [μ], jac=true, Wfact=true)
+van = ODEFunction(de, [x,y], [μ], jac=true)
 
 """
 Van der Pol Equations
@@ -97,7 +97,7 @@ eqs = [D(y₁) ~ -k₁*y₁+k₃*y₂*y₃,
        D(y₂) ~ k₁*y₁-k₂*y₂^2-k₃*y₂*y₃,
        D(y₃) ~ k₂*y₂^2]
 de = ODESystem(eqs)
-rober = ODEFunction(de, [y₁,y₂,y₃], [k₁,k₂,k₃], jac=true, Wfact=true)
+rober = ODEFunction(de, [y₁,y₂,y₃], [k₁,k₂,k₃], jac=true)
 
 """
 The Robertson biochemical reactions: (Stiff)
@@ -174,7 +174,7 @@ eqs = [D(y₁) ~ I₁*y₂*y₃,
        D(y₂) ~ I₂*y₁*y₃,
        D(y₃) ~ I₃*y₁*y₂]
 de = ODESystem(eqs)
-rigid = ODEFunction(de, [y₁,y₂,y₃], [I₁,I₂,I₃], jac=true, Wfact=true)
+rigid = ODEFunction(de, [y₁,y₂,y₃], [I₁,I₂,I₃], jac=true)
 
 """
 Rigid Body Equations (Non-stiff)
@@ -334,7 +334,7 @@ eqs = [D(y1) ~ -p1*y1 + p2*y2 + p3*y3 + p4,
 de = ODESystem(eqs)
 hires = ODEFunction(de, [y1,y2,y3,y4,y5,y6,y7,y8],
                         [p1,p2,p3,p4,p5,p6,p7,p8,p9,p10,p11,p12],
-                        jac=true, Wfact=true)
+                        jac=true)
 
 u0 = zeros(8)
 u0[1] = 1
@@ -374,8 +374,7 @@ eqs = [D(y1) ~ p1*(y2+y1*(1-p2*y1-y2)),
        D(y3) ~ p3*(y1-y3)]
 de = ODESystem(eqs)
 jac = calculate_jacobian(de)
-ModelingToolkit.calculate_factorized_W(de)
-orego = ODEFunction(de, [y1,y2,y3], [p1,p2,p3], jac=true, Wfact=true)
+orego = ODEFunction(de, [y1,y2,y3], [p1,p2,p3], jac=true)
 
 """
 Orego Problem (Stiff)