Simbad code is used to study the combined effect of external non linearities and space charge non linearities on the beam stability using a simple FODO lattice. Gaussian and parabolic particle distribution are used for these simulations and results are compared with Möhl and Metral theoretical results.