Hi again,
On Dec 17, 11:05 am, Simon King <[email protected]> wrote:
--snip--
> Namely, one can do the following:
> sage: n=2
> sage: A=Matrix([[var('a_%d_%d'%(i,j)) for j in range(n)] for i in range
> (n)])
> sage: solve((A*A.transpose()).list(), A.list())
Oops, I forgot that I wanted orthogonal matrices. But Helas:
sage: solve((A*A.transpose()-MatrixSpace(QQ,2)(1)).list(), A.list())
---------------------------------------------------------------------------
ValueError Traceback (most recent call
last)
...
ValueError: Unable to solve [a_0_1^2 + a_0_0^2 - 1, a_0_1*a_1_1 +
a_0_0*a_1_0, a_0_1*a_1_1 + a_0_0*a_1_0, a_1_1^2 + a_1_0^2 - 1] for
([a_0_0, a_0_1, a_1_0, a_1_1],)
Bad luck indeed.
Best regards,
Simon
--~--~---------~--~----~------------~-------~--~----~
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to
[email protected]
For more options, visit this group at
http://groups.google.com/group/sage-support
URLs: http://www.sagemath.org
-~----------~----~----~----~------~----~------~--~---
