Computes breakthrough curve for continuous (step) contaminant source in 1D with advection, dispersion, and retardation.