A cromodinâmica quântica (QCD) é a teoria das interações fortes, uma força fundamental que descreve a interação entre quarks e glúons que, por sua vez, constituem os hádrons como os prótons, nêutrons e píons. O modelo de Nambu–Jona-Lasinio (NJL) é uma teoria efetiva para a Cromodinâmica Quântica (QCD). O objetivo do projeto é o estudo do diagrama de fases da QCD sob campos magnéticos e elétricos intensos utilizando o modelo NJL. Para o nosso estudo, vamos analisar a matéria hadrônica em estrelas de nêutron, pois os seus campos elétrico e magnético são muito intensos. A partir das equações de Klein-Gordon, podemos encontrar o campo elétrico da estrela. Como as equações podem ter fontes não-lineares, resolvemos essas equações numericamente. Nisso, métodos de integração numérica e de solução de equações matriciais foram estudados.