Loading and displaying a table (MATLAB)

With the MATLAB routines GetBoundaryTables_RWZ and EvaluateKernelFrequency we can load and display a specific table located in the directory TablesHeunRBC (see here) for the Regge-Wheeler equation with \(\rho_b = 60\), \(\ell=5\), \(M=1\). Using the following MATLAB code, we load the table into memory and evaluate it according to \( \widehat{\Xi}_{12}(s,120M) = \sum_{k=1}^{12} (\gamma_k/(2M))/(s - \beta_k/(2M))\).

% Settings
M=1;
ELL=5;
Rho = 60; % r/(2M)
RBC_Dir = 'Absolute_Path_Containing_TablesHeunRBC'; % Ex: /home/uname/Desktop/temp

% Load table
[kap,mu,qR,qI,mR,mI,nmu,nmR] = GetBoundaryTables_RWZ(Rho,ELL,M,'rw',RBC_Dir);

% Evaluate on imaginary axis 
y = linspace(-1,1,10000); % Laplace frequency taken to be s=i*y
[full_kernel,pair_kernel,sing_kernel] = EvaluateKernelFrequency(kap,mu,qR,qI,mR,mI,y);

% Plot
set(0, 'defaultaxesfontsize',16,'defaultaxeslinewidth',1.2,...
       'defaultlinelinewidth',1.2,'defaultpatchlinewidth',0.7,...
       'defaulttextfontsize',18,'defaultTextInterpreter', 'latex'); 
plot(y,real(full_kernel),'black')
hold on
plot(y,imag(full_kernel),'blue--')
xlabel('Im(s)')
title('RBC kernel: RW equation with $\ell = 5$ and $\rho = 60$','Interpreter', 'Latex')
legend({'Re $\widehat{\Xi}_{12}$','Im $\widehat{\Xi}_{12}$'},'Interpreter', 'Latex')