The aim of this paper is to demonstrate, by simple numerical simulations, the main transport properties of disordered electron systems. These systems undergo the metal insulator transition when either Fermi energy crosses the mobility edge or the strength of the disorder increases over critical value. We study how disorder affects the energy spectrum and spatial distribution of electronic eigenstates in the diffusive and insulating regime, as well as in the critical region of the metal-insulator transition. Then, we introduce the transfer matrix and conductance, and we discuss how the quantum character of the electron propagation influences the transport properties of disordered samples. In the weakly disordered systems, the weak localization and anti-localization as well as the universal conductance fluctuation are numerically simulated and discussed. The localization in the one dimensional system is described and interpreted as a purely quantum effect. Statistical properties of the conductance in the critical and localized regimes are demonstrated. Special attention is given to the numerical study of the transport properties of the critical regime and to the numerical verification of the single parameter scaling theory of localization. Numerical data for the critical exponent in the orthogonal models in dimension 2 < d ≤ 5 are compared with theoretical predictions. We argue that the discrepancy between the theory and numerical data is due to the absence of the self-averaging of transmission quantities. This complicates the analytical analysis of the disordered systems. Finally, theoretical methods of description of weakly disordered systems are explained and their possible generalization to the localized regime is discussed. Since we concentrate on the one-electron propagation at zero temperature, no effects of electronelectron interaction and incoherent scattering are discussed in the paper.