function [ Pb ] = BudalLimit(S,f,V)

rho = 1025;
g = 9.81;

c0 = rho*pi*g/4;

cinf = 245;
  
df = f(2)-f(1);

Pow = 2*df*S.*(f.^2);
Pb = 2*V*c0*sqrt(sum(Pow));


end
