Abstract
Mathematical models of earthquake disasters and responses must be useful for disaster prevention and its systems design. As a fundamental mathematical model this paper presents a 3-D simultaneous 1st-order differential equation which corresponds to a network with three kinds of sound, damaged and irrecoverable nodes and links combining them. Furthermore, the author composes a metamorphosis and evolution network consisting of six variations of the fundamental mathematical model in earthquake disaster areas and about non-systematic responses. Solving these equations qualitatively by isocline method, ruin or stability of the disaster and prevention systems is made clear by observing phase plane orbits with some coefficients. Finally, an immunity system is implied in order to follow the bigger earthquake disasters and systematic responses.