Добавил к вашему варианту пару строк в умножении, в итоге получилось 3 варианта для сравнения: как мне кажется точность их не очень различается, возможно где-то я допустил ошибку. Прикладываю исходники и графики зависимостей ошибки от степени полинома в логарифмическом масштабе.
Ну дисперсия поменьше будет, уже хорошо. Какой степени полиномы вы планируете интерполировать реализуемым методом?