Various numerical methods based on the theory of distribution of sources, of distribution of doublets and on the theory of potential of acceleration, were advanced for the study of unstationary supersonic movements around thin airfoils. The paper presents an adaptation of “Harmonic Gradient Method” of sources distribution theory.