PUH-BrianMartell puh_rebound_asymmetry_factor_simulation.py- Updated New Py Code

puh_rebound_asymmetry_factor_simulation.py

import numpy as np
import matplotlib.pyplot as plt

# PUH v25: Rebound Asymmetry Factor Sim — Dipole Amplitude vs Asymmetry Factor
f_asym = np.linspace(1, 1.2, 500)  # Asymmetry factor arb.
dipole_amp = (f_asym - 1) * 0.01  # Amplitude toy linear

plt.figure(figsize=(10,6))
plt.plot(f_asym, dipole_amp, label='Dipole Amplitude vs Asymmetry', color='cyan', lw=2)
plt.axvline(1.06, color='gold', ls='--', label='Observed Asymmetry ~1.06')
plt.axhline(0.001, color='red', ls='--', label='Observed Dipole ~10^{-3}')
plt.xlabel('Rebound Asymmetry Factor f_{asym}'); plt.ylabel('Dipole Amplitude (arb.)')
plt.title('PUH v25: Rebound Asymmetry Factor Sim')
plt.legend(); plt.grid(alpha=0.3)
plt.tight_layout()
plt.savefig('puh_rebound_asymmetry_factor_simulation.png', dpi=300)
plt.show()

print("Linear dipole from asymmetry — converge observed at ~1.06.")