var('p,q,r')
r=matrix([[p,q],[q,r]])
print r
s=r*r;s
det(s)
print "r is positive definite for det not equal to zero"

<ananthalakshmi>, <student>, <prakasam>