A numerical approach to evaluate the describing function (DF)
A numerical approach to evaluate the DF based on the Gaussian quadrature of Chebyshev is given. This method has the advantage of being applicable to functions in conditions that are less restrictive than the polynomial expansion or Fourier series expansion.