diff --git a/examples/create-a-model-with-coordinates.ipynb b/examples/create-a-model-with-coordinates.ipynb
index f2b12eed..e84c21b9 100644
--- a/examples/create-a-model-with-coordinates.ipynb
+++ b/examples/create-a-model-with-coordinates.ipynb
@@ -150,7 +150,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "m.solve(solver_name=\"highs\")"
+    "m.solve(solver_name=\"highs\", output_flag=False)"
    ]
   },
   {
diff --git a/examples/create-a-model.ipynb b/examples/create-a-model.ipynb
index a158e0cf..b6fc9705 100644
--- a/examples/create-a-model.ipynb
+++ b/examples/create-a-model.ipynb
@@ -215,7 +215,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "m.solve(solver_name=\"highs\")"
+    "m.solve(solver_name=\"highs\", output_flag=False)"
    ]
   },
   {
diff --git a/examples/manipulating-models.ipynb b/examples/manipulating-models.ipynb
index 81106ab3..6903386b 100644
--- a/examples/manipulating-models.ipynb
+++ b/examples/manipulating-models.ipynb
@@ -48,9 +48,9 @@
     "con2 = m.add_constraints(5 * x + 2 * y >= 3 * factor, name=\"con2\")\n",
     "\n",
     "m.add_objective(x + 2 * y)\n",
-    "m.solve(solver_name=\"highs\")\n",
+    "m.solve(solver_name=\"highs\", output_flag=False)\n",
     "\n",
-    "m.solve(solver_name=\"highs\")\n",
+    "m.solve(solver_name=\"highs\", output_flag=False)\n",
     "sol = m.solution.to_dataframe()\n",
     "sol.plot(grid=True, ylabel=\"Optimal Value\")"
    ]
@@ -95,7 +95,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "m.solve(solver_name=\"highs\")\n",
+    "m.solve(solver_name=\"highs\", output_flag=False)\n",
     "sol = m.solution.to_dataframe()\n",
     "sol.plot(grid=True, ylabel=\"Optimal Value\")"
    ]
@@ -137,7 +137,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "m.solve(solver_name=\"highs\")\n",
+    "m.solve(solver_name=\"highs\", output_flag=False)\n",
     "sol = m.solution.to_dataframe()\n",
     "sol.plot(grid=True, ylabel=\"Optimal Value\")"
    ]
@@ -190,7 +190,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "m.solve(solver_name=\"highs\")\n",
+    "m.solve(solver_name=\"highs\", output_flag=False)\n",
     "sol = m.solution.to_dataframe()\n",
     "sol.plot(grid=True, ylabel=\"Optimal Value\")"
    ]
@@ -242,7 +242,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "m.solve(solver_name=\"highs\")\n",
+    "m.solve(solver_name=\"highs\", output_flag=False)\n",
     "sol = m.solution.to_dataframe()\n",
     "sol.plot(grid=True, ylabel=\"Optimal Value\")"
    ]
@@ -276,7 +276,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "m.solve(solver_name=\"highs\")\n",
+    "m.solve(solver_name=\"highs\", output_flag=False)\n",
     "sol = m.solution.to_dataframe()\n",
     "sol.plot(grid=True, ylabel=\"Optimal Value\")"
    ]
@@ -326,7 +326,7 @@
     "# Penalize activation of z in the objective\n",
     "m.objective = x + 3 * y + 10 * z\n",
     "\n",
-    "m.solve(solver_name=\"highs\")"
+    "m.solve(solver_name=\"highs\", output_flag=False)"
    ]
   },
   {
@@ -346,7 +346,7 @@
    "source": [
     "m.variables.binaries.fix()\n",
     "m.variables.binaries.relax()\n",
-    "m.solve(solver_name=\"highs\")\n",
+    "m.solve(solver_name=\"highs\", output_flag=False)\n",
     "\n",
     "# Dual values are now available on the constraints\n",
     "m.constraints[\"con1\"].dual"
diff --git a/examples/transport-tutorial.ipynb b/examples/transport-tutorial.ipynb
index b42e67e8..cd5cdccd 100644
--- a/examples/transport-tutorial.ipynb
+++ b/examples/transport-tutorial.ipynb
@@ -417,7 +417,7 @@
    "outputs": [],
    "source": [
     "# Solve the model\n",
-    "m.solve(solver_name=\"highs\")"
+    "m.solve(solver_name=\"highs\", output_flag=False)"
    ]
   },
   {