Neste trabalho, realizamos um estudo numérico para o modelo cosmológico de Friedmann partindo de conceitos Newtonianos. Originalmente, esse modelo foi construído utilizando o ferramental matemático da Relatividade Geral, que possui grande nível de complexidade e não é acessível para os estudantes dos níveis iniciais da graduação. Sendo assim, a abordagem desse modelo através de conceitos Newtonianos torna acessível a compreensão de alguns conceitos básicos da Cosmologia. Além disso, também é possível que o estudante utilize esse material como uma introdução aos estudos dos métodos numéricos e computacionais. Na descrição de um Universo expansivo, as equações diferenciais são resolvidas numericamente pelo Método de Runge-Kutta de quarta ordem. Os resultados apresentados pelos gráficos das soluções numéricas condizem com os dados observacionais, descrevendo qualitativamente as eras pelas quais o Universo passou. Como motivação para o estudante, disponibilizamos as rotinas computacionais em C, Octave e Simulink.