A otimização de experimentos é um processo de reduz tempo, solventes, uso de equipamentos e melhora os resultados a partir de análise fatorial. A produção de géis fitoterápicos necessita de avaliação de estabilidade para poder ser utilizada de forma segura. Sendo assim, este estudo apresenta como objetivo otimizar por meio da análise fatorial do tipo Box-Behnken a produção de uma formulação de gel fitoterápico utilizando como variáveis dependentes os testes de densidade, espalhabilidade e pH obtidos durante a análise de estabilidade preliminar. As variáveis dependentes foram os percentuais do polímero de carbopol 940®, de óleo de manga e do extrato de chá verde. Após a execução do fatorial, foram desenvolvidas 13 formulações em triplicatas, submetidas a centrifugação e 6 ciclos de estresse térmicos, no qual demonstraram caraterísticas físico-químicas aceitáveis, em seguida, observou-se utilizando gráficos de superfície de resposta baixos níveis de aderência dos ativos na formulação quando haviam altos índices de densidade e espalhamento. O ponto ótimo da formulação foi composto pelos níveis polímero -1, extrato -1 e óleo +1, com desejabilidade 0,67. Este ponto se mostrou reprodutível na prática. Portanto, sugere-se que a análise fatorial é uma ferramenta importante para a otimização e garantia de formulações fitoterápicas de qualidade e estáveis.