Численное интегрирование методом Гаусса (Квадратурная формула Гаусса)

//---------------------------------------------------------------------------
 
#include <vcl.h>
#pragma hdrstop
#include <math.h>
 
#include "Unit1.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm1 *Form1;
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
        : TForm(Owner)
{
}
//---------------------------------------------------------------------------
 
void __fastcall TForm1::Button1Click(TObject *Sender)
{
double a=StrToFloat(Edit1->Text);
double b=StrToFloat(Edit2->Text);
double *t;
t=new double [4];
double *A;
A=new double [4];
double *integ;
integ=new double [4];
 
t[0]=0.4442940564280;
t[1]=-0.4442940564280;
t[2]=1.026908012297201;
t[3]=-1.026908012297201;
A[0]=1.5008445;
A[1]=0.18250558;
A[2]=0.36003490;
A[3]=-0.04338504;
 
 
double integ1=0;
for(int i=0;i<4;i++)
{
 integ[i]=A[i]*sin((b+a)/2+(b-a)*t[i]/2);
}
 
for(int i=0;i<4;i++)
{
 integ1=integ1+integ[i];
}
integ1=(b-a)/2*integ1;
 Label1->Caption=FloatToStr(integ1);
}
//---------------------------------------------------------------------------