||'''Book'''|| Linear Algebra ||
 ||'''Author'''|| Gilbert Strang ||
 ||'''Edition'''||  ||

/* code_begins */

{{{
print "Given equations are :\nx+3y+3z=1\n2x+6y+9z=5\n-x-3y+3z=5"
print "Converting these into matrix form we obtain"
A=matrix([[1,3,3],[2,6,9],[-1,-3,3]])
print "Coefficient Matrix A=\n",A
var('x,y,z')
X=matrix([[x],[y],[z]])
print "Matrix X=\n",X
b=matrix([[1],[5],[5]])
print "Matrix b=\n",b
print "Solving given equations using solve_right()function"
K=A.solve_right(b)
print "Obtained values of x,y and z are:"
x=K[0]
y=K[1]
z=K[2]
print "x=",x
print "y=",y
print "z=",z
C=matrix([[1,3,1,2],[2,6,4,8],[0,0,2,4]])
print "Given matrix C=\n",C
e=matrix([[1],[3],[1]])
print "Given matrix e=\n",e
var('p,q,r,s')
D=matrix([[p],[q],[r],[s]])
print "Given matrix D=\n",D
print "Solving these using the function solve_right()"
print "Let M be the resultant matrix"
M=C.solve_right(e)
print "Matrix M=\n",M
print "Assigning these values to their relative values in D matrix"
p=M[0]
q=M[1]
r=M[2]
s=M[3]
print "Obtained values are:"
print "p=",p
print "q=",q
print "r=",r
print "s=",s


Result:

   	

Given equations are :
x+3y+3z=1
2x+6y+9z=5
-x-3y+3z=5
Converting these into matrix form we obtain
Coefficient Matrix A=
[ 1  3  3]
[ 2  6  9]
[-1 -3  3]
Matrix X=
[x]
[y]
[z]
Matrix b=
[1]
[5]
[5]
Solving given equations using solve_right()function
Obtained values of x,y and z are:
x= (-2)
y= (0)
z= (1)
Given matrix C=
[1 3 1 2]
[2 6 4 8]
[0 0 2 4]
Given matrix e=
[1]
[3]
[1]
Given matrix D=
[p]
[q]
[r]
[s]
Solving these using the function solve_right()
Let M be the resultant matrix
Matrix M=
[1/2]
[  0]
[1/2]
[  0]
Assigning these values to their relative values in D matrix
Obtained values are:
p= (1/2)
q= (0)
r= (1/2)
s= (0)


}}}

/* code_ends */

 * '''Solution by''': 
   * <T.R.N.Karthik>, <Student>, <Gitam University>