Max freq response

idmlin.py

@@ -55,6 +55,10 @@ class IDMLin:
         # w,v = sparse.linalg.eigs(a)
         return all(real(w)<1e-10)
 
+    def maxresponse(self):
+        _, mag, _ = self.sys().bode()
+        return max(mag)
+
     def os(self):
         f,g,h = self.f, self.g, self.h
         return 4*f-(g+h)**2