Using a Mathematica notebook, we plot the Hill function for the exponent n ranging from 1 to 10:
For those wishing to use delay differential equation in Mathematica, there is a small package developped and available here. It is require for example to model the equations in Number of active transcription factor binding sites is essential for the Hes7 oscillator by S. Zeiser et al.