In this paper, we investigate the dynamics of a mathematical model of obesity population within fluctuating social environment. A stochastic differential equation model is constructed by perturbing two social related parameters of the deterministic model with white noise terms characterized by Gaussian distribution having zero mean and unit spectral density. In order to compute the numerical solution of the stochastic models Euler-Maruyama numerical method is used. Confidence intervals for the overweight and obesity childhood population are computed using Monte Carlo method. Analysis of the numerical results reveals that small perturbations on the parameters are not a major driving force for dynamical transitions from the underlying deterministic model. In addition, numerical results indicate a close relationship between the amplitude of the fluctuation of the social environment parameters and the variability of forecasts for the incidence of the obesity in the population.