Программа "Построение поверхности"

Код программы:

unit Unit1;

interface

uses

 Windows, Messages, SysUtils, Classes, Graphics, Controls, Forms, Dialogs,
 StdCtrls, ExtCtrls, Spin;

type

 TForm1 = class(TForm)
   PaintBox1: TPaintBox;
   SpinEdit1: TSpinEdit;
   SpinEdit2: TSpinEdit;
   Label1: TLabel;
   Label2: TLabel;
   procedure PaintBox1Paint(Sender: TObject);
   procedure FormCreate(Sender: TObject);
   procedure SpinEdit1Change(Sender: TObject);
   procedure SpinEdit2Change(Sender: TObject);
   procedure PaintBox1MouseDown(Sender: TObject; Button: TMouseButton;
     Shift: TShiftState; X, Y: Integer);
   procedure PaintBox1MouseUp(Sender: TObject; Button: TMouseButton;
     Shift: TShiftState; X, Y: Integer);
   procedure PaintBox1MouseMove(Sender: TObject; Shift: TShiftState; X,
     Y: Integer);
 private
   { Private declarations }
 public
   { Public declarations }
 end;

var

 Form1: TForm1;

implementation {$R *.DFM} var O, A, B, C, OA, OB, OC, Mouse : TPoint;

 iP, R : integer;
 xmin, xmax, ymin, ymax, zmin, zmax : double;
 nx, ny : word;

// Функция инициализации структуры; function Init(x,y : Longint) : TPoint; var S : TPoint; begin

 S.x:=x;
 S.y:=y;
 Init:=S;

end; // Функция получения вектора по его началу и концу; function Vect(A,B : TPoint) : TPoint; var S : TPoint; begin

 S.x:=B.x-A.x;
 S.y:=B.y-A.y;
 Vect:=S;

end; // Функция проверяющая попадание курсора мыши в окресности заданной точки экрана function InCircle(A: TPoint; R : double): boolean; begin

 if sqr(Mouse.x-A.x)+sqr(Mouse.y-A.y) < R*R then InCircle:=True
 else InCircle:=False;

end; // Вычисление функции z function f(x,y : double) : double; begin

 f:=cos(x)*sin(y)/(1+y*y);

end; //Процедура созданиям формы procedure TForm1.FormCreate(Sender: TObject); begin

 xmin:=-2; xmax:=4;
 ymin:=-2; ymax:=4;
 zmin:=-2; zmax:=2;
 nx:=SpinEdit1.Value;
 ny:=SpinEdit2.Value;
 O:=Init(100, 200);
 A:=Init(20, 300);
 B:=Init(300, 200);
 C:=Init(100, 20);
 R:=5;

end; // Процедура рисования поверхности procedure TForm1.PaintBox1Paint(Sender: TObject);

 procedure MoveTo(A : TPoint);
 begin
   PaintBox1.Canvas.MoveTo(A.x, A.y);
 end;
 procedure LineTo(A : TPoint);
 begin
   PaintBox1.Canvas.LineTo(A.x, A.y);
 end;
 procedure Line(A, B : TPoint);
 begin
   PaintBox1.Canvas.MoveTo(A.x, A.y);
   PaintBox1.Canvas.LineTo(B.x, B.y);
 end;
 procedure Circle(A: TPoint; R : word);
 begin
   PaintBox1.Canvas.Ellipse(A.x-R,A.y-R,A.x+R,A.y+R);
 end;

var x, y, dx, dy, z, Qa, Qb, Qc : double;

 P : TPoint;
 kx, ky : integer;

begin

 OA:=Vect(O,A);
 OB:=Vect(O,B);
 OC:=Vect(O,C);
 //Вычисление шагов
 dx:=(xmax-xmin)/nx;
 dy:=(ymax-ymin)/ny;
 with PaintBox1 do begin
   //Фон изображения
   Canvas.Brush.Color:=clWhite;
   Canvas.Rectangle(0,0,Width, Height);
   //Рисование осей
   Circle(O,R);
   Line(O,A); Circle(A,R);
   Line(O,B); Circle(B,R);
   Line(O,C); Circle(C,R);
 end;
 //Рисование поверхности
 for kx:=0 to nx do begin
   x:=xmin+kx*dx;
   for ky:=0 to ny do begin
     y:=ymin+ky*dy;
     z:=f(x,y);
     Qa:=(x-xmin)/(xmax-xmin);
     Qb:=(y-ymin)/(ymax-ymin);
     Qc:=(z-zmin)/(zmax-zmin);
     P.x:=round(O.x+Qa*OA.x+Qb*OB.x+Qc*OC.x);
     P.y:=round(O.y+Qa*OA.y+Qb*OB.y+Qc*OC.y);
     if y=ymin then MoveTo(P) else LineTo(P);
   end;
 end;
 for ky:=0 to ny do begin
   y:=ymin+ky*dy;
   for kx:=0 to nx do begin
     x:=xmin+kx*dx;
     z:=f(x,y);
     Qa:=(x-xmin)/(xmax-xmin);
     Qb:=(y-ymin)/(ymax-ymin);
     Qc:=(z-zmin)/(zmax-zmin);
     P.x:=round(O.x+Qa*OA.x+Qb*OB.x+Qc*OC.x);
     P.y:=round(O.y+Qa*OA.y+Qb*OB.y+Qc*OC.y);
     if x=xmin then MoveTo(P) else LineTo(P);
   end;
 end;

end; //Количество точек по оси х procedure TForm1.SpinEdit1Change(Sender: TObject); begin

 nx:=SpinEdit1.Value;
 PaintBox1Paint(Sender);

end; //Количество точек по оси у procedure TForm1.SpinEdit2Change(Sender: TObject); begin

 ny:=SpinEdit2.Value;
 PaintBox1Paint(Sender);

end; //Определение номера выбранной точки procedure TForm1.PaintBox1MouseDown(Sender: TObject; Button: TMouseButton;

 Shift: TShiftState; X, Y: Integer);

begin

 iP:=-1;
 Mouse:=Init(X,Y);
 if InCircle(O, R) then iP:=0
 else if InCircle(A, R) then iP:=1
 else if InCircle(B, R) then iP:=2
 else if InCircle(C, R) then iP:=3;

end; //Процедура вычисления новых координат выбранной базисной точки и перерисовывания граффика procedure TForm1.PaintBox1MouseUp(Sender: TObject; Button: TMouseButton;

 Shift: TShiftState; X, Y: Integer);

begin

 case iP of
   0: O:=Init(X,Y);
   1: A:=Init(X,Y);
   2: B:=Init(X,Y);
   3: C:=Init(X,Y);
 end;
 PaintBox1Paint(Sender);

end; //Процедура обработки движения мыши procedure TForm1.PaintBox1MouseMove(Sender: TObject; Shift: TShiftState; X,Y: Integer); begin

 Mouse:=Init(X,Y);
 if InCircle(O,R) or InCircle(A,R) or InCircle(B,R) or InCircle(C,R)
 then Cursor:=crHandPoint
 else Cursor:=crDefault;

end; end.