import numpy as np import matplotlib.pyplot as plt # 1. Define observational and model data (Simulated based on NGC 2903 context) R_obs = np.array([0.5, 1.0, 2.0, 3.0, 5.0, 8.0, 12.0, 16.0, 20.0, 25.0]) V_obs = np.array([120, 160, 195, 210, 215, 212, 210, 208, 205, 202]) # Density-Corrected V_FRCFD (where the core tracks the high baryonic baseline) V_FRCFD = np.array([140, 185, 230, 250, 260, 245, 210, 205, 202, 200]) # 2. Calculate Residuals (Observed - Predicted) residuals = V_obs - V_FRCFD # 3. Generate Plot plt.figure(figsize=(10, 5)) plt.scatter(R_obs, residuals, color='red', label='Residuals (V_obs - V_FRCFD)') plt.axhline(0, color='black', linestyle='--', alpha=0.5) # 4. Apply Labels and Scaling plt.title('Residuals — Density-Corrected Blind FRCFD — NGC 2903') plt.xlabel('Radius R (kpc)') plt.ylabel('Residual (V_obs - V_FRCFD) [km/s]') # Setting requested y-axis limits plt.ylim(-150, 150) plt.grid(True, alpha=0.3) plt.legend() plt.show()

Popular posts from this blog

THE GOLDEN BALLROOM/BUNKER

ICE PROUDLY ANNOUNCES NEW “ELITE” TASK FORCE COMMANDER JEREMY DEWITTE

Conceptual Summary #2: (∂t2​S−c2∇2S+βS3)=σ(x,t)⋅FR​(C[Ψ])