По коду необходимо создать EXE-файл

#include <vcl.h>
#pragma hdrstop
 
#include "Unit1.h"
#include <math.h>
#define sz 11
#define cap 12
//---------------------------------------------------------------------
#pragma package(smart_init)
#pragma link "CSPIN"
#pragma resource "*.dfm"
TForm1 *Form1;
void Add(long double,long double,long double);
void Graph();
void Polygraph(int n,double cc[]);
void line();
void poly();
void exp();
void Axis();
//=====================================================================
float xy[2][sz]={  {2 ,2.5,3.4,4.4,4.7,5.3,5.7,6.5,7.2,7.6,8.1},
                {27.325,35.558, 30.607, 48.505, 55.289, 57.7120, 68.755,62.211,70.334,78.661,80.55}};
long double c[cap];
long double mtr[cap][cap];
//---------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
        : TForm(Owner)
{
}
void __fastcall TForm1::Button1Click(TObject *Sender)
{
line();
}
void __fastcall TForm1::FormShow(TObject *Sender)
{
Form1->StringGrid1->Rows[0]->Strings[1]="X";
Form1->StringGrid1->Rows[0]->Strings[2]="Y";
Form1->StringGrid1->Rows[0]->Strings[3]="Q";
Form1->StringGrid1->Rows[0]->Strings[4]="y - Q";
}
void __fastcall TForm1::Button2Click(TObject *Sender)
{
poly();}
 
