x=[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30];
y=[67 88 80 68 31 27 31 26 16 20 6 9 2 5 12 8 3 4 1 5 7 4 1 5 7 3 4 8 6 6];
x0=[32];
y1=[8];
for i=2:5
p=polyfit(x,y,5);
y0=polyval(p,x0)
%plot(x0,y0,'ro',x0,y1,'b+')