## page was renamed from Kreyszig-18.9 ex1
 ||'''Book'''|| Advanced Engineering Mathematics ||
 ||'''Author'''|| Erwin Kreyszig ||
 ||'''Edition'''|| 8th Edition ||

/* code_begins */

{{{
A=A0=matrix([[6,4,1,1],[4,6,1,1],[1,1,5,2],[1,1,2,5]])
print A
print A0
v11=0 
a21= A0[1,0]
print A
print a21
a31=A0[2,1]
a41=A0[3,1]
print a31,a41
x1= a21*a21 + a31*a31 + a41*a41
print x1
s1= x1^0.5
print s1
x2=(1+a21/s1)/2
print x2
v21= x2^0.5
print v21
v31=(a31*a21)/(2*v21*s1)
print v31
v41=(a41*a21)/(2*v21*s1)
print v41
v1=matrix([[0],[v21],[v31],[v41]])
print v1
I= identity_matrix(4)
print I
p1= I-(2*(v1)*(v1.transpose()))
print p1
A1= p1*A0*p1
print A1
}}}

/* code_ends */

 * '''Solution by''': 
   * V.J.Ramchand,Student, GITAM University
   * G.Vamsi krishna, Student, GITAM University