数值分析 最小二乘拟合

function charpt3
x=-1:0.5:2;
y=[-4.447 -0.452 0.551 0.048 -0.447 0.549 4.552];
n=3;               %三次拟合多项式
p=polyfit（xyn）;  %三次拟合多项式的系数
y1=polyval（px）;   %拟合多项式条件下给定节点处的函数值
r=（y1-y）*（y1-y）‘; %平方误差 y是实验所得函数值  
x0=-1.0:0.01:2.0;
y0=polyv