Created vectorized function to calculate P_ul

function p_ul = power_upper_limit(max_power,n_ps,prob)

%POWER_UPPER_LIMIT Upper limit on signal power at given prob for max_power

%measured

%   Detailed explanation goes here

    parser=inputParser;

    parser.KeepUnmatched=true;

    addRequired(parser,'prob',@(x) isreal(x) && x<=1);

    %p_ul = fzero(@(x) ncx2cdf(max_power,2*n_ps,x)-prob,max_power);

    p_ul = arrayfun(@(max_power) fzero(@(x) ncx2cdf(max_power,2*n_ps,x)-prob,max_power),max_power);

end