Unpolarized QED parton distribution functions in NLO

Perturbative solutions for unpolarized QED parton distribution and fragmentation functions are presented explicitly in the next-to-leading logarithmic approximation. The scheme of iterative solution of QED evolution equations is described in detail. Terms up to (alpha 3L2) are calculated analytically, where L=ln(mu F2/me2) is the large logarithm which depends on the factorization energy scale mu F >> m e . The results are process independent and relevant for future high-precision experiments.