Thread Subject: matrix multiplication in bsxfun way

I would like to do something like:

a=rand(2,2);
b=rand(2,2,3);
c=bsxfun(@times,a,b)

except I want to do MATRIX MULTIPLICATION instead of @times (which does POINT-WISE MULTIPLICATION).

the result I want should be a 2 x 2 x 3 array where
c(:,:,1) = a*b(:,:,1)
c(:,:,2) = a*b(:,:,2)
c(:,:,3) = a*b(:,:,3)

I vaguely recall ARRAYFUN being able to do something similar, but upon looking at the documentation for ARRAYFUN i couldn't come up with anything.

Thank you,
Juliette

Pure Matlab solution:
http://www.mathworks.com/matlabcentral/fileexchange/24260-multiple-same-size-linear-solver

MEX solution (highly recommended):
http://www.mathworks.com/matlabcentral/fileexchange/25977-mtimesx-fast-matrix-multiply-with-multi-dimensional-support

Bruno

http://www.mathworks.com/matlabcentral/fileexchange/25977-mtimesx-fast-matrix-multiply-with-multi-dimensional-support

ABSOLUTELY BRILLIANT!

THANK YOU BRUNO AND MATT , AND MOST OF ALL JAMES TURSA.