import math import numpy rot = numpy.empty( (2,2) ) angle = 30. * math.pi/180. rot[0][0] = math.cos(angle) rot[0][1] = math.sin(angle) rot[1][0] = -math.sin(angle) rot[1][1] = math.cos(angle) print rot print " " a = numpy.zeros( (18,2) ) angle = 0. incr = 2*math.pi/ 18. for i in range(18): a[i][0] = math.cos(angle) a[i][1] = math.sin(angle) angle += incr print a print " " print numpy.round(a, 5) print " " b = numpy.dot(a,rot) # multiply the points array by the rotation matrix b = 2. * b # scale the points by 2. print b print " " #print numpy.round(b, 5) print b.round(3)