A sub-beam approximation step for approximating a particle beam with a set of a plurality of sub-beams each having a Gaussian distribution, and simulating a state in which each sub-beam of the plurality of sub-beams is deflected by a scanning device and progresses, thereby making each sub-beam a patient A sub-beam dose distribution calculating step for calculating a dose distribution formed inside the patient by calculating each sub-beam dose distribution formed inside and integrating the calculated sub-beam dose distributions is provided.粒子線を、それぞれガウス分布を有する複数のサブビームの集合で近似するサブビーム近似ステップと、複数のサブビームのそれぞれのサブビームが走査装置により偏向されて進行する状態をシミュレーションすることにより、それぞれのサブビームが患者内部に形成するそれぞれのサブビーム線量分布を演算して、演算したそれぞれのサブビーム線量分布を積算することにより粒子線が患者内部に形成する線量分布を求めるサブビーム線量分布演算ステップとを備えるようにした。
