> I am to trying to learn Maxima and want to solve a transformation
> matrix from the following eq:
>
> Rdx = T.Rrl
>
> Rdx and Rrl are defined like :
>
> Rdx=matrix([sh*sp*sr
> +ch*cr,cp*sr,ch*sp*sr-cr*sh],[cr*sh*sp-ch*sr,cp*cr,sh*sr
> +ch*cr*sp],[cp*sh,-sp,ch*cp])
>
Write Rdx:matrix(...) instead of Rdx=matrix(...)
then inverse(Rdx) reports no problem:
/************* begin Maxima session ****************/
(%i11) Rdx:matrix([sh*sp*sr
+ch*cr,cp*sr,ch*sp*sr-cr*sh],[cr*sh*sp-ch*sr,cp*cr,sh*sr
+ch*cr*sp],[cp*sh,-sp,ch*cp]);
[ sh sp sr + ch cr cp sr ch sp sr - cr sh ]
[ ]
(%o11) [ cr sh sp - ch sr cp cr sh sr + ch cr sp ]
[ ]
[ cp sh - sp ch cp ]
(%i12) invert(Rdx).Rdx, factor;
[ 1 0 0 ]
[ ]
(%o12) [ 0 1 0 ]
[ ]
[ 0 0 1 ]
/************* end Maxima session ****************/
--
Mario Rodriguez Riotorto
www.biomates.net