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

uses graphABC;
function f(x:real):real;
begin
f:=cos(2*x)+cos(0.5*x);
end;
function y(x:real):real;
begin
if(f(x)<0)and(x>=pi)and(x<=1.5*pi)then y:=1
else y:=f(x);
end;
var xn,xk,x,dx,mx,my:real;
    x0,y0,i,n:integer;
    s1,s2:string;
begin
setwindowsize(800,500);
centerwindow;
//интервал
xn:=0;
xk:=2*pi;
//шаг для функции f(x)
dx:=pi/18;
//таблица значений функции f(x)
setfontsize(10);
textout(60,0,'Таблица значений');
textout(10,20,'аргумента и функции f(x) шаг=2pi/36');
setpencolor(clBlack);
line(10,40,230,40);
textout(30,42,'x');
textout(90,42,'f(x)');
textout(145,42,'x');
textout(200,42,'f(x)');
line(10,60,230,60);
x:=xn;
n:=60;//шаг по вертикали в таблице
while x<=xk/2 do
 begin
  str(x:0:2,s1);
  str(f(x):0:3,s2);
  textout(20,n,s1);
  textout(80,n,s2);
  line(10,n,120,n);
  inc(n,20);
  x:=x+dx;
 end;
line(10,40,10,n);
line(65,40,65,n);
line(120,40,120,n);
line(10,n,120,n);
x:=pi+dx;
n:=60;
while x<=xk do
 begin
  str(x:0:2,s1);
  str(f(x):0:3,s2);
  textout(130,n,s1);
  textout(185,n,s2);
  line(120,n,230,n);
  inc(n,20);
  x:=x+pi/18;
 end;
line(120,n,230,n);
line(10,40,10,n);
line(175,40,175,n);
line(230,40,230,n);
line(10,n,120,n);
//построение графика
//начало координат
x0:=250+(windowwidth-250) div 10;
y0:=windowheight div 2;
mx:=(windowwidth-x0-20)/2/pi;
my:=(y0-20)/2;//максимальное значение функции
line(x0+round(xn*mx),y0,x0+round(xk*mx),y0);
textout(x0+round(xk*mx),y0-20,'X');
line(x0,0,x0,2*y0);
textout(x0+5,0,'Y');
textout(x0+5,y0+10,'0');
//засечки и подписи по оси Х
for i:=1 to 4 do
 begin
  line(x0+round(i*pi*mx/2),y0-3,x0+round(i*pi*mx/2),y0+3);
  line(x0-round(i*pi*mx/2),y0-3,x0-round(i*pi*mx/2),y0+3);
  case i of
  1:s1:='0.5pi';
  2:s1:='pi';
  3:s1:='1.5pi';
  4:s1:='2pi';
  end;
  textout(x0+round(i*pi*mx/2),y0+10,s1);
 end;
//засечки и подписи по оси Y
for i:=1 to 20 do
if i mod 2=0 then
 begin
  line(x0-3,y0-round(i*my/10),x0+3,y0-round(i*my/10));
  line(x0-3,y0+round(i*my/10),x0+3,y0+round(i*my/10));
  textout(x0-25,y0-round(i*my/10)-5,floattostr(i/10));
  textout(x0-30,y0+round(i*my/10)-5,floattostr(-i/10));
 end;
//график функции y(f(x))
setpencolor(clBlue);
x:=xn;
while x<=xk do
 begin
  if x=xn then moveto(x0+round(x*mx),y0-round(y(x)*my))
  else lineto(x0+round(x*mx),y0-round(y(x)*my));
  x:=x+0.01;
 end;
textout(x0+150,20,'График функции y=f(x)');
end.