void __fastcall TForm1::Button3Click(TObject *Sender)
{
Form1->PaintBox1->Canvas->Brush->Color=clWhite;
Form1->PaintBox1->Canvas->FillRect(ClientRect);
}
void __fastcall TForm1::Button4Click(TObject *Sender)
{
exp();
}
//Функция, добавляющая значения функции и аппроксимирующей функции в таблицу на экране
void Add(long double x,long double y, long double Q){
int cnt= Form1->StringGrid1->RowCount;
Form1->StringGrid1->RowCount++;
Form1->StringGrid1->Rows[cnt]->Strings[1]=FloatToStr(x);
Form1->StringGrid1->Rows[cnt]->Strings[2]=FloatToStr(y);
Form1->StringGrid1->Rows[cnt]->Strings[3]=FloatToStr(Q);
 
Form1->StringGrid1->Rows[cnt]->Strings[4]= FloatToStr(y-Q);
 
}
void Axis(){
Form1->StringGrid1->FixedRows=1;
//-------------draw axis------------------
Функция чертит оси координат в поле вывода
Form1->PaintBox1->Canvas->MoveTo(0,Form1->PaintBox1->Height-10);//Y
Form1->PaintBox1->Canvas->LineTo(Form1->PaintBox1->Width,Form1->PaintBox1->Height-10);
 
Form1->PaintBox1->Canvas->MoveTo(Form1->PaintBox1->Left+10,0); //X
Form1->PaintBox1->Canvas->LineTo(Form1->PaintBox1->Left+10,Form1->PaintBox1->Height);
}
//Функция строит график по табличным значениям и рассчитывает среднеквадратическое отклонение
void Graph(){
double x,y,Q,sx,sy,boty,QQ=0,xx=0;
sx=50;
sy=5;
double sig=0;
Axis();
boty=Form1->PaintBox1->Height-10;
 Form1->PaintBox1->Canvas->Brush->Color=clGreen;
 
for (int i=1;i<Form1->StringGrid1->RowCount;i++){
 x=StrToFloat(Form1->StringGrid1->Rows[i]->Strings[1]);
 y=StrToFloat(Form1->StringGrid1->Rows[i]->Strings[2]);
 Q=StrToFloat(Form1->StringGrid1->Rows[i]->Strings[3]);
 
 sig+=pow(StrToFloat(Form1->StringGrid1->Rows[i]->Strings[4]),2);
 //--------------------draw---------------------
 if ((QQ==0)&&(xx==0)) {QQ=Q;xx=x;}
 Form1->PaintBox1->Canvas->Pen->Color=clGreen;
 Form1->PaintBox1->Canvas->MoveTo(xx*sx,boty-QQ*sy);Form1->PaintBox1->Canvas->LineTo(x*sx,boty-Q*sy);
  Form1->PaintBox1->Canvas->Pen->Color=clBlue;
  Form1->PaintBox1->Canvas->Ellipse(x*sx-2,boty-y*sy-2,x*sx+2,boty-y*sy+2);
 QQ=Q;xx=x;
}
sig=sqrt(sig/Form1->StringGrid1->RowCount);
Form1->StatusBar1->Panels->Items[1]->Text="6= "+FloatToStr(sig);
}
//Функция построения графика полинома. Рассчитывает значения полинома в промежуточных точках и чертит график. Входные параметры степень полинома, указатель на массив коэффициентов.
void Polygraph(int n,long double cc[]){
double x,y,Q,sx,sy,boty,QQ=0,xx=0;
sx=50;
sy=5;
double sig=0;
Axis();
boty=Form1->PaintBox1->Height-10;
 Form1->PaintBox1->Canvas->Brush->Color=clGreen;
 
for (int i=1;i<Form1->StringGrid1->RowCount;i++){
 x=StrToFloat(Form1->StringGrid1->Rows[i]->Strings[1]);
 y=StrToFloat(Form1->StringGrid1->Rows[i]->Strings[2]);
 Q=StrToFloat(Form1->StringGrid1->Rows[i]->Strings[3]);
 
 sig+=pow(StrToFloat(Form1->StringGrid1->Rows[i]->Strings[4]),2);
 //--------------------draw---------------------
 if ((QQ==0)&&(xx==0)) {QQ=Q;xx=x;}
  Form1->PaintBox1->Canvas->Pen->Color=clBlue;
  Form1->PaintBox1->Canvas->Ellipse(x*sx-2,boty-y*sy-2,x*sx+2,boty-y*sy+2);
 QQ=Q;xx=x;
}
 Form1->PaintBox1->Canvas->Pen->Color=clGreen;
//(((((((((((((((((((((((((((((((((((((((((
QQ=0;
for (float x=1;x<10;x+=0.01) {
Q=0;
    for (int i=0;i<=n;i++){
      Q+=cc[i]*pow(x,i);
    }
    if (QQ==0) Form1->PaintBox1->Canvas->MoveTo(x*sx,boty-sy*Q);
    Form1->PaintBox1->Canvas->LineTo(x*sx,boty-sy*Q); QQ=Q;
    }
    //)))))))))))))))))))))))))))))))))))))))))
sig=sqrt(sig/Form1->StringGrid1->RowCount);
Form1->StatusBar1->Panels->Items[1]->Text="6= "+FloatToStr(sig);
}
//Функция чертит график степной функции и рассчитывает ее среднеквадратичное отклонение
void expgraph(long double c0,long double c1){
double x,y,Q,sx,sy,boty,QQ=0,xx=0;
sx=50;
sy=5;
double sig=0;
Axis();
boty=Form1->PaintBox1->Height-10;
 Form1->PaintBox1->Canvas->Brush->Color=clGreen;
 
for (int i=1;i<Form1->StringGrid1->RowCount;i++){
 x=StrToFloat(Form1->StringGrid1->Rows[i]->Strings[1]);
 y=StrToFloat(Form1->StringGrid1->Rows[i]->Strings[2]);
 Q=StrToFloat(Form1->StringGrid1->Rows[i]->Strings[3]);
 
 sig+=pow(StrToFloat(Form1->StringGrid1->Rows[i]->Strings[4]),2);
 //--------------------draw---------------------
 if ((QQ==0)&&(xx==0)) {QQ=Q;xx=x;}      
  Form1->PaintBox1->Canvas->Pen->Color=clBlue;
  Form1->PaintBox1->Canvas->Ellipse(x*sx-2,boty-y*sy-2,x*sx+2,boty-y*sy+2);
 QQ=Q;xx=x;
}
 Form1->PaintBox1->Canvas->Pen->Color=clGreen;
//(((((((((((((((((((((((((((((((((((((((((
QQ=0;
for (float x=1;x<10;x+=0.01) {
Q=0;
         Q=c0*pow(x,c1);
 
        if (QQ==0) Form1->PaintBox1->Canvas->MoveTo(x*sx,boty-sy*Q);
    Form1->PaintBox1->Canvas->LineTo(x*sx,boty-sy*Q); QQ=Q;
    }
    //)))))))))))))))))))))))))))))))))))))))))
sig=sqrt(sig/Form1->StringGrid1->RowCount);
 Form1->StatusBar1->Panels->Items[0]->Text="Y(x)= "+FloatToStr(c0)+"*x^ "+FloatToStr(c1);
 Form1->StatusBar1->Panels->Items[1]->Text="6= "+FloatToStr(sig);
}
//=====================================================================
//Реализация метода Гаусса для решения СЛАУ. Входные параметры, размерность, матрица значений, вектора аргумента и свободных членов.
void gauss(int h,long double a[][cap],long double *x,long double *b){
long double tmp;
for (int k=0;k<h-1;k++){
 for (int i=k+1;i<h;i++) {
 tmp=a[i][k]/a[k][k];
 b[i]-=b[k]*tmp;
    for (int j=0;j<h;j++){
      a[i][j]-=a[k][j]*(tmp);
    } } }//i j k
 for(int i=h-1;i>=0;i--){
 x[i]=b[i];
    for (int j=i+1;j<h;++j)
    {
      x[i]-=a[i][j]*x[j];
    }
    x[i]/=a[i][i];}
}
//Функция  рассчитывает коэффициенты  для аппроксимирующей 
void line(){ степенной  функции  методом сведения к линейной.
int lvl=4,i,j;
Form1->StringGrid1->RowCount=1;
//==-swap to linear aproximation-==
long double a[2][12],b[2],c[2];
for (int i=0;i<2;i++) for (int j=0;j<2;j++){ a[i][j]=0;b[i]=0;c[i]=0;}
for (int i=0;i<sz;i++){
      a[0][0]+=pow(log(xy[0][i]),2); a[0][1]+=(log(xy[0][i]));     b[0]+=(log(xy[0][i]))*log(xy[1][i]);
     a[1][0]+=(log(xy[0][i]));               a[1][1]+=1.0;                  b[1]+=log(xy[1][i]);
}
gauss (2,a,c,b);
for (int i=0;i<sz;i++) Add(xy[0][i],xy[1][i],exp(c[1])*pow(xy[0][i],c[0]));
 
Graph();
Form1->StatusBar1->Panels->Items[0]->Text="Y(x)= "+FloatToStr(exp(c[1]))+"*x^ "+FloatToStr(c[0]);
}
//Функция  состовляет матрицу и вектора для дальнейшего расчета коэффициентов  полинома по методу Гаусса.
void poly(){
Form1->StringGrid1->RowCount=1;
int n=Form1->CSpinEdit1->Value;
long double a[12][12],b[12];
for (int i=0;i<12;i++)for(int j=0;j<12;j++) {a[i][j]=0;b[i]=0;c[i]=0;}
for(int i=0;i<=n;i++)
        for (int j=0;j<=n;j++)
                for (int m=0;m<sz;m++){
                        a[i][j]+=pow(xy[0][m],i+j);
                        if (j==0)b[i]+=xy[1][m]*pow(xy[0][m],i);
                }
gauss(n+1,a,c,b);
long double Q=0;
 
Form1->StatusBar1->Panels->Items[0]->Text="Y(x)= ";
for (int i=0;i<sz;i++){
        for (int j=0;j<=n;j++) Q+=c[j]*pow(xy[0][i],j);
 Add(xy[0][i],xy[1][i],Q); Q=0;
 }
 for (int j=0;j<=n;j++){
 Form1->StatusBar1->Panels->Items[0]->Text+=FloatToStr(c[j]);
 if (j>0)Form1->StatusBar1->Panels->Items[0]->Text+="*x^"+IntToStr(j);
 if(j<n) Form1->StatusBar1->Panels->Items[0]->Text+=" + ";
 }
Polygraph(n,c);
}
//Функция рассчитывает значение первого (второго) уравнения системы.
double func1(double x,double y){
  double deb =0;
  for (int i=0;i<sz;i++) deb+=(xy[1][i]-x*pow(xy[0][i],y))*pow(xy[0][i],y);
  return deb;
}
double func2(double x,double y){
  double deb= 0;
  for (int i=0;i<sz;i++)
  deb+=(xy[1][i]-x*pow(xy[0][i],y))*
  x*pow(xy[0][i],y)*log(xy[0][i]);
  return deb;
}
//Функция рассчитывает и составляет  матрицу коэффициентов (Якобиан) для системы уравнений.
void jac(long double c0,long double c1,long double res[][cap]){
   long double a11=0,a12=0,a21=0,a22=0;
 
 for (int i=0;i<sz;i++){
 a11+=-pow(xy[0][i],2*c1);a12+=pow(xy[0][i],c1)*log(xy[0][i])*(xy[1][i]-2*c0*pow(xy[0][i],c1));
 a21=a12;                a22+=c0*powl(logl(xy[0][i]),2)*powl(xy[0][i],c1)*(xy[1][i]-c0*2*powl(xy[0][i],c1));
 }
 res[0][0]=a11;res[0][1]=a12;
 res[1][0]=a21;res[1][1]=a22;
}
//Функция рассчитывает коэффициенты степной функции в явной форме по методу Ньютона.
 
void exp(){
Form1->StringGrid1->RowCount=1;
        long double x0=1;long double y0=1;
   long double xx[2],rr[3];
   double eps=1e-5;
   long double b,cf[4];
   int num=0;
   x0=15;y0=1;
   do{
   xx[0]=-func1(x0,y0);xx[1]=-func2(x0,y0);
   jac(x0,y0,mtr);
   gauss(2,mtr,rr,xx);
   x0+=rr[0];y0+=rr[1];
   if( fabs(rr[0])>fabs(rr[1]) )
    b=fabs(rr[0]);else b=fabs(rr[1]);
    num++;
   }while(b>=eps);
   double x,y,Q,sx,sy,boty,QQ=0;
sx=50;
sy=5;
double sig=0;
   for (int i=0;i<sz;i++){
        Q=x0*pow(xy[0][i],y0);
 Add(xy[0][i],xy[1][i],Q);
 }
    expgraph(x0,y0);
}