f=300:50;3400;
a=[295 376 484 582 661 733 794 851 902 955 1000];
b=[1035 1072 1109 1122 1109 1072 1035 1000];
c=[977 955 928 905 881 861 842 824];
d=[807 791 775 760 745 732 720 708];
e=[698 689 679 670 661 652 643 634];
f=[626 617 607 698 590 580 571 562];
g=[553 543 534 525];
l=[501 473 444 412];
h=[a b c d e f g];
i=.001 .* h;
m=.001 .*l;
mm=m.^2;
k=i.^2;
q=length(i);
x=sum(k);
y=sum(mm);
w=(5/310)*x+(1/31)*y
w0=(5/310)*q+4/31
z=10*log10(w/(w0))