﻿import numpy as np  # モジュールNumpyをnpという名前でインポート
from scipy import stats  # 高度な統計分析を行うためのサブモジュール
z = stats.norm.ppf(0.95, loc = 0, scale = 1)
z  # z_0.05を表示
mu = 1.0
print(stats.norm.sf(z, loc=np.sqrt(33)*mu/2, scale=1))  # 検出力(n=33)
print(stats.norm.sf(z, loc=np.sqrt(34)*mu/2, scale=1))  # 検出力(n=34)
print(stats.norm.sf(z, loc=np.sqrt(35)*mu/2, scale=1))  # 検出力(n=35)
print(stats.norm.sf(z, loc=np.sqrt(36)*mu/2, scale=1))  # 検出力(n=36)