Effective simulations of semiconductor superlattices are presented in the paper. The simulations have been based on the Wannier function method approach where a new algorithm, inspired by Büttiker probes, has been incorporated into determining the Green function procedure. The program is of a modular structure, and its modules can either work independently, or interact with each other following a predefined algorithm. Such structuring not only accelerates simulations and makes the transport parameters possible to initially assess, but also enables accurate analysis of quantum phenomena occurring in semiconductor superlattices. In this paper, the capabilities of type I superlattice simulator, developed earlier, are presented, with particular emphasis on the new block where the Fermi levels are determined by applying Büttiker probes. The algorithms and methods used in the program are briefly described in the further chapters of our work, where we also provide graphics illustrating the results obtained for the simulated structures known from the literature.