From e623e21f2694f1176992424abe8793e38407d232 Mon Sep 17 00:00:00 2001
From: Denys Dutykh <dutykh@users.noreply.github.com>
Date: Tue, 20 Feb 2024 12:10:49 +0400
Subject: [PATCH] Polynomial eigenvalue problem solution in multiple precision

---
 findeigshp.m | 23 +++++++++++++++++++++++
 1 file changed, 23 insertions(+)
 create mode 100644 findeigshp.m

diff --git a/findeigshp.m b/findeigshp.m
new file mode 100644
index 0000000..dbda6ad
--- /dev/null
+++ b/findeigshp.m
@@ -0,0 +1,23 @@
+% The multi-precision version
+
+clear
+close all
+format longG
+
+addpath('../advanpix/'); % please set the path to advanpix repository
+
+mp.Digits(200);
+
+M0 = mp.read('M0_200.mat');
+M1 = mp.read('M1_200.mat');
+M2 = mp.read('M2_200.mat');
+
+e = polyeig(M0, mp('1i')*M1, M2);
+
+figure(1);
+plot(real(e), imag(e), 'bx'), grid off, hold off
+xlabel('$\mathrm{Re}\Omega$', 'interpreter', 'LaTeX', 'fontsize', 14);
+ylabel('$\mathrm{Im}\Omega$', 'interpreter', 'LaTeX', 'fontsize', 14);
+axis([-6 6 -5 1]);
+title('Worm hole QNMs');
+set(gcf, 'color', 'w');
\ No newline at end of file