import numpy as npDF=4.66920160910299; PI=np.pig=(1-np.sqrt(1-4*DF*PI*np.sqrt(2)/90))/2; z=g*(1-g)# P1 = E1の代数的帰結 (厳密証明)assert abs(np.pi/(90*z)**2*DF**2 - 1/(2*PI)) < 1e-14 # 厳密成立# 射影変換定理C = np.sqrt(2)*PI**5/8100inv_a = (90*z)**2/PI Cprint(f"α_ARK·δ_F² = {np.pi/(90*z)**2*DF**2:.15f}")print(f"1/(2π) = {1/(2*PI):.15f}")print(f"1/α = {inv_a:.8f} (誤差: {(inv_a-137.035999084)/137.035999084*100:.2e}%)")