Linear Systems of Equations

xxxxxxxxxx

 

1

var('x1,x2,x3,x4')

2

x=vector(SR,[x1,x2,x3,x4])

3

A=matrix(SR,[[8,-5,-2,-2],[-2,4,-4,1],

4

             [-2,1,4,4],[4,1,3,10]])

5

b=vector(SR,[35,-8,-11,10])

6

Ab=A.augment(b,subdivide=True)

7

sAb=solve(list(A*x-b),x1,x2,x3,x4)

8

ppl=[['Decision Methods:'],

9

     [Ab.rref()],[sAb],

10

     [A.solve_right(b)],[A.inverse()*b],

11

     ['For the Linear System of Equations: '],[Ab]]

12

table(ppl)

Language: Sage Gap GP HTML Macaulay2 Maxima Octave Python R Singular

---

Messages