We have developed a code to simulate equilibrium properties in strongly correlated superconducting systems in real space. In the code the tJ model is solved by applying a mean-field approximation together with the Gutzwiller approximation. The model includes spin degrees of freedom and can accordingly be used to describe amongst other things the interplay of magnetism and superconductivity at the edges of superconductors.