Lipari and Stanev developed a method for range fluctuation of high energy muons, stressing the importance of accounting for the fluctuations of the energy loss in radiative processes in 1991 and, now, their method has become the basement for the energy determination of high energy muons through the measurement of the Cherenkov light yields due to those muons in KM3 physics. Once, Takahashi et al. developed a method for the investigation on the depth intensity relation of high energy muons in which all the stochastic processes concerned are taken into account exactly (1983). Now, we make the method by Takahashi et al. revival for the same purpose of the application to the analysis of future KM3 physics. In the present paper, our concern is restricted to the introduction to the fundamental of our method and some subsequent results thereby in which the real simulated behaviors of high energy muons from $10^{12}$eV to $10^{18}$eV, the survival probabilities of high energy and so on are included. The discussion around the practical application of our method to the KM3 physics is entrusted in the subsequent papers. As far as the survival probability of high energy muons is concerned, our method gives nearly the same results to Lipari and Stanev's in some regions and gives the deviated results from theirs in another ones. Thus, we examine the application limit of their method and clarify the reason why, comparing with our method. The most distinct difference in the both methods may become apparent in the treatment on the Cherenkov light yields spectrum by which one may estimate the energies of the muons concerned. We will mention to them in subsequent papers.

The validity of our Monte Carlo simulation procedure (the time sequential procedure) had been verified by the corresponding analytical procedure which is independent of our method methodologically. Also, the results obtained by our procedure are compared with those obtained by the different Monte Carlo simulation procedure (the v_{cut} procedure) which has been exclusively utilized by different authors and the agreement between ours and theirs are found to be good. By utilizing our Monte Carlo procedure, the validity of which is guaranteed in two different procedures, we investigate not only the fluctuation of high energy muons themselves but also the fluctuation of the various quantities related to the energy losses by the muons, which are difficult to obtain by the v_{cut} procedure. Namely, we obtain the fluctuation on energy losses of muons in the present paper(Part1). In a subsequent paper (Part2), based on the idea developed in Part1, we will examine the correlation between the primary energies of muons and their energy deposits and the correlation between the energy deposits of the muons concerned and their Cherenkov light yields in the KM3 detector in detail.

The validity of our Monte Carlo simulation procedure ({\it the integral method}) had been verified by the corresponding analytical procedure of which is quite independent of our method methodologically. Also, the results obtained by our procedure are compared with those obtained by the different Monte Carlo simulation procedure ({\it the differential method}) which have been exclusively utilized by the different authors and the agreement between them are found to be well. By utilizing our Monte Carlo procedures, the validity of which is guaranteed in two different procedures, we investigate not only the fluctuation of high energy muons themselves but also fluctuation of the various quantities related to the energy losses by the muons, which are difficult to obtain by {\it the differential method}. Namely, we obtain fluctuation on energy losses of the muons, fluctuation on Cherenkov lights due to the accompanied cascade showers initiated by the muon and the correlations between them . Finally, we obtain the transition curves for Cherenkov lights in KM3 detector, taking into account of all possible fluctuations in the stochastic processes and point out the difficulty of the reliable estimation of the energy of the muons which are resultants of muon neutrino events in the KM3 detectors.

We calculate high and (ultra-)high energy upward-going muon neutrino propagation through the Earth and the induced muon energy distribution near the one cubic kilometer detector using the Monte Carlo simulation, due to both charged current interaction and neutral one. The initiated neutrino energies on the surface of the Earth are 1PeV, 1EeV and 1ZeV. The mean free paths of (ultra-)high energy neutrino events generated by the deep inelastic scattering may be comparable with the diameter of the Earth or less than it. Therefore, the induced muon production distribution is influenced by the change of the densities interior to the Earth. Furthermore, in such situation, the contribution from the neutral current neutrino interaction to the induced muon production distribution cannot be neglected. We report several examples of the deep inelastic scattered depth of muon neutrino in the Earth and the induced muon energy distribution near the detector.

By using the integral method in the muon propagation through water, we calculate the range fluctuation of high and ultra high energy muons. Many authors divide all radiative processes into two parts, namely, the continuous part and radiative part in their Monte Carlo simulation in order to consider the fluctuation in the both ranges and energies of the muons, while we treat all stochastic processes as exactly as possible, without the introduction of the continuous parts in all stochastic processes. The validity of our Monte Carlo method is checked by the corresponding analytical method which is methodologically independent on the Monte Carlo procedure. Accompanied cascade showers are generated by the direct electron pair production, bremsstrahlung and photo-nuclear interaction. These showers are calculated by the exact Monte Carlo Method in one dimensional way. We report survival probabilities, range distributions and examples of individual muon behavior.