  Using repeated Laplace transform, We find an analytical solution for DGLAP
evolution equations for extracting the pion, kaon and proton Fragmentation
Functions (FFs) at NLO approximation. We also study the symmetry breaking of
the sea quarks Fragmentation Functions, $D_{\bar q}^h (z,Q^2)$ and simply
separated them according to their mass ratio. Finally, we calculate the total
Fragmentation Functions of these hadrons and compare them with experimental
data and those from global fits. Our results show a good agreement with the FFs
obtained from global parameterizations as well as with the experimental data.
