Kernel smoothing for multivariate observations

Hello,

May I ask if there is a R-code to generate the kernel density estimates in fields?  I need to estimate the surface with "epanechnikov" kernel for multidimensional data, but all I can find is the intensity estimation for point pattern which is different. Your answer are highly appreciated!

Thanks!

Sincerely yours,

Juan