Расчет значений тригонометрической функции на заданном интервале

function F(x: real): real;
begin
  F := 3 * sin(x) + cos(x);
end;
 
function y(x: real): real;
begin
  if (x >= 0) and (x <= 2 * pi) then
    y := f(x)
  else
  if f(x) < 0 then
    if (x >= pi) and (x <= 1.5 * pi) then
      y := -1;
end;
 
var
  x: real;
  g: text;
 
begin
  assign(g, 'data.txt');
  rewrite(g);
  x := 0;
  while x <= 2 * pi do
  begin
    writeln(g, x:12:3, f(x):12:3);
    x := x + 2 * pi / 36;
  end;
  close(g);
end.

uses GraphABC;
 
const dx = 2*pi/36;     // шаг вычислений
      // граница вычислений по оси Х
      xmin = 0;
      xmax = 2*pi;
 
var x1, x2, y1, y2 : real;  // переменные для построения графика
 
//------ Функция вычисления F(x)-------------//
 function Solution(x: real): real;
  begin
   result:=3*sin(x) + cos(x);
  end;
 
//---- Записываем вычисления в файл ------//
 procedure WriteDataToFile;
  var
   f: text;
   x: real;
  begin
   x:=xmin;
   Assign(f,'c:\d.txt');
   Rewrite(f);
   writeln(f,'----------------------');
   writeln(f,'|    x     |    y      |');
   writeln(f,'----------------------');
   while x <= xmax do
    begin
     writeln(f,'|  ', x:4:2,'  | ',Solution(x):6:2, '  |');
     x:=x+dx;
    end;
   writeln(f,'----------------------');
   Close(f);
  end;
 
 
BEGIN
      WriteDataToFile;
      
      setwindowsize(220, 280);         //устанавливаем размер графического окна
      line(10, 0, 10, 280);            //ось ордината
      line(0, 140, 220, 140);          //ось абсцисс
 
// Отмечаем точками интервал [п, 1.5*п]
      Circle(110,140,3);               //Пи
      Circle(160,140,3);               //1.5*Пи
 
{-----------------------------------------------------------------------------
Для построения графика воспользуемся ф-цией:
Line(X1,Y1,X2,Y2 : integer) - рисует линию от точки (X1,Y1) до точки (X2,Y2);
-----------------------------------------------------------------------------}
      x1:=xmin;
      x2:=xmin;
      while x2 <= xmax do
      begin
           y1:=-Solution(x1);
           x2:=x1+dx;
           while ((x2>=pi) and (x2 <= 1.5*pi)) do
           begin
                 if y1 >= 0 then y2:=1;
                 Line(10+round(x1*100/pi),140+round(y1*100/3),10+round(x2*100/pi),140+round(y2*100/3));
                 x1:=x2;
                 x2:=x1+dx;
           end;
           y2:=-Solution(x2);
           Line(10+round(x1*100/pi),140+round(y1*100/3),10+round(x2*100/pi),140+round(y2*100/3));
           x1:=x2;
      end;
END.

xmin = 0;

function Solution(x: real): real;
  begin
   result:=3*sin(x) + cos(x);
  end;

x1:=xmin;

y1:=-Solution(x1)