Gradient density estimation in arbitrary finite dimensions using the method of stationary phase