diff --git a/README.md b/README.md
index 3cabf97..ceac784 100644
--- a/README.md
+++ b/README.md
@@ -12,10 +12,9 @@ function (t,u) ↦ 1 using Adams' method.
 
 ```rust
 use sundials::CVode;
-let f = |t, u: &[f64; 1], du: &mut [f64; 1]| { *du = [1.] };
-let ode = CVode::adams(f, 0., &[0.])?;
-let mut u1 = [f64::NAN];
-ode.solve(1., &mut u1);
+let mut ode = CVode::adams(0., &[0.], |t, u, du| *du = [1.])?;
+let (u1, _) = ode.solution(1.);
+assert_eq!(u1[0], 1.);
 ```
 
 Then `u[0]` contains the solution u at time t=1 and `u[1]` the
diff --git a/src/lib.rs b/src/lib.rs
index 02ecc1c..9ee2011 100644
--- a/src/lib.rs
+++ b/src/lib.rs
@@ -13,8 +13,7 @@
 //! use sundials::CVode;
 //! # fn main() -> Result<(), Box<dyn std::error::Error>> {
 //! let mut ode = CVode::adams(0., &[0.], |t, u, du| *du = [1.])?;
-//! let mut u1 = [f64::NAN];
-//! ode.solve(1., &mut u1);
+//! let (u1, _) = ode.solution(1.);
 //! assert_eq!(u1[0], 1.);
 //! # Ok(()) }
 //! ```