Application of the Marquardt least-squares method (MLSM) to the estimation of non-linear parameters of functionsused for representing various lightning current waveshapes is presented in this paper. Parameters are determined for the Pulse,Heidler’s and DEXP function representing the first positive, first and subsequent negative stroke currents as given in IEC62305-1 Standard Ed.2, and also for some other fast- and slow-decaying lightning current waveshapes. The results prove theability of the MLSM to be used for the estimation of parameters of the functions important in lightning discharge modeling.