The root finding somewhat mysteriously always seems to find the unstable fix point....

517f065b · Frank Hellmann · 317d82bc · 517f065b
Commit 517f065b authored Mar 15, 2017 by Frank Hellmann
--- a/PSD/complex_current_and_nodes.py
+++ b/PSD/complex_current_and_nodes.py
@@ -208,7 +208,6 @@ if __name__ == "__main__":
    root_rhs = define_root_rhs_omega(2*len(Y), rhs)

    guess = np.ones(7, dtype=np.float64) * 0.5
-    guess[-1] = 0.

    print(root_rhs(guess))