Masaharu Nagayama (one of investigators of the present research project) has devised a computer code that can analyze bifurcation structures in a neighborhood of double bifurcation points. This code deals with bifurcation phenomena of pulse solutions to mono-stable reaction-diffusion systems, and is equipped with the following two functions : (1)It can find a critical point and construct its bifurcation branch, (2)It can extend existing bifurcation branches. In order to devise the code, we consider the reaction-diffusion system on a finite interval (0,L) subject to the periodic boundary condition where L is a fixed large positive number. From the phase condition we obtain the equation that determines the propagating velocity of traveling pulse, and by the Keller method we express the dependency on a parameter p included in the equation systems. The problem formularized as in the above is numerically solved by the Newton method in the computer code. We note that a solution is a set of {solutions to reaction-diffusion systems, c, p}. When a traveling pulse bifurcates from a standing pulse, there appear two zero-eigenvalues, one of which is a trivial one trivial one corresponding to parallel translation. Our code applies to not only this case but also the cases where two crucial zero-eigenvalues exist except the trivial one.The head investigator have dealt with standing and traveling combustion pulses of a mathematical model for self-propagating high-temperature syntheses including both the cooling effect and raw material supply system. Employing a piece-wise constant function for the reaction term, we have studied the existence of pulse solutions in a mathematically rigorous way, and also the collision dynamics of combustion pulses on a circle domain.