`dnorm.dual` will make straightfoward a computation (in R), that
works both with numeric or dual objects. `dnorm` will call `dnorm.dual`
if any of the objects is of class dual, or `stats::dnorm` is all objects
are of class numeric. As `stats::dnorm` is in written in C it is factor.
If you care for performance, use `stats::dnorm` directly for non dual numbers,
and `dnorm.dual` for dual numbers.