We perform Monte Carlo simulations of a supersymmetric matrix model, which
is obtained by dimensional reduction of 4D SU(N) super Yang-Mills theory. T
he model can be considered as a four-dimensional counterpart of the IIB mat
rix model. We extract the space-time structure represented by the eigenvalu
es of bosonic matrices. In particular we compare the large-N behavior of th
e space-time extent with the result obtained from a low-energy effective th
eory. We measure various Wilson loop correlators which represent string amp
litudes and we observe a non-trivial universal scaling in N. We also observ
e that the Eguchi-Kawai equivalence to ordinary gauge theory does hold at l
east within a finite range of scale. Comparison with the results for the bo
sonic case clarifies the role of supersymmetry in the large-N dynamics. It
does affect the multi-point correlators qualitatively, but the Eguchi-Kawai
equivalence is observed even in the bosonic case.