import numpy as np

# 1. Parámetros físicos unificados del modelo VCSD-I
G, c, rho_0, s_0, sigma_0, pasos = 1.0, 1.0, 1.0, 0.5, 0.1, 300
a_c = 1.252838  
factores_escala = np.linspace(3.0, a_c, pasos)

# 2. Bucle en una sola línea continua (Imposible de fallar por sangrado)
hist_H2 = [(8 * np.pi * G / 3) * (rho_0 / a**4) - (16 * np.pi**2 * G**2 / (3 * c**2)) * (s_0**2 / a**6) + (sigma_0**2 / a**6) for a in factores_escala]
hist_densidad = [rho_0 / a**4 for a in factores_escala]

# 3. Despliegue de los resultados reales de la inversión métrica
print("\n" + "=" * 75)
print("     DINÁMICA PRIMORDIAL VALIDADA (VCSD-I) - RESULTADOS NUMÉRICOS")
print("=" * 75)
print(f"Inicio Colapso  -> Factor Escala: {factores_escala[0]:.4f}  | Densidad: {hist_densidad[0]:.4f} | H²: {hist_H2[0]:.4f}")
print(f"Paso Intermedio -> Factor Escala: {factores_escala[150]:.4f}  | Densidad: {hist_densidad[150]:.4f} | H²: {hist_H2[150]:.4f}")
print(f"Paso de Rebote  -> Factor Escala: {factores_escala[-1]:.4f}  | Densidad: {hist_densidad[-1]:.4f} | H²: {hist_H2[-1]:.4f}")
print("=" * 75)