Mam wiele próbek(y_i, (a_i, b_i, c_i)) gdziey zakłada się, że zmienia się jako wielomiana,b,c do pewnego stopnia. Na przykład dla danego zestawu danych i stopnia 2 mogę wyprodukować model

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

Można to zrobić za pomocą najmniejszych kwadratów i jest to niewielkie rozszerzenie procedury polyfit numpy. Czy gdzieś w ekosystemie Pythona istnieje standardowa implementacja?