diff --git a/examples/1protationsymmetry/doc/main.md b/examples/1protationsymmetry/doc/main.md
index ebb097d694f29f2e3b07f8fc8879dfe7e2c450aa..8fafcb7edbb6aa2fa92f31e10f592cd49f25f366 100644
--- a/examples/1protationsymmetry/doc/main.md
+++ b/examples/1protationsymmetry/doc/main.md
@@ -36,7 +36,7 @@ and compute the convergence rates.
 #include <dumux/common/parameters.hh> // for getParam
 #include <dumux/common/integrate.hh>  // for integrateL2Error
 
-#include <dumux/linear/seqsolverbackend.hh> // for UMFPackBackend
+#include <dumux/linear/seqsolverbackend.hh> // for ILU0BiCGSTABBackend
 #include <dumux/linear/pdesolver.hh>        // for LinearPDESolver
 #include <dumux/assembly/fvassembler.hh>
 #include <dumux/assembly/diffmethod.hh>
@@ -133,7 +133,7 @@ solution and stores the result therein.
     using Assembler = FVAssembler<TypeTag, DiffMethod::analytic>;
     auto assembler = std::make_shared<Assembler>(problem, gridGeometry, gridVariables);
 
-    using LinearSolver = UMFPackBackend;
+    using LinearSolver = ILU0BiCGSTABBackend;
     auto linearSolver = std::make_shared<LinearSolver>();
     LinearPDESolver<Assembler, LinearSolver> solver(assembler,  linearSolver);
     solver.setVerbose(false); // suppress output during solve()
diff --git a/examples/1protationsymmetry/main.cc b/examples/1protationsymmetry/main.cc
index 6db3dcd673af56b49e1ff26f77c30b423f2431e3..e25029a781aa041bb83b2b8be4f6596236b009d3 100644
--- a/examples/1protationsymmetry/main.cc
+++ b/examples/1protationsymmetry/main.cc
@@ -33,7 +33,7 @@
 #include <dumux/common/parameters.hh> // for getParam
 #include <dumux/common/integrate.hh>  // for integrateL2Error
 
-#include <dumux/linear/seqsolverbackend.hh> // for UMFPackBackend
+#include <dumux/linear/seqsolverbackend.hh> // for ILU0BiCGSTABBackend
 #include <dumux/linear/pdesolver.hh>        // for LinearPDESolver
 #include <dumux/assembly/fvassembler.hh>
 #include <dumux/assembly/diffmethod.hh>
@@ -122,7 +122,7 @@ int main(int argc, char** argv) try
     using Assembler = FVAssembler<TypeTag, DiffMethod::analytic>;
     auto assembler = std::make_shared<Assembler>(problem, gridGeometry, gridVariables);
 
-    using LinearSolver = UMFPackBackend;
+    using LinearSolver = ILU0BiCGSTABBackend;
     auto linearSolver = std::make_shared<LinearSolver>();
     LinearPDESolver<Assembler, LinearSolver> solver(assembler,  linearSolver);
     solver.setVerbose(false); // suppress output during solve()