a=polyfit(SNR1,10*log10(S1),3);
mymin=min(SNR1);
mymax=max(SNR1);
xx=linspace(mymin,mymax,100);
yy=polyval(a,xx);
plot(SNR1,10*log10(S1),'o',xx,yy);
hold on

a=polyfit(SNR2,10*log10(S2),3);
mymin=min(SNR2);
mymax=max(SNR2);
xx=linspace(mymin,mymax,100);
yy=polyval(a,xx);
plot(SNR2,10*log10(S2),'o',xx,yy);
hold on

a=polyfit(SNR3,10*log10(S3),3);
mymin=min(SNR3);
mymax=max(SNR3);
xx=linspace(mymin,mymax,100);
yy=polyval(a,xx);
plot(SNR3,10*log10(S3),'o',xx,yy);
hold on

a=polyfit(SNR4,10*log10(S4),3);
mymin=min(SNR4);
mymax=max(SNR4);
xx=linspace(mymin,mymax,100);
yy=polyval(a,xx);
plot(SNR4,10*log10(S4),'o',xx,yy);
hold on

a=polyfit(SNR5,10*log10(S5),3);
mymin=min(SNR5);
mymax=max(SNR5);
xx=linspace(mymin,mymax,100);
yy=polyval(a,xx);
plot(SNR5,10*log10(S5),'o',xx,yy);
hold on

a=polyfit(SNR6,10*log10(S6),3);
mymin=min(SNR6);
mymax=max(SNR6);
xx=linspace(mymin,mymax,100);
yy=polyval(a,xx);
plot(SNR6,10*log10(S6),'o',xx,yy);
hold on

a=polyfit(SNR7,10*log10(S7),3);
mymin=min(SNR7);
mymax=max(SNR7);
xx=linspace(mymin,mymax,100);
yy=polyval(a,xx);
plot(SNR7,10*log10(S7),'o',xx,yy);
hold on


a=polyfit(SNR,10*log10(Shannon),3);
mymin=min(SNR);
mymax=max(SNR);
xx=linspace(mymin,mymax,100);
yy=polyval(a,xx);
plot(SNR,10*log10(Shannon),'o',xx,yy);
hold on
%axis([2,8,-120,-50])