Построить график функции ln(x)

uses graph;
function F(x:real):real;
begin
F:=ln(x);
end;
var xn,xk,x,y,mx,my,dx:real;
    x0,y0,i,px,py:integer;
    s:string;
begin
x0:=0;
initgraph(x0,y0,'');
x0:=30;{правый край экрана}
y0:=getmaxY div 2;{центр}
xn:=0.1;xk:=20;{интервал по Х}
mx:=(getmaxX-2*x0)/(xk-xn);{масштаб по Х}
my:=(y0-20)/(F(xk));{по У}
setbkcolor(15);{цвет фона}
setcolor(8);{цвет линий}
line(0,y0,getmaxX,y0);{оси}
line(x0,0,X0,getmaxY);
for i:=1 to round(xk)+1 do{максимальное количество засечек в одну сторону}
  begin
    line(x0-3,y0+round(i*my),x0+3,y0+round(i*my));{засечки на оси У}
    line(x0+round(i*mx),y0-3,x0+round(i*mx),Y0+3); {засечки на оси Х}
    line(x0-round(i*mx),y0-3,x0-round(i*mx),Y0+3);
    str(i,s);
    {подпись оси У}
    outtextXY(x0-15,y0+round(i*my),'-'+s);{соответственно засечкам}
    outtextXY(x0-15,y0-round(i*my),s);
    {подпись оси Х}
    outtextXY(x0+round(i*mx),y0+10,s);
    outtextXY(x0-round(i*mx),y0+10,'-'+s);
  end;
{центр}
outtextXY(x0+5,y0+10,'0');
{подписи концов осей}
outtextXY(getmaxX-10,y0-10,'X');
outtextXY(x0-10,10, 'Y');
setcolor(8);
{график}
x:=xn;
dx:=0.001;
while x<=xk do
   begin
     x:=x+dx; {наращиваем х}
     px:=x0+round(x*mx); {изменение по оси Х}
     y:=F(x); {график в реальных координатах}
     py:=y0-round(y*my);{в экранных}
     putpixel(px,py,8);
   end;
readln
end.