Numerical evaluations for shock and discontinuity velocity related to similarity solutions of quasi-linear hyperbolic systems