w=[50 100 150 200 250 300 350 400 450 500 550 600 650 700 750 800 850 900 950 1000 1050 1100 1150 1200 1250 1300 1350 1400 1450 1500 1550 1600 1650 1700 1750 1800 1850 1900 1950 2000 2050 2100 2150 2200 2250 2300 2350 2400 2450 2500 2550 2600 2650 2700 2750 2800 2850 2900 2950 3000 3100 3200 3300 3400 3500 3600 3700 3800 3900 4000 4100 4200 4300 4400 4500 4600 4700 4800 4900 5000];
p=[.71 8.91 35.5 89.1 178 295 376 484 582 661 733 794 851 902 955 1000 1035 1072 1109 1122 1109 1072 1035 1000 977 955 928 905 881 861 842 824 807 791 775 760 745 732 720 708 698 689 679 670 661 652 643 634 626 617 607 598 590 580 571 562 553 543 534 525 501 473 444 412 376 335 292 251 214 178 144.5 116 92.3 72.4 56.2 43.7 33.9 26.3 20.4 15.9];
p1=p.^2;
w1=w(1,6:75);p2=p(1,6:75);p3=p2.^2;
s=0
for i=1:length(w1)-1
s=s+(w1(i+1)-w1(i))*(p3(i)+p3(i+1))*.5;
end
s1=(1122^2)*3100;
r=10*log10(s1/s);
 s2=1024^2*3100;
r1=10*log10(s2/s)
s=0
for i=1:length(w)-1
s=s+(w(i+1)-w(i))*(p1(i)+p1(i+1))*.5;
end
r3=10*log10(s2/s)