Eu tenho muitas amostras(y_i, (a_i, b_i, c_i)) Ondey presume-se que varie como um polinômio ema,b,c até certo ponto. Por exemplo, para um dado conjunto de dados e grau 2 eu poderia produzir o modelo

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

Isso pode ser feito usando mínimos quadrados e é uma pequena extensão da rotina polyfit numpy. Existe uma implementação padrão em algum lugar no ecossistema Python?