| Author |
Comment/Response |
Peter Pein
|
 04/25/12 11:03am
Hi,
are you looking for sth. like PseudoInverse[]?
In[1]:= coeff=RandomInteger[10,{2,3}]
Out[1]= {{9,1,6},{3,4,4}}
In[2]:= eqsys=And@@Thread[coeff.{x,y,z}=={a,b}]
Out[2]= 9x + y + 6z == a && 3x + 4y + 4z == b
In[3]:= psinv=PseudoInverse[coeff]
Out[3]= {{204/1813,-(141/1813)},{-(179/1813),417/1813},{26/1813,142/1813}}
In[4]:= soln=psinv.{a,b}
Out[4]= {(204 a)/1813-(141 b)/1813,-((179 a)/1813)+(417 b)/1813,(26 a)/1813+(142 b)/1813}
In[5]:= eqsys/.Thread[{x,y,z}->soln]//Expand
Out[5]= True
hth,
Peter
URL: , |
|
