Introduction ============ **pyEQUIB** package is a collection of `Python `_ programs developed to perform plasma diagnostics and abundance analysis using emission line fluxes measured in ionzed nebulae. It uses the `AtomNeb Python Package `_ to read collision strengths and transition probabilities for collisionally excited lines (CEL), and recombination coefficients for recombination lines (RL). This Python package can be used to determine interstellar extinctions, electron temperatures, electron densities, and ionic abundances from the measured fluxes of emission lines. It mainly contains the follwing API functions written purely in Python: Collisional Excitation Unit --------------------------- **API functions for collisionally excited lines (CEL)** have been developed based on the algorithm of the FORTRAN program `EQUIB `_ written in FORTRAN by `Howarth & Adams (1981) `_. The program EQUIB calculates atomic level populations and line emissivities in statistical equilibrium in multi-level atoms for different physical conditions of the stratification layers where the chemical elements are ionized. Using the Python implementation of the program `EQUIB `_, electron temperatures, electron densities, and ionic abundances are determined from the measured fluxes of collisionally excited lines. Recombination Unit ------------------ **API functions for recombination lines (RL)** have been developed based on the algorithm of the recombination scripts by X. W. Liu and Y. Zhang included in the FORTRAN program `MOCASSIN `_. These API functiosn are used to determine ionic abundances from recombination lines for some heavy element ions. Reddening Unit -------------- **API functions for reddening and extinctions** have been developed according to the methods of the reddening law functions from `STSDAS IRAF Package `_, which are used to obtain interstellar extinctions and deredden measured fluxes based on different reddening laws.