ABSTRACT: A general method of simulation of processes in dusty based on special programs is presented here. It is pos-sible to prepare the modeling of the dusty in volcano like the dust sound waveguides. Dusty is in state of the plasma .Waveguides are formed by the distribution of dusty particles with various masses m = m(x) in trans-verse coordinate. The dust sound waves propagate along the longitudinal z-direction. In the case of contact of dusty plasma with a semi-infinite dielectric, there exists the dust acoustic mode that possesses the negative group velocity (backward wave) in the specified interval of wave numbers. For analysis it is necessary to use the special numerical methods of calculation of the equations with boundary conditions. Simulation of ion sound wave propagation shows a new dispersion between frequency and wave vector. In some region of pa-rameters of dusty the negative dispersion of wave takes place. This means that the phase and group velocities of wave are opposite (negative dispersion). This phenomenon takes place, when the mass of dust particles has the maximum in the center of the waveguide. The negative dispersion caused the instability in dusty, which open the possibility to create a new phenomenon in dusty including the high temperature and the flame.