Функция pow и график

//---------------------------------------------------------------------------
 
#include <vcl.h>
#include <math.h>
#pragma hdrstop
 
#include "Unit1.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma link "CSPIN"
#pragma resource "*.dfm"
#define gbp_cad 38
#define approx_point 30
#define prognoz_point 7
#define max_stepen_pol 5
TForm1 *Form1;
double A[gbp_cad][2], *a=NULL, *b=NULL, **sums=NULL, approx[approx_point+1], prognoz[prognoz_point+1];
int pol, N, K, r=0; //pol - степень аппроксимируещего полинома
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
        : TForm(Owner)
{
}
//---------------------------------------------------------------------------
void Refresh() //освобождаем массивы 
{for(int i=0; i<pol; i++){
       delete [] sums[i];
   }
   delete [] a;
   delete [] b;
   delete [] sums;
} 
 
//---------------------------------------------------------------------------
void memory_for_massiv() //выделяем память для массивов
{int i,j,k;
 a = new double[pol];
 b = new double[pol];
 sums = new double*[pol];
 for(i=0; i<pol; i++)
     {sums[i] = new double[pol];
      a[i]=0;
          b[i]=0;
      for(j=0; j<pol; j++)
          sums[i][j] = 0;
     }
}
//---------------------------------------------------------------------------
void set_koef() //заполняем коэфиценты системы уравнений и столбец свободных членов
{for(int i=0; i<pol; i++)
     {for(int j=0; j<pol; j++)
          {sums[i][j] = 0;
           for(int k=0; k<N; k++)
               sums[i][j] += pow(A[k][0], i+j);   }
     }
 for(int i=0; i<pol; i++)
     {for(int k=0; k<N; k++)
         b[i] += pow(A[k][0], i) * A[k][1];
     }
}
 
//---------------------------------------------------------------------------
void approximation() //аппроксимируем
{for (int i=0;i<N;i++)
     {approx[i]=0;
      for(int j=0; j<pol; j++)
          approx[i]+=a[j]*pow(A[i][0],j);
}
}
//---------------------------------------------------------------------------
void extrapolation() //прогнозируем
{for (int i=0;i<prognoz_point;i++)
     {prognoz[i]=0;
      for (int j=0;j<pol;j++)
          prognoz[i]+=a[j]*pow(A[i+N][0],j);
     }
}
//---------------------------------------------------------------------------
void Calculations() //все вычисления
{int i=0,j=0,k=0;
 double M,s;
 memory_for_massiv();
 set_koef();
 for(k=0; k<pol; k++)  // применяем метод гаусса для приведения к треугольному виду
     {for(i=k+1; i<pol; i++)
          {if(sums[k][k]==0)
               {ShowMessage("Solution is not exist.");
                return;
               }
           M = sums[i][k] / sums[k][k];
           for(j=k; j<pol; j++)
               {sums[i][j] -= M * sums[k][j]; 
               }
           b[i] -= M*b[k];
          }
     }
 for(i=K; i>=0; i--)
     {s = 0;
      for(j = i; j<pol; j++)
          s = s + sums[i][j]*a[j];
      a[i] = (b[i] - s) / sums[i][i];
     }
}
 
//---------------------------------------------------------------------------
void __fastcall TForm1::FormCreate(TObject *Sender)
{
CSpinEdit1->MinValue=1;
CSpinEdit1->MaxValue=max_stepen_pol;
Button3->Enabled=false;
Memo1->Lines->LoadFromFile("data.txt");
if (Memo1->Lines->Count<gbp_cad) {ShowMessage("Ваш файл содержит меньше 38 точек"); exit(1);}
Memo1->Visible=false;
for (int n=0;n<+37;n++)
{
A[n][0]=n;
A[n][1]=StrToFloat(Memo1->Lines->Strings[n]); //заполняем массив из txt файла
}
for (int i=0;i<=30;i+=2)
Chart1->Series[0]->AddXY(A[i][0],A[i][1]);
N=approx_point+1;
}
//---------------------------------------------------------------------------
 
void __fastcall TForm1::Button1Click(TObject *Sender) //кнопка аппроксимации
{
Button3->Enabled=true;
 {Chart1->Series[1]->Clear();
      Chart1->Series[2]->Clear();
          Chart1->Series[3]->Clear();
          Refresh();
     }
 K=CSpinEdit1->Value;
 pol=K+1;
Calculations(); 
approximation();
for (int i=0;i<N;i++) 
     Chart1->Series[1]->AddXY(A[i][0],approx[i]);
r=1;
}
//---------------------------------------------------------------------------
 
void __fastcall TForm1::Button2Click(TObject *Sender) //скользящая средняя
{
float k=(StrToFloat(Memo1->Lines->Strings[26])+StrToFloat(Memo1->Lines->Strings[28])+StrToFloat(Memo1->Lines->Strings[30]))/3;
float p=(StrToFloat(Memo1->Lines->Strings[28])+StrToFloat(Memo1->Lines->Strings[30])+k)/3;
float o=(StrToFloat(Memo1->Lines->Strings[30])+k+p)/3;
float n=(k+p+o)/3;
float z=(o+n)/2;
Chart1->Series[2]->AddXY(30,StrToFloat(Memo1->Lines->Strings[30]));
Chart1->Series[2]->AddXY(32,k);
Chart1->Series[2]->AddXY(34,p);
Chart1->Series[2]->AddXY(36,o);
Chart1->Series[2]->AddXY(37,z);
}
//---------------------------------------------------------------------------
void __fastcall TForm1::Button3Click(TObject *Sender) //дорисовываем прогноз
{
 {extrapolation();
Chart1->Series[3]->AddXY(A[approx_point][0],approx[approx_point]);
 
for (int i=0;i<prognoz_point;i++)
     Chart1->Series[3]->AddXY(A[i+N][0],prognoz[i]);
}
 
}
//---------------------------------------------------------------------------