replace removed np.trapz by sp.trapezoid

bnavigator · bnavigator · commit 843d747646d4 · 2026-02-01T19:30:51.000+01:00
diff --git a/examples/cds110-L8a_maglev-limits.ipynb b/examples/cds110-L8a_maglev-limits.ipynb
@@ -253,7 +253,7 @@
     "omega = np.linspace(0, 1e6, 100000)\n",
     "for name, sys in zip(['C1', 'C2', 'C3'], [magS1, magS2, magS3]):\n",
     "    freqresp = ct.frequency_response(sys, omega)\n",
-    "    bodeint = np.trapz(np.log(freqresp.magnitude), omega)\n",
+    "    bodeint = sp.trapezoid(np.log(freqresp.magnitude), omega)\n",
     "    print(\"Bode integral for\", name, \"=\", bodeint)\n",
     "\n",
     "print(\"pi * sum[ Re(pk) ]\", pi * np.sum(magP.poles()[magP.poles().real > 0]))"
