Graph: построить семейство одинаковых кругов, центры которых расположены на тангенсоиде - Free Pascal

program tangensoid;
{$mode objfpc}{$H+}
{решаем задачу с [url]http://www.cyberforum.ru/free-pascal/thread1583591.html[/url]
 по зпдпнию: Построить семейство одинаковых кругов, центры которых
 расположены на тангенсоиде (Y = C*sin(X)/cos(X) на [–ПЂ/2: ПЂ/2]).
 Нарисовать оси координат.}
 
uses
  graph, math;
 
const
  R = 2;      // радиус "точки"
  DX = R*2;   //шаг по X
var
  gDriver, gMode: smallint;
  cx, cy: integer;  //координаты центра экрана
  c: single = 1.0;  //коэфициент для формулы
  i: integer;
  stepx: Single;    //шаг по X в радианах
  sf:    Single;
  f:    Single;     //значение тангенсойда
 
function tangensoid(const ax: single): single; inline;
begin
    Result:= sin(ax)/cos(ax);
end;
 
begin
  gDriver := detect;
  initgraph(gdriver, gmode, ' '); // Графический модуль
  cx := getmaxx div 2;
  cy := getmaxy div 2;
 
  //запрашивем коэф. C (0.1 - хорошее значение)
  write('Enter (C): ');
  ReadLn(c);
 
  //оси координат
  SetColor(white);
 
  Line(0, cy, GetMaxX, cy);
  Line(cx, 0, cx, GetMaxY);
 
 
  //рисуем "точки"
 
  //высчитываем шаг в радианах
  stepx:= pi / GetMaxX * DX;
  sf:=    GetMaxY / pi; // scale factor - кэф. чтобы привести значения ф-ции
                        //в диапазоне -pi/2 .. pi/2
                        // к диапазону разрешения экрана 0.. 768 (например)
 
  //начальное значение X
  SetColor(red);
  for i:= 0 to GetMaxX div DX do begin
    f:= sf * c * tangensoid(-(pi/2) + i*stepx);
    FillEllipse(i*DX, cy - trunc(f), R, R);
  end;
 
 
  Readln;
  CloseGraph;
end.