У меня много образцов(y_i, (a_i, b_i, c_i)) гдеy Предполагается, что варьируется как полином вa,b,c  до определенной степени. Например, для заданного набора данных и степени 2 я мог бы создать модель

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

Это может быть выполнено с использованием метода наименьших квадратов и является небольшим продолжением рутины nimpy 'polyfit. Есть ли стандартная реализация где-нибудь в экосистеме Python?