This paper provides FORTRAN subroutines for the calculation of the partially-Correlated quadratic-Speed-Dependent Hard-Collision (pCqSDHC) profile and of its two limits: the quadratic-Speed-Dependent Voigt (qSDV) and the quadratic-Speed-Dependent Hard-Collision (qSDHC) profiles. Numerical tests successfully confirm the analytically derived fact that all these profiles can be expressed as combinations of complex Voigt probability functions. Based on a slightly improved version of the CPF subroutine Humlicek. J Quant Spectrosc Radiat Transfer 1979;21:309 for the calculation of the complex probability function, we show that the pCqSDHC, qSDHC and qSDV profiles can be quickly calculated with an accuracy better than 10−4. •FORTRAN subroutines of some speed-dependent profiles are provided.•These profiles can be expressed as combinations of the complex probability function.•They can be accurately calculated with reasonable computer cost.