% 曲線擬合程序
function [xi,yi]=dpoly(n)
x=0:0.1:1;
y=[-.477 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.30 11.2];
p=polyfit(x,y,n);
xi=linspace(0,1,100);
yi=polyval(p,xi);
plot(x,y,'-o',xi,yi,'--');
xlabel('x');
ylabel('y=f(x)');
title('n order curve fitting');
end