Assuming Gaussian initial conditions, we show that in a certain, angle-averaged squeezed limit the N-point function of matter is related to the response of the matter power spectrum to a long-wavelength density perturbation, P− 1 d nP (k| δ L)/dδ L n| δ L= 0, with n= N− 2. By performing N-body simulations with a homogeneous overdensity superimposed on a flat Friedmann-Robertson-Lemaȋtre-Walker (FRLW) universe using the separate universe approach, we obtain measurements of the nonlinear matter power spectrum response up to n= 3, which is equivalent to measuring the fully nonlinear matter 3-to 5-point function in this squeezed limit. The sub-percent to few percent accuracy of those measurements is unprecedented. We then test the hypothesis that nonlinear N-point functions at a given time are a function of the linear power spectrum at that time, which is predicted in an Einstein-de Sitter (EdS) universe by standard perturbation theory (SPT) and its variants that are based on the ideal pressureless fluid equations. Specifically, we compare the responses computed from the separate universe simulations and simulations with a rescaled initial (linear) power spectrum amplitude. We find discrepancies of 10% at k≃ 0.2–0.5 h Mpc− 1 for 5-to 3-point functions at z= 0. The discrepancy occurs at higher wavenumbers at z= 2. Thus, theoretical predictions that are insensitive to the growth history, such as SPT and its variants assuming EdS, even when carried out to arbitrarily high order, are guaranteed to fail to describe matter N-point functions (N> 2) around that scale.