Tengo muchas muestras(y_i, (a_i, b_i, c_i)) dóndey se presume que varía como un polinomio ena,b,c hasta cierto punto Por ejemplo, para un conjunto dado de datos y grado 2 podría producir el modelo

y = a^2 + 2ab - 3cb + c^2 +.5ac

Esto se puede hacer usando mínimos cuadrados y es una ligera extensión de la rutina polyfit de numpy. ¿Hay alguna implementación estándar en algún lugar del ecosistema de Python?