I'm trying to use LDA for fault diagnosis but I'm get problems to runs it

code

load d11_te.dat

d11 = zscore(d11); size 480*50

d11_te = zscore(d11_te) size 960*50

num1=size(d11,1);

num2=size(d11_te,1);

l=num1+num2;

X=[d11; d11_te];

Y=[ones(n1,1);2*ones(n2,1)];

% Y=Y';

%%%%%%%%%%%%%%%%%%%%%% Computing LFDA solution

[T,Z]=LFDA(X,Y,1);

function

for c=unique(Y')

Xc=X(:,Y==c);

nc=size(Xc,2);

error

Index exceeds matrix dimensions.

Error in LFDA (line 48)

Xc=X(:,Y==c);

Error in demo_LFDA (line 20)

[T,Z]=LFDA(X,Y,1);