Numerical analysis of 3D space-fractional neutral-type delayed reaction–diffusion equations using a high-order difference technique