diff --git a/docs/src/example.md b/docs/src/example.md
index 29caa2a..5a18074 100644
--- a/docs/src/example.md
+++ b/docs/src/example.md
@@ -51,7 +51,7 @@ function simulate_with_tunables(tunables)
   coeffs = reshape(tunables[6:end], size(p.coeffs))
   newp = Parameters(subpars, coeffs)
   newprob = remake(prob; p = newp)
-  sol = solve(prob, Tsit5())
+  sol = solve(newprob, Tsit5())
   return sum(sol.u[end])
 end
 ```