loops over a 3d array

[Marco Caliari]

Dear users,

I'm wondering if it is possible to vectorize the following code

Hx = rand(n); % n natural number
Hy = rand(n);
Hz = rand(n);
u = rand(n,n,n);
for k = 1:n
  us(:,:,k) = zeros(n);
  for j = 1:n
    us(:,:,k) = us(:,:,k)+Hz(k,j)*(Hy*u(:,:,j)*Hx');
  end
end

Thanks,

Marco