Ich habe viele Proben(y_i, (a_i, b_i, c_i)) wohery wird als Polynom angenommena,b,c bis zu einem gewissen Grad. Zum Beispiel könnte ich für einen gegebenen Datensatz und Grad 2 das Modell erstellen

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

Dies kann mit den kleinsten Quadraten durchgeführt werden und ist eine leichte Erweiterung der Polyfit-Routine von numpy. Gibt es irgendwo im Python-Ökosystem eine Standardimplementierung?