In this paper, by using our improved plane wave decomposition method, we study the scars in the eigenfunctions of the stadium billiard from a very low state to as high as about the 1 millionth state. In the systematic searching for scars of various types, we have used the approximate criterion based on the quantization of the classical action along the unstable periodic orbit supporting the scar. We have analysed the profile of the integrated probability density along the orbit. We found that the maximal integrated intensity of different types of scars scales in different way with the , which confirms qualitatively and quantitatively the existing theories of scars such as that of Bogomolny and Robnik.