[Maxima-lang-es] Duda sobre integral sobre matriz

[Nicolas G. <nic...@gm...>]

Hola,

Estaba calculando unas matrices de rigidez (Para verificacion del
método de elementos finitos)  y obtuve un error, que creo se debe a
memoria.

Primero defino unas matrices que son.

(%i1) N4:[1/4*(1+r)*(1+s),
   1/4*(1-r)*(1+s),
   1/4*(1-r)*(1-s),
   1/4*(1+r)*(1-s)];

(%i2) H4: genmatrix(lambda([i,j], 0), 2, 8)$
for i:1 thru 4 do(
    H4[1,2*i-1]:N4[i],
    H4[2,2*i]:N4[i]
);
H4;


Luego quiero realizar la integral de la siguiente matriz y no funciona:

(%i3) M4:rho*integrate(integrate(transpose(H4).H4,r,-1,1),s,-1,1)$


context: too many contexts.
 -- an error. To debug this try: debugmode(true);


Si remplazo la integral de la matriz por un par de ciclos en el que
integro elemento a elemento la matriz
no tengo ese problema:

(%i4) M4:transpose(H4).H4;

(%i5) for i:1 thru 8 do(
  for j:1 thru 8 do(
      M4[i][j]:rho*integrate(integrate(M4[i][j],s,-1,1),r,-1,1)
  )
);

Antes había hecho la integral de la matriz sin problema en la versión
5.20 de Maxima pero leugo de pasar
a la 5.24 o 5.25 no me funciona.

Esto lo probé en Windows 7.


¿Saben a qué se debe este error?

Buena suerte y gracias,

-- 
Nicolás Guarín Zapata
Ingeniero Físico
Especialista en Mecánica Computacional
Grupo de Mecánica Aplicada
http://mecanica.eafit.edu.co/
Universidad EAFIT
(574) 2619500 ext. 9676

Linux user #514760