This talk summarizes recent efforts and results in the simulation of supersymmetric gauge theories using numerical lattice simulations. I will show how non-perturbative properties like the low energy bound state spectrum and the phase transitions can be determined from simulations of N=1 supersymmetric Yang-Mills theory. These show some distinguished properties of this theory. I will also discuss the approach for the simulation of N=4 supersymmetric Yang-Mills theory and possible future developments.