print 'finding inverse of A1'
print '---------------------'
A1=matrix([[0,2],[3,0]])
x1=(A1[0,0]*A1[1,1])-(A1[0,1]*A1[1,0])
print '---A1---'
print A1
print '---det of A1---'
print x1
x2= matrix ([[A1[1,1],-A1[0,1]],[-A1[1,0],A1[0,0]]])
print '---inverse---'
print (1/x1)*x2
print 'finding inverse of A2'
print '---------------------'
A2=matrix([[2,0],[4,2]])
x3=(A2[0,0]*A2[1,1])-(A2[0,1]*A2[1,0])
print '---A2---'
print A2
print '---det of A2---'
print x3
x4= matrix ([[A2[1,1],-A2[0,1]],[-A2[1,0],A2[0,0]]])
print '---inverse---'
print (1/x3)*x4
print 'inverse of A3'
print '--------------'
var('theta')
A3=matrix([[cos(theta),-sin(theta)],[sin(theta),cos(theta)]])
x5=(A3[0,0]*A3[1,1])-(A3[0,1]*A3[1,0])
print '---A3---'
print A3
print '---det of A3---'
print x5
x6= matrix ([[A3[1,1],-A3[0,1]],[-A3[1,0],A3[0,0]]])
print '---inverse---'
print (1/x5)*x6

Strang-1.6-1-U (last edited 2010-12-18 06:53:07 by jejiramchand)