In this paper we develop a model of an integrated production and distribution system comprised of multisite batch plants and several warehouse. An integration of production and distribution system may generate significant savings by trading off the costs associated with the whole, rather than minimizing production and distribution costs separately. The mathematical approach is presented in this research. The problem consists of two or more multiproduct batch plants producing several products to supply several warehouse. The demands of each product at warehouses are treated as uncertain parameters described by normal distribution functions. A mixed integer nonlinear programming(MINLP) formulation is proposed. The effectiveness of the application of the suggested approach for the production and distribution of multisite batch plants will be illustrated with an example problem. The problem is solved with GAMS/DICOPT++.