やっぱり計算を間違えていたかも知れない。 修正後は こうなった。違いは散乱位置を d = -log(rnd()) / σt ただし 0 < d < distance としていたのを d = -log(1-rnd()* (1-exp(-σt * distance))) / σt とした。当然だけど確率分布(probability distribution)…