Построение графика функций

uses GraphABC;
 
var 
 
  R001: real;
  Hs, Hr, Lg, Lr, Ls, teta, r,  Le, A001, Ap, F,  av, ah, kkp, akp, yr, kh, kv, B, B1, Re, y, q, dzeta, hi,  fi, p, v: double;
  i,x0,y0,x1,x2,y1,k: integer;
  mx,my,dx,x:real;
   w, n:integer;
    s:string;
 
begin
 
  
B:=0:
A001:=18.031;
 
x1:=-3;//начало и конец графика в реальных значениях
x2:=15;
 
 
begin
  x0:=windowwidth div 2;//начало координат по оси Х
  y0:=windowheight div 2; //по оси Y
  mx:=(x0-30)/x2;//масштаб по Х
  my:=(y0-50)/x2;//масштаб по Y
  n:=12;//количество делений по осям в 1 сторону
  line(20,y0,windowwidth-20,y0); //ось Х
  line(x0,20,x0,windowheight-20); //ось У
  for i:=1 to n do
    begin
      line(x0-3,y0-round(i*my/2),x0+3,y0-round(i*my/2));//засечки на оси У
      line(x0-3,y0+round(i*my/2),x0+3,y0+round(i*my/2));
      line(x0+round(i*mx/2),y0+3,x0+round(i*mx/2),y0-3); //засечки на оси Х
      line(x0-round(i*mx/2),y0+3,x0-round(i*mx/2),y0-3);
      str(i/2:0:1,s);
      //подпись оси У
      textout(x0-20,y0-round(i*my/2),s);{соответственно засечкам}
      textout(x0-25,y0+round(i*my/2),'-'+s);
      //подпись оси Х
      textout(x0+round(i*mx/2),y0+10,s);
      textout(x0-round(i*mx/2),y0+10,'-'+s);
    end;
      //центр
      textout(x0+5,y0+10,'0');
      //подписи концов осей
      textout(windowwidth-30,y0-20,'X');
    textout(x0+10,10, 'Y');
    //график
      x:=0.01;
      while x<=0.1 do
        begin
          Ap:=(A001*exp(-(0.655+0.033*ln(x)-0.045*ln(A001)-B*(1-x)*sin(teta*pi/180))*ln(x/0.01)));
           setpixel(x0+round(x*mx),y0-Ap*my,clRed);
           x:=x+0.01; //наращиваем х
        end;
 
end;
 
 
 
end.