偏最小二乘Matlab源代码求助-经管之家官网！

本人下载了一份偏最小二乘Matlab源代码，但运行中程序报错，求助各位高手帮助排解问题！！！！
function [t,q,w,wh,f0,FF]=fun717(px,py,C)
%px自变量的输入个数
%py输入因变量的个数。
% C输入的自变量和因变量组成的矩阵
%t提取的主成分
%q为回归系数。
%w最大特征值所对应的特征向量。
%wh处理后的特征向量
%f0回归的标准化的方程系数
%FF原始变量的回归方程的系数
c=norm1(C); %norm1为标准化函数
y=c(:,px+1:px+py); %截取标准化的因变量
E0=c(:,1:px);
F0=c(:,px+1:px+py);
A=E0'*F0*F0'*E0;
w(:,1)=maxdet(A); %求最大特征向量
t(:,1)=E0*w(:,1); %提取主成分
E(:,1:px)=E0-t(:,1)*(E0'*t(:,1)/(t(:,1)'*t(:,1)))'; % 获得回归系数
p(:,1:px)=(E0'*t(:,1)/(t(:,1)'*t(:,1)))';
for i=0:px-2
B(:,px*i+1:px*i+px)=E(:,px*i+1:px*i+px)'*F0*F0'*E(:,px*i+1:px*i+px);
w(:,i+2)=maxdet(B(:,px*i+1:px*i+px));
% maxdet为求最大特征值的函数
t(:,i+2)=E(:,px*i+1:px*i+px)*w(:,i+2);
p(:,px*i+px+1:px*i+2*px)=(E(:,px*i+1:px*i+px)'*t(:,i+2)/(t(:,i+2)'*t(:,i+2)))';
E(:,px*i+px+1:px*i+2*px)=E(:,px*i+1:px*i+px)-t(:,i+2)*(E(:,px*i+1:px*i+px)'*t(:,i+2)/(t(:,i+2)'*t(:,i+2)))';
end
for s=1:px
q(:,s)=p(1,px*(s-1)+1:px*s)';
end
[n,d]=size(q);
for h=1:px
iw=eye(d);
for j=1:h-1
iw=iw*(eye(d)-w(:,j)*q(:,j)');
end

wh(:,h)=iw*w(:,h);
end
forj=1:py
zr(j,:)=(regress1(y(:,j),t))'; %求回归系数
endfor j=1:px
fori=1:py %生成标准化变量的方程的系数矩阵
w1=wh(:,1:j);
zr1=(zr(i,1:j))';
f0(i,:,j)=(w1*zr1)';
end

[normxy,meanxy,covxy]=norm1(C); %normxy标准化后的数据矩阵
%meanxy每一列的均值
%covxy每一列的方差
ccxx=ones(py,1)*meanxy(1,1:px);
ccy=(covxy(1,px+1:px+py))'*ones(1,px);
ccx=ones(py,1)*(covxy(1,1:px));
ff=ccy.*f0(:,:,j)./ccx;
fff=-(sum((ccy.*ccxx.*f0(:,:,j)./ccx)')-meanxy(1,px+1:px+py))';
FF(:,:,j)=[fff,ff]; %生成原始变量方程的常数项和系数矩阵
end
当程序运行至粉红色字段时，程序报错，试了很多方法，均不行。