О.О.П / lab16 / Unit1

//---------------------------------------------------------------------------

#include <vcl.h>
#pragma hdrstop

#include "Unit1.h"
#include "Math.h"
double fun(double);
double Mn_Lagr(double*,double,int);
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm1 *Form1;
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
        : TForm(Owner)
{
}
//---------------------------------------------------------------------------
void __fastcall TForm1::Button1Click(TObject *Sender)
{
double x,h,h1, a, b, *mas_x, *mas_y_t;
int i,n,m;
a = StrToFloat(Edit3->Text);
b = StrToFloat(Edit4->Text);
m = StrToInt(Edit1->Text);
n = StrToInt(Edit2->Text);
h = (b-a)/(m-1);
h1 = (b-a)/(n-1);
mas_x = new double[m+1];
mas_y_t = new double[n+1];
for(x=a, i=0; i<m; i++)
{
  mas_x[i] = x;
  x+=h;
}
Memo1->Lines->Add("---- Мнoгoчлeн Лaгpaнжa ---");
Memo1->Lines->Add("Пoлyчили " + IntToStr(n) + " знaчeний:");
for(x=a, i=0; i<n; i++, x+=h1)
{
  mas_y_t[i] = Mn_Lagr(mas_x,x,m);
  Memo1->Lines->Add(" x = "+FloatToStrF(x,ffFixed,8,2)+"  f*(x) = "+FloatToStrF(mas_y_t[i],ffFixed,8,4));
}
Chart1->Series[0]->Clear();
Chart1->Series[1]->Clear();
for(x=a-0.1; x<b+0.1; x+=0.01)
  Chart1->Series[0]->AddXY(x,fun(x));
for(x=a,i=0; i<n; i++,x+=h1)
  Chart1->Series[1]->AddXY(x,mas_y_t[i]);
delete []mas_x;
delete []mas_y_t;
}
//---------------------------------------------------------------------------
double fun(double x) {
    	return pow(x,2)-10*pow(sin(x),2);
}
//---------------------------- Мнoгoчлeн Лaгpaнжa -------------------------------------------
double Mn_Lagr(double *x, double xt, int kol)  {
int i, k;
double e, p=0;
for(k=0; k<kol; k++)
{
  e=1.;
  for (i=0;i<kol;i++)
    if (i!=k)
      e *= ((xt-x[i])/(x[k]-x[i]));
  p += e*fun(x[k]);
}
return p;
}