  We extend the modular invariance property of the trace functions of vertex
operator algebra on the set of irreducible modules (Zhu's theory) to the case
of trace functions of intertwining operators.
