ベータ分布のスクリプトをscipy.statsで書きなおし