なお演算時間は半日かかる模様