No presente estudo foi estimada a curva de lactação média de 7239 lactações de vacas da raça Caracu, com base em 299.949 controles semanais. As funções ajustadas aos dados médios da curva de lactação foram Linear Hiperbólica e Quadrática Logarítmica, como lineares, e Gama Incompleta e Polinomial Inversa, como não-lineares. Os parâmetros foram estimados por meio de processos iterativos, método de Gauss-Newton. O Coeficiente de Determinação ajustado (R²A), o Teste de Durbin-Watson (DW) e os desvios (%) entre as produções totais (PTL) observadas e estimadas foram usados como medida da qualidade do ajuste. Os resíduos estimados, apesar de pequenos, apresentaram correlação serial positiva, pois o teste DW foi significativo para as quatro funções estudadas, indicando violação da pressuposição de independência entre estes resíduos. Os R²A foram superiores a 97% e os desvios entre as PTL observadas e estimadas, mínimos, indicando que, de acordo com estes critérios, qualquer função poderia ser usada para ajustar a curva de lactação média desse rebanho. Quanto aos componentes da curva, a melhor aproximação foi obtida pela função Polinomial Inversa, para produção inicial, tempo e produção no pico. As funções Gama Incompleta e Quadrática Logarítmica apresentaram ajustes semelhantes, porém a primeira apresentou melhor definição dos parâmetros e componentes da curva de lactação.