diff --git a/Eigen/src/IterativeLinearSolvers/ConjugateGradient.h b/Eigen/src/IterativeLinearSolvers/ConjugateGradient.h
index ced3e310c..35b89404a 100644
--- a/Eigen/src/IterativeLinearSolvers/ConjugateGradient.h
+++ b/Eigen/src/IterativeLinearSolvers/ConjugateGradient.h
@@ -120,7 +120,7 @@ struct traits<ConjugateGradient<_MatrixType,_UpLo,_Preconditioner> >
   * SparseMatrix<double> A(n,n);
   * // fill A and b
   * ConjugateGradient<SparseMatrix<double> > cg;
-  * cg(A);
+  * cg.compute(A);
   * x = cg.solve(b);
   * std::cout << "#iterations:     " << cg.iterations() << std::endl;
   * std::cout << "estimated error: " << cg.error()      << std::endl;