Em 1998, observações de supernovas do tipo Ia constataram que o nosso Universo está em expansão acelerada. Para explicar esse fenômeno, foi introduzida uma componente exótica com pressão negativa chamada energia escura. Essa inclusão é feita no contexto do Modelo Cosmológico Padrão, assumindo o Princípio Cosmológico e a teoria da Relatividade Geral de Einstein. Neste cenário, o objetivo deste trabalho é reconstruir a equação de estado da energia escura, w(z)w(z), assumindo-a como um fluido barotrópico, usando uma abordagem paramétrica e independente de modelo, de modo que apenas os dados determinem a reconstrução. O método utilizado baseia-se em splines cúbicas, comumente aplicadas em interpolação polinomial. As análises foram divididas em três etapas, incorporando dados observacionais cada vez mais abrangentes e com mudanças no intervalo de redshift. Na terceira etapa, aplicamos o método de Monte Carlo via Cadeias de Markov (MCMC) para obter estimativas mais confiáveis das incertezas, superando as limitações da matriz de Fisher. Como resultado, para redshifts baixos (0 ≤ z ≤ 3), observamos um comportamento bimodal para w(z)w(z), com picos em w=−1w=−1 e em 0<w<1/30<w<1/3. Esse comportamento é compatível com a literatura, que também aponta valores próximos de w=−1w=−1, como w^=−1.087±0.03592w^=−1.087±0.03592, consistentes com a constante cosmológica.