Построить График функции Arccosx в паскале

uses graphABC;
const e=0.001;//точность вычисления корней
function f(x:real):real;
begin
if x=0 then f:=pi/2
else if (x>0)and(x<=1) then f:=arctan(sqrt(1-x*x)/x)
else if(x>=-1)and(x<0) then f:=pi+arctan(sqrt(1-x*x)/x);
end;
 
var xn,xk,x,dx,mx,my,x1,x2,x3,y:real;
    x0,y0,n,i:integer;
    s:string;
begin
setwindowsize(500,300);
centerwindow;
//начало координат
x0:=windowwidth div 2;
y0:=windowheight-50;
//оси
line(0,y0,2*x0,y0);
textout(2*x0-20,y0-20,'X');
line(x0,0,x0,y0);
textout(x0+5,0,'Y');
textout(x0+5,y0+10,'0');
xn:=-1;
xk:=1;
mx:=(x0-20);
my:=(y0-20)/pi;
//подписи по оси Х
for i:=1 to 10 do
 begin
  line(x0+round(i*mx/10),y0-3,x0+round(i*mx/10),y0+3);
  line(x0-round(i*mx/10),y0-3,x0-round(i*mx/10),y0+3);
  textout(x0+round(i*mx/10),y0+10,floattostr(i/10));
  textout(x0-round(i*mx/10),y0+10,floattostr(-i/10));
 end;
//подписи по оси Y
line(x0-5,y0-round(pi*my/2),x0+5,y0-round(pi*my/2));
line(x0-5,y0-round(pi*my),x0+5,y0-round(pi*my));
textout(x0-25,y0-round(pi*my/2)-5,'pi/2');
textout(x0-30,y0-round(pi*my)-5,'pi');
setpencolor(clBlue);
x:=xn;
dx:=0.01;
while x<=xk do
 begin
  if x=xn then moveto(x0+round(x*mx),y0-round(f(x)*my))
  else lineto(x0+round(x*mx),y0-round(f(x)*my));
  x:=x+dx;
 end;
end.