Matlab - operaciones con arreglos de 3 dimensiones

Hola que tal. Quisiera saber cual es la forma más eficiente de realizar las siguientes operaciones

A=rand(2,3,15);
B=rand(3,5,15);
C=rand(3,5);
D=rand(15):

%opepacion1

for i=1:15
E(:,:,i)=A(:,:,i)*B(:,:,i);
end

%operacion2

for i=1:15
F(:,:,i)=A(:,:,i)*C;
end

%operacion3

for i=1:15
G(:,:,i)=A(:,:,i)*D(i);
end