Subr_ComputeIgMLEParams(G)

function [Mu,Lambda,Theta] = Subr_ComputeIgMLEParams(G)
%
%
global M N
%
Lambda = zeros(1,M);
Mu = zeros(1,M);
Theta = zeros(1,M);
%
for ii = 1:M
    
    z = zeros(1,N);
    z(1,1:N) = G(ii,:);
    
    %   Compute Mu
    Mu(1,ii) = mean(z);                        % compute mean
    
    %   Compute Lambda
    InvSum = sum(1./z) - N/mean(z);
    Lambda(1,ii) = N/InvSum;
    
    %   Compute Theta
    Theta(1,ii) = Lambda(1,ii)/Mu(1,ii);
    clear  InvSum z
    
end