Numerical algorithms for eveluation of higher order moments for semi-Markov rewards processes are presented. Results of numerical experiments are given and commented.