Universal scaling behavior of directed percolation around the upper critical dimension

In this work we consider the steady state scaling behavior of directed percolation around the upper critical dimension. In particular we determine numerically the order parameter, its fluctuations as well as the susceptibility as a function of the control parameter and the conjugated field. Additionally to the universal scaling functions, several universal amplitude combinations are considered. We compare our results with those of a renormalization group approach.