Коэффициенты ряда Тейлора в МНК методе

[grandmix]

В общем суть такова . Мне надо сделать аппроксимацию методом наименьших квадратов , используя полином n-ой степени и вывести график этого полинома. Степень выбирается пользователем.
как известно полином имеет след.вид
a0+a1*x+a2*(x^2)+... и т.д.
я написал часть программы, которая находит коэффициенты перед иксами a0,a1..an
теперь надо вывести график. Я хочу это сделать, предварительно задав функцию, которая будет вычислять значения в зависимости от аргумента и по точечно строить его

float plot(double x) {
    return a0+a1*x+a2*(x^2)+..
}

тонкость в том, что степень полинома выбирается пользователем и следовательно число слагаемых a0+a1*x+a2*(x^2)+ будет всегда разным

Пример

Выбираю степень полинома 1, тогда функция должна возвращать

float plot(double x) {
    return a0+a1*x;
}

Выбираю степень полинома 3, тогда функция должна возвращать

float plot(double x) {
    return a0+a1*x+a2*(x^2)+a3*(x^3);
}

Вопрос : Как , имея заданный набор a0,a1,...an, число которых зависит от выбора пользователя, организовать соответсвующий вывод функции ???

float plot(double x) {
    return a0+a1*x+a2*(x^2)+..an+x^n
}

Помогите пожалуйста разобраться

[p51x]

res = 0;
for (int i = 0; i < n; ++i){
  res += a[i]*pow(x, i);
}
return res;

[grandmix]

спасибо большое !!