Numerical study on ergodic properties of triangular billiards