infradius_hnorm.m

%{
%{
 % @file infradius_hnorm.m
 % @brief Helper function for infinite radius computation with normal CDF transformation.
%}
%}
 
%{
%{
 % @brief Helper function for infinite radius computation with normal CDF transformation.
 % @fn infradius_hnorm(x, L, N, alpha)
 % @param x Normal CDF transformation parameters.
 % @param L Service demand matrix.
 % @param N Population vector.
 % @param alpha Load-dependent rate matrix.
 % @return y Evaluated function value for integration.
%}
%}
function y=infradius_hnorm(x,L,N,alpha)
M = size(L,1);
Nt = sum(N);
beta = N/Nt;
t = normcdf(x,0,1);
tb = sum(beta.*t,2);
h = @(x) pfqn_gld(sum(L.*repmat(exp(2*pi*1i*(t-tb)),M,1),2),Nt,alpha) * prod(normpdf(x,0,1),2);
 
y = zeros(size(x,1),1);
for i=1:size(x,1)
    y(i) = real(h(x(i,:)));
end
end