def METHOD_NAME(self, a0, a1, vertices):
    l = rangen.uniform(self.multMin, self.multMax)
    theta = rangen.uniform(0.0, 2.0*pi)
    phi = rangen.uniform(0.0, pi)
    T = l*rotationMatrix(Vector(cos(theta)*sin(phi),
                                sin(theta)*sin(phi),
                                cos(phi)))
    verts = vector_of_Vector()
    for x in vertices:
        verts.append(T*x)
    return T*a0, T*a1, Polyhedron(verts), T