Pkmi_Calculate.m

k_max=100;
m_max=100;
global Pkmi_array
Pkmi_array=zeros(k_max+2,m_max,2*m_max+1);
for k=-1:k_max
    k
    for m=max([1,ceil(k/2)]):m_max
        for i=-m:m
            Pkmi_array(k+2,m,i+m+1)=Pkmi_closed_form(k,m,i);
        end
    end
end

save Pkmi_array Pkmi_array