function B = partial_least_squares(X,Y,l)     

[~,~,~,~,B] = plsregress(X',Y',l);
B = B(2:end,:)';

end

 