  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.
