Линеаризация функции полиномом н-ого порядка

[DreamForLives]

Мне нужно написать программу: "Линеаризация функции 2ух переменных полиномом н-ого порядка".

Но если честно, никак не могу понять как к ней подойти ...даже сам ход мыслей её написания.

п.с. пока это задание смогла написать только для полинома 2ого,3тьего порядка...но от меня всё равно требуют н-ого(

[DreamForLives]

#include "stdafx.h"
#include "conio.h"
#include <cmath>

float a1,a2,a3,a4;
int const j=1;

float massiv[j];

int _tmain(int argc, _TCHAR* argv[])
{
int l,h,i=0,n1,n2,n3,n4,k1,k2,k3,k4;
printf ("vvedite kolichestvo tochek, ne menee 5 dlia x^2 i ne menee 8 dlia x^3 \n");
scanf ("%d", &h);
printf ("dannie tochki\n");
for (int l=0;l<h;l++)
{
scanf ("%f", &massiv[l]);
}
printf ("\nintervali koefficientov\n\n");
printf ("a1 ot "); scanf ("%d", &n1); printf ("a1 do "); scanf ("%d", &k1);
printf ("\na2 ot "); scanf ("%d", &n2); printf ("a2 do "); scanf ("%d", &k2);
printf ("\na3 ot "); scanf ("%d", &n3); printf ("a3 do "); scanf ("%d", &k3);
printf ("\na4 ot "); scanf ("%d", &n4); printf ("a4 do "); scanf ("%d", &k4);
for (float a1=n1;a1<k1;a1++)
{
for (float a2=n2;a2<k2;a2++)
{
for (float a3=n3;a3<k3;a3++)
{
if
(pow(massiv[i]-(a1*i*i+a2*i+a3),2)<0.1
&&
pow(massiv[i+1]-(a1*(i+1)*(i+1)+a2*(i+1)+a3),2)<0.1 &&
pow(massiv[i+2]-(a1*(i+2)*(i+2)+a2*(i+2)+a3),2)<0.1
&&
pow(massiv[i+3]-(a1*(i+3)*(i+3)+a2*(i+3)+a3),2)<0.1
&&
pow(massiv[i+4]-(a1*(i+4)*(i+4)+a2*(i+4)+a3),2)<0.1
&&
a1>0 || a1<0)
{
printf ("\npolinom 2 stepeni a1*x*x+a2*x*a3\n\n");
printf ("a1 %f \na2 %f \na3 %f\n", a1,a2,a3);
printf ("\n");
}
if
(pow(massiv[i]-(a1*i*i+a2*i+a3),2)<0.1
&&
pow(massiv[i+1]-(a1*(i+1)*(i+1)+a2*(i+1)+a3),2)<0.1
&&
pow(massiv[i+2]-(a1*(i+2)*(i+2)+a2*(i+2)+a3),2)<0.1
&&
pow(massiv[i+3]-(a1*(i+3)*(i+3)+a2*(i+3)+a3),2)<0.1
&&
pow(massiv[i+4]-(a1*(i+4)*(i+4)+a2*(i+4)+a3),2)<0.1
&&
a1==0)
{
printf ("\npolinom 1 stepeni a1*x+a2\n");
printf ("a1 %f \na2 %f\n", a2,a3);
printf ("\n");
}
for (float a4=n4;a4<k4;a4++)
{
if
(pow(massiv[i]-(a1*i*i*i+a2*i*i+a3*i+a4),2)<0.1
&&
pow(massiv[i+1]-(a1*(i+1)*(i+1)*(i+1)+a2*(i+1)*(i+1 )+a3*(i+1)+a4),2)<0.1
&&
pow(massiv[i+2]-(a1*(i+2)*(i+2)*(i+2)+a2*(i+2)*(i+2 )+a3*(i+2)+a4),2)<0.1
&&
pow(massiv[i+3]-(a1*(i+3)*(i+3)*(i+3)+a2*(i+3)*(i+3 )+a3*(i+3)+a4),2)<0.1
&&
pow(massiv[i+4]-(a1*(i+4)*(i+4)*(i+4)+a2*(i+4)*(i+4 )+a3*(i+4)+a4),2)<0.1
&&
pow(massiv[i+5]-(a1*(i+5)*(i+5)*(i+5)+a2*(i+5)*(i+5 )+a3*(i+5)+a4),2)<0.1
&&
pow(massiv[i+6]-(a1*(i+6)*(i+6)*(i+6)+a2*(i+6)*(i+6 )+a3*(i+6)+a4),2)<0.1
&&
pow(massiv[i+7]-(a1*(i+7)*(i+7)*(i+7)+a2*(i+7)*(i+7 )+a3*(i+7)+a4),2)<0.1
&&
a1>0 || a1<0)
{
printf ("\npolinom 3 stepeni a1*x*x*x+a2*x*x+a3*x+a4\n\n");
printf ("a1 %f \na2 %f \na3 %f\na4 %f \n", a1,a2,a3,a4); printf ("\n");
}
}
}
}
}
getch();
return 0;
}