Numerical study of scars in a chaotic billiard

We study numerically the scaling properties of scars in stadium billiard. Using the semiclassical criterion, we have searched systematically the scars of the same type through a very wide range, from ground state to as high as the 1 millionth state. We have analyzed the integrated probability density along the periodic orbit. The numerical results confirm that the average intensity of certain types of scars is independent of $\hbar$ rather than scales with $\sqrt{\hbar}$. Our findings confirm the theoretical predictions of Robnik (1989).