disp(' ');
disp(' add_dB.m  ver 1.0  February 9, 2010');
disp(' by Tom Irvine ');
%
disp(' ');
disp(' This script assumes that the dB values represent ');
disp(' a sound pressure level or a power spectral density ');
%
clear a;
clear sum;
clear length;
clear C;
%
disp(' ');
disp(' Enter vector name of dB values ');
a=input(' ');
n=length(a);
%
C=10;
sum=0;
for(i=1:n)
  sum=sum+10^(a(i)/C);
end
sum=C*log10(sum);
%
out1=sprintf('\n overall level = %8.4g dB \n',sum);
disp(out1);