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

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

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.