function [dg g m] = delta_gain( pave);

g = [5 10 15 20 25 30 35 40 45 50 55];

for k = 1:length(pave)-1;
    dg(k) = 10*log10(pave(k+1)/pave(1) );
end

stem(g, dg); grid on

[m] = polyfit(g(1:5),dg(1:5),1)
hold all;
plot(g,g*m(1)+m(2));

xlabel('Commanded Gain (dB)')
ylabel('Actual Gain Relative 0 dB Gain Output(dB)')

return