Интерполяционный многочлен Лагранжа

Program grafik;
Uses graph, crt;
const
x1=-10;
x2=10;
var
  gm, gd, i:integer;
  xa, ya, px, py:integer;
  x, y, mx, my, dx:real;
  s,s1:string;
begin
  gd:=detect;
  initgraph(gd,gm,'');
  xa:=GetMaxX div 2;
  ya:=GetMaxY div 2;
  mx:=2*xa/(x2-x1);
  my:=2*ya/10;
  Line(0,ya,GetMaxX,ya);
  Line(xa,GetmaxY,xa,0);
  For i:=1 to 10 do
   begin
     line(xa-3,ya-round(i*my),xa+3,ya-round(i*my));
     line(xa-3,ya+round(i*my),xa+3,ya+round(i*my));
     line(xa+round(i*mx),ya-3,xa+round(i*mx),ya+3);
     line(xa-round(i*mx),ya-3,xa-round(i*mx),ya+3);
     str(i,s);
     outtextXY(xa-25,ya-round(i*my),s);
     outtextXY(xa-25,ya+round(i*my),'-'+s);
     str(i,s);
     outtextXY(xa+round(i*mx),ya+10,s);
     outtextXY(xa-round(i*mx),ya+10,'-'+s);
   end;
  outtextXY(xa+5,ya+10,'0');
  outtextXY(GetMaxX-10,ya-10,'X');
  outtextXY(xa-10,10,'Y');
  dx:=0.001;
  x:=0;
  while x<=x2 do
    begin
     px:=xa+round(x*mx);
     if abs(x)>0.001 then
       begin
         y:=ln(1+sqrt(x))-cos(x);
         py:=ya-round(y*my);
         putpixel(px,py,4);
       end;
     x:=x+dx;
    end;
  Readln;
  CloseGraph;
End.