function [p] = p_wu_integral(lam,d,alp,a,P,T)
%P_WU_INTEGRAL Summary of this function goes here
%   Detailed explanation goes here
    [~,~, p, ~, ~] = wake_bounds(lam,d,alp,a,P,T);
end