simultaneous diagonalization of two matrices

How is simultaneous diagonalization performed on two matrices ( say A and B ) using matlab ?

Answer 1

[Matt J]

If A and B are the same size

out = bsxfun(@times, cat(3,A,B), eye(size(A)) );