Hello,
I am trying to solve a system of matrix equations of the form:
transpose(B)*(Mi)*B==Ni
where i is in range(p) for some fixed p > 9, Mi and Ni have integer
entries, and all matrices are of size p by p. I want to solve for B
(it need not have integer entries).
I have tried a brute force approach. I can make a list of equations
corresponding to each entry, and in test cases (p=2) it will solve for
B. In general, I would like to define a new command that inputs p and
uses the "solve( , )" function. Something like the following
def matrix_solve( , ,p)
solutions=solve(list of equations from entries, x1,x2, ... xp)
return solutions
The first input is a list of equations which I can create and input
into solve. The trouble is that the second input must be of the form
x1,x2, ... xp and cannot be a list, string, or tuple.
Question 1: Is there a clever way to input x1,x2, ... xp into the
solve function given some p?
Question 2: Is there a more clever/efficient way to solve systems of
matrix equations?
Thank you in advance for any help,
M. Miller
--~--~---------~--~----~------------~-------~--~----~
To post to this group, send email to sage-support@googlegroups.com
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://sage.math.washington.edu/sage/ and http://sage.scipy.org/sage/
-~----------~----~----~----~------~----~------~--~---
