The fragmentation functions for or quark to or meson are derived up to QCD next-to-leading order (NLO). They are further computed numerically and presented in figures precisely. In order to reach a higher accuracy, we have also tried properly to apply them to estimating the and production at a -factory (an collider running at the energy of the -boson pole).