In this paper the approximating capabilities of fuzzy systems with overlapping Gaussian concepts are considered. A new method for the computation of the system coefficients is provided, showing that it guarantees uniform approximation of the derivatives of the target function. Moreover, the connection with radial basis functions approximations is highlighted.