We describe in detail two numerical simulation methods valid to study systems whose thermostatistics is described by generalized entropies, such as Tsallis. The methods are useful for applications to non-trivial interacting systems with a large number of degrees of freedom, and both short-range and long-range interactions. The first method is quite general and it is based on the numerical evaluation of the density of states with a given energy. The second method is more specific for Tsallis thermostatistics and it is based on a standard Monte Carlo Metropolis algorithm along with a numerical integration procedure. We show here that both methods are robust and efficient. We present results of the application of the methods to the one-dimensional Ising model both in a short-range case and in a long-range (non-extensive) case. We show that the thermodynamic potentials for different values of the system size N and different values of the non-extensivity parameter q can be described by scaling relations which are an extension of the ones holding for the Boltzmann-Gibbs statistics (q = 1). Finally, we discuss the differences in using standard or non-standard mean value definitions in the Tsallis thermostatistics formalism and present a microcanonical ensemble calculation approach of the averages.