A method to evaluate nonlinear centroid correlation functions is presented that is amenable to simple numerical computation. It can be implemented with the centroid molecular dynamics method for approximate quantum dynamics with no additional assumptions. Two nonlinear correlation functions are evaluated for a model potential using this scheme and compared with results from exact quantum calculations.