Построение графика в декартовой системе координат

//---------------------------------------------------------------------------
 
#include <vcl.h>
#pragma hdrstop
 
#include "Unit1.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm1 *Form1;
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
                         : TForm(Owner)
{
}
#include "math.h"
float f(float x){
return 2*sin(x)*exp(x/5);
}
//---------------------------------------------------------------------------
 
void __fastcall TForm1::OnPaint(TObject *Sender)
{
float x1,x2,y1,y2; //границы изменения аргумента и значения функции
float x,y; //аргумент функции и значение функции в точке х
float dx; //приращение аргумента
int l,b; //левый нижний угол области вывода
int w,h; //ширина и высота области вывода
int x0,y0; //начало координат
float mx, my; //масштаб по осям Х и У
// область вывода графика
l=10;
b=Form1->ClientHeight-20;
h=Form1->ClientHeight-40;
w=Form1->Width-20;
x1=0; //нижняя граница диапазона аргумента
x2=25; //верхняя граница диапазона аргумента
dx=0.05;
y1=f(x1); //минимально значение функции
y2=f(x2); //максимально значение функции
 
do
{
y=f(x);
if(y<y1) y1=y;
if(y>y2) y2=y;
x+=dx;
}
while(x<=x2);
my=(float)h/(abs)(y2-y1); //масташб по оси У
my=w/(abs)(x2-x1); //масштаб по оси Х
 
//оси
x0=1+abs(x2-x1);
y0=b-abs(y1*my);
Canvas->MoveTo(x0,b); Canvas->LineTo(x0,b-h);
Canvas->MoveTo(1,y0); Canvas->LineTo(1+w,y0);
x=x1;
//Построение графика
do
{
y=f(x);
Canvas->Pixels[x0+x*mx][y0+y*my]=clRed;
x+=dx;
}
while(x<=x2);
}