Untitled

1861 days ago by matrix

A=matrix(4,6,[1,1,1,0,6,0,0,0,2,1,7,2,0,-1,0,2,1,4,0,-1,2,2,7,1]) B=A.rref() c=vector([2,1,3,1]) d=vector([-9/2,6,-3/2,5]) print B print print 3*c+8*d 
       
[   1    0    0    0    2 -9/2]
[   0    1    0    0    1    6]
[   0    0    1    0    3 -3/2]
[   0    0    0    1    1    5]
(-30, 51, -3, 43)
[   1    0    0    0    2 -9/2]
[   0    1    0    0    1    6]
[   0    0    1    0    3 -3/2]
[   0    0    0    1    1    5]
(-30, 51, -3, 43)