Por medio de los parámetros de densidad:
$$
\Omega_{DM}=0.25\approx 25\%
$$$$
\Omega_{B}=0.05\approx 5\%
$$$$
\Omega_{\nu}=10^{-3} \approx 0.1\%
$$$$
\Omega_{\gamma}=10^{-4} \approx 0.01\%
$$$$
\Omega_{\Lambda}=0.6989 \approx 69.89\%
$$$$
H_{0}=1
$$Los valores de estos parámetros ya se pueden actualizar a los datos 2015 de Planck.
Para que Julia trabaje en un modo compilación (JIT), debemos de considerar en lo más mínimo variables globales. Entonces definimos una función que tenga como entradas: $a,b \in \mathbb{R}$ que son nuestros puntos iniciales y finales. $N \in \mathbb{N}$ es el número total de puntos a graficar. Se tiene que las ecuaciones a graficar son: $$H=\left(x_0^2 e^{-3N}+y_0^2 e^{-3N}+l_0^2 e^{-4N}+m_0^2 e^{-4N}+C\right)^{\frac{1}{2}}$$ $$H(0)=1=x_0^2 +y_0^2 +l_0^2 +m_0^2 +C// C= 1 -x_0^2+y_0^2+l_0^2+m_0^2=z_0^2 $$
$$H=\left(x_0^2 e^{-3N}+y_0^2 e^{-3N}+l_0^2 e^{-4N}+m_0^2 e^{-4N}+z_0^2\right)^{\frac{1}{2}}$$y para cada componente es: $$ x=\dfrac{x_0}{H}e^{-\frac{3}{2}N}\\ y=\dfrac{y_0}{H}e^{-\frac{3}{2}N}\\ z=x=\dfrac{z_0}{H}\\ l=\dfrac{l_0}{H}e^{-2N}\\ m=\dfrac{m_0}{H}e^{-2N} $$
$$ x^2=\dfrac{x_0^2}{H^2}e^{-3N}=\Omega_{DM}\\ y^2=\dfrac{y_0^2}{H^2}e^{-3N}=\Omega_B\\ z^2=x=\dfrac{z_0^2}{H^2}=\Omega_\Lambda\\ l^2=\dfrac{l_0^2}{H^2}e^{-4N}=\Omega_\nu\\ m^2=\dfrac{m_0^2}{H^2}e^{-4N}=\Omega_\gamma $$Defino la función para $H$.
In [18]:
0.268+0.048+1e-3+1e-4+0.6989
Out[18]:
In [19]:
function Hache(N::Float64)
x0=0.268
y0=0.048
l0=1e-3
m0=1e-4
z0=0.685
H0=1.0
H = x0*exp(-3.0*N)+ y0*exp(-3.0*N) + l0*exp(-4.0*N) + m0*exp(-4.0*N) + z0
return H
end
Out[19]:
In [20]:
1-(0.268+0.048+1e-3+1e-4+0.685)
Out[20]:
Ahora una función para cada una de mis componentes
Se define para $\Omega_B$ y $\Omega_{DE}$ y C0 determina la entrada del valor de la componentes
In [21]:
function f1(fH,C0,a,b,p)
x=Array{Float64,1}[]
Nt=(Array{Float64,1})[]
for i in a:p:b
x_n = C0*exp(-3.0i)
x_nn = x_n/fH(i)
push!(x,[x_nn])
# push!(Nt,[i])
end
x
end
Out[21]:
Se define para la $\Omega_{rad}$ y $\Omega_\nu$, y también C0 determina el valor de la respectiva componente
In [22]:
function f2(fH,C0,a,b,p)
x=Array{Float64,1}[]
Nt=(Array{Float64,1})[]
for i in a:p:b
x_n = C0*exp(-4.0*i)/fH(i)
x_nn = x_n
push!(x,[x_nn])
# push!(Nt,[i])
end
x
end
Out[22]:
Se define la funcón que depende de $\Omega_\Lambda$
In [23]:
function f3(fH,C0,a,b,p)
x=Array{Float64,1}[]
Nt=(Array{Float64,1})[]
for i in a:p:b
x_n = C0/fH(i)
x_nn = x_n
push!(x,[x_nn])
# push!(Nt,[i])
end
x
end
Out[23]:
In [24]:
using PyPlot
Se gráfican las 5 componentes del Universo
In [26]:
a=-20.0; b=5.0; p=0.1;
DMs=f1(Hache,0.268,a,b,p);
Bs=f1(Hache,0.048,a,b,p);
Rs=f2(Hache,1e-3,a,b,p);
Ns=f2(Hache,1e-4,a,b,p);
CCs=f3(Hache,0.685,a,b,p);
ts=[a:p:b];
plt1 = plot(ts, DMs,label=L"$\Omega_{DM}$");
plt2 = plot(ts, Bs,label=L"$\Omega_B$");
plt3 = plot(ts, CCs,label=L"$\Omega_\Lambda$");
plt4 = plot(ts, Rs,label=L"$\Omega_\gamma$");
plt5 = plot(ts,Ns,label=L"$\Omega_\nu$")
ylim(0,1.01)
xlim(-10,2)
grid("on");
ylabel(L"$\Omega$");
xlabel(L"$N=log(a)$")
legend(handles=[plt1,plt2,plt3,plt4,plt5],loc=6);
axvline(0,color="gray");
show()
In [ ]:
In [ ]:
In [ ]: