Added ``power_analysis`` function