График функции

{$R *.dfm}
procedure TForm1.FormCreate(Sender: TObject);
begin
 align:=alClient;
end;
 
procedure TForm1.FormPaint(Sender: TObject);
const al = pi/4;
 // x^2/a^2 + y^2/b^2 - z^2/c^2 = 1
 // однополостный гиперболоид
 // для заданных x, y точки на поверхности может не существовать
 // - функция возвращает false
 // или существует 2 точки - z1, z2
 function f(a,b,c:integer; x,y:real; var z1, z2 :real):boolean;
 var m:real;
 begin
  z1 := 0; z2 := 0;
  m := ((x*x)/(a*a)+(y*y)/(b*b) -1) * (c*c);
  result := (m >= 0) and (m < 1e37);
  if not result then exit;
  z1 := sqrt(m);
  z2 := -z1;
 end;
var a, b, c:intege;
    x1, y1, x0, y0, m:integer;
    x, y, dx, dy, z1, z2:real;
begin
// Задание параметров поверхности
a := 10; //strtoint(inputbox('окно ввода','введи a',' '));
b := 10; //strtoint(inputbox('окно ввода','введи b',' '));
c := 15; //strtoint(inputbox('окно ввода','введи c',' '));
// центр системы координат - в центре формы
x0 := width div 2;
y0 := height div 2;
// Рисуем систему координат
canvas.MoveTo(x0-10, y0); canvas.lineTo(width-50, y0); canvas.textout(width-40,y0,'X');
canvas.MoveTo(x0, y0+10); canvas.lineTo(x0, 0+50);     canvas.textout(x0, 0+40, 'Z');
canvas.MoveTo(x0+trunc(10*sin(al)), trunc(y0-10*cos(al)));
x1 := x0-trunc(y0*sin(al));
y1 := y0+trunc(y0*cos(al));
canvas.lineTo(x1, y1);
canvas.textout(x1, y1, 'Y');
m := 2;   // масштабный коэффициент
dx := 1;
dy := 1;
x := -100;
while x <= 100 do
 begin
  y := -100;
  while y <= 100 do
   begin
    if f(a, b, c, x, y, z1, z2) then
     begin
      x1 := trunc(x0 + m*x - m*y*cos(al));
      y1 := trunc(y0 + m*y*cos(al)- m*z1);
      canvas.Pixels[x1,y1] := clpurple;
      x1 := trunc(x0 + m*x - m*y*cos(al));
      y1 := trunc(y0 + m*y*cos(al)- m*z2);
      canvas.Pixels[x1,y1] := clpurple;
     end;
    y := y + dy;
   end;
  x := x + dx;
 end;
end;

function f(a,b,c:integer; x,y:real; var z1, z2 :real):boolean;
 var m:real;
 begin
  z1 := 0; z2 := 0;
  m := ((x*x)/(a*a)+(y*y)/(b*b) -1) * (-c*c);//вот здесь добавить минус
  result := (m >= 0) and (m < 1e37);
  if not result then exit;
  z1 := sqrt(m);
  z2 := -z1;
 end;

// Задание параметров поверхности
a := 100; //strtoint(inputbox('окно ввода','введи a',' '));
b := 60; //strtoint(inputbox('окно ввода','введи b',' '));
c := 80; //strtoint(inputbox('окно ввода','введи c',' '));

m := ((x*x)/(a*a)+(y*y)/(b*b) -1) * (c*c);