clearclcX=-1011;Y=125X^2+1;dx0= 00739644970414201;dxn=-00739644970414201;S=csfitXYdx0dxn;x1=-1001-05;y1=polyvalS1x1-X 1;x2=-050010;y2=polyvalS2x2-X2;x3=000105; y3=polyvalS3x3-X3;x4=050011;y4=polyvalS4

clear,clc：清空命令窗口和工作空间中的变量。

X=-1:0.1:1：定义一个行向量X，包含从-1到1，间隔为0.1的所有数。

Y=1./(25*X.^2+1)：定义一个行向量Y，每个元素为对应X元素的函数值。

dx0= 0.0739644970414201;dxn=-0.0739644970414201：定义两个常数，用于限定边界条件。

S=csfit(X,Y,dx0,dxn)：使用csfit函数拟合X和Y，限定边界条件为dx0和dxn，返回拟合结果。

x1=-1:0.01:-0.5;y1=polyval(S(1,:),x1-X (1))：定义一个行向量x1，包含从-1到-0.5，间隔为0.01的所有数；定义一个行向量y1，每个元素为拟合结果S的第一行多项式在对应x1元素处的函数值。

x2=-0.5:0.01:0;y2=polyval(S(2,:),x2-X(2))：定义一个行向量x2，包含从-0.5到0，间隔为0.01的所有数；定义一个行向量y2，每个元素为拟合结果S的第二行多项式在对应x2元素处的函数值。

x3=0:0.01:0.5; y3=polyval(S(3,:),x3-X(3))：定义一个行向量x3，包含从0到0.5，间隔为0.01的所有数；定义一个行向量y3，每个元素为拟合结果S的第三行多项式在对应x3元素处的函数值。

x4=0.5:0.01:1;y4=polyval(S(4,:),x4-X(4))：定义一个行向量x4，包含从0.5到1，间隔为0.01的所有数；定义一个行向量y4，每个元素为拟合结果S的第四行多项式在对应x4元素处的函数值。

plot(x1,y1,x2,y2,x3,y3,x4,y4, X,Y,'.')：绘制多个曲线和散点图，其中x1和y1为一组，x2和y2为一组，x3和y3为一组，x4和y4为一组，X和Y为一组，用不同颜色区分不同的曲线和散点图