Matlab - Multiplicar N matrices

%% http://www.lawebdelprogramador.com/foros/Matlab/1316740-Multiplicar_N_matrices.html

% Hola soy nuevo en este foro y llevo ya unos dos o tres añitos con matlab 

% pero a bajo nivel.

% 

% He estado haciendo un programa para sacar la factorizacion LU de una

% matriz, aunque se que el propio programa tiene una funcion para ello,

% mi duda es:

% 

% Con la funcion eval (interna del programa) he conseguido crear N 

% matrices, que he llamado M, a las cuales les he dado valores dentro

% de un bucle; mi duda es, ahora necesito multiplicar todas esas variables

% (Matrices) M, que reciben los nombres (M1,M2,M3,...,Mk,...,MN) y para

%  ello necesito llamarlas a todas, entonces, ¿existe una funcion que

% las llame sucesivamente o como programo en for, while o if para que me

% multiplique esas matrices???; porque no se como hacer para

% llamar a una matriz generica Mk para hacer algo con ella...

% Si no me he explicado bien pregunten que me aclaro mejor.

clc

clear all

i=1;

n_matrices=3;

%generando  matrices aleatorias

while i<=n_matrices

    R=round(rand(3,3))+1;

    v=genvarname(['M',num2str(i)]);

    eval([v, '=R']);

    i=i+1;

end

%multiplicando matrices generadas aleatoriamente

a=whos('M*');

Producto=eval(a(1).name);

for i=2:n_matrices

Producto=Producto*eval(a(i).name);

end

%producto de las matrices generadas aleatoriamente

display(Producto)

EJECUCION

>> genvarname3

M1 =

     2     2     1

     1     1     2

     2     2     1

M2 =

     1     2     1

     2     1     2

     1     1     1

M3 =

     1     1     1

     2     2     1

     1     2     1

Producto =

    28    35    21

    20    25    15

    28    35    21

>> genvarname3

M1 =

     1     2     2

     1     2     1

     2     1     1

M2 =

     1     1     1

     1     1     1

     1     1     2

M3 =

     2     1     2

     2     1     2

     1     2     1

Producto =

    27    24    27

    21    18    21

    21    18    21

clc,clear all

for k = 1:5

v = genvarname(['A',num2str(k)]);

eval([v ' = k+3'])

end

who