Calculation of neutron response functions in complex geometries with the MCNP code