检查下列代码并改正function S = csfitX Y dx0 dxnN = lengthX-1;H = diffX;D = diffYH;A = H1N-2;B = 2H1N-1+H2N;C = H2N-1;U = 6diffD;B1 = B1 - H12;U1 = U1 - 3D1;BN-1 = BN-1 - HN2;UN-1 = UN-1 - 3-DN;for k = 2N-1 - 教育

在第二行的分号后面多了一个空格，应该删除。
在第三行的分号后面多了一个中文分号，应该改为英文分号。
在第四行的分号后面多了一个空格，应该删除。
在第五行的分号后面多了一个空格，应该删除。
在第六行的分号后面多了一个空格，应该删除。
在第七行的分号后面多了一个空格，应该删除。
在第八行的分号后面多了一个空格，应该删除。
在第九行的分号后面多了一个空格，应该删除。
在第十行的分号后面多了一个空格，应该删除。
在第十一行的分号后面多了一个空格，应该删除。
在第十二行的分号后面多了一个空格，应该删除。
在第十三行的分号后面多了一个空格，应该删除。
在第十四行的分号后面多了一个空格，应该删除。
在第十五行的分号后面多了一个空格，应该删除。
在第十七行的变量M没有初始化，应该在第八行添加M=zeros(N,1)。
在第十七行的下标越界，应该改为M(N)=U(N-1)/B(N-1)。
在第十九行的变量M没有初始化，应该在第八行添加M=zeros(N,1)。
在第二十行的下标越界，应该改为M(k+1) = (U(k) - C(k)*M(k+2))/B(k+1)。
在第二十二行的变量M没有初始化，应该在第八行添加M=zeros(N,1)。
在第二十三行的下标越界，应该改为M(1) = 3*(D(1)-dx0)/H(1) - M(3)/2;。
在第二十四行的下标越界，应该改为M(N+1) = 3*(dxn-D(N))/H(N) - M(N-1)/2;。
在第二十五行的变量S应该改为S=zeros(N-1,4)。
在第三十行的字符串中，单引号和双引号不匹配，应该改为plot(x1,y1,x2,y2,x3,y3,x4, y4, 'o', X,Y,'.')