Optimizing code: various distance

Hi all,

Is there a way to improve the performance of the following code?

D=zeros(size(X,1),size(Y,1));
for l=1:size(X,1)
    D(l,:)=sum((abs(bsxfun(@minus,Y,X(l,:))).^p),2).^(1/p);
end

Thank you in advance