To propose a method and an analysis program for predicting the molecular motion and alleviation/flow phenomenon of high molecular solution in which solid particles are dispersed.

As for tube configuring elements configuring respective high molecules and particles (solid particles) generated by a system generation means 114, a system motion calculation means 115 in a molecular motion analysis device 1 performs a numerical analysis process for temporally developing position coordinates, alignment vectors, the number of monomers in a tube configuring elements and the other party tube configuring elements forming tangle according to a Langevin equation. The solid particles can be preliminarily arranged in a certain place of this system, or generated by using the phase separation of high molecules to be successively calculated by the numerical analysis process. The position of the tube configuring element of the solid particle section in the created system is temporally developed according to Formula 9.