График функции Гаусса

uses graphABC;
function F(d,x:real):real;
begin
F:=d*exp(-x*x);
end;
var a,b,mx,my,x,dx,y1,y2,y3,y4:real;
    x0,y0,px,py,i:integer;
    s:string;
begin
a:=-5;
b:=5;
x0:=windowwidth div 2;{начало координат}
y0:=windowheight div 2;
{коэффициенты по осям для перевода реальных координат в экранные}
mx:=(windowwidth-60)/(b-a);{длина в пикселах по оси Х/реальная длина}
my:=(y0-80)/F(2,0);{длина оси по У/(макс}
line(0,y0,windowwidth,y0); {ось Х}
line(x0,windowwidth,x0,0); {ось У}
for i:=1 to round(b)+1 do   {максимальное количество засечек в одну сторону}
   begin
    line(x0-3,y0-round(i*my),x0+3,y0-round(i*my));{засечки на оси У}
    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);
    {подпись оси У}
    textout(x0-20,y0-round(i*my),s);{соответственно засечкам}
    textout(x0-25,y0+round(i*my),'-'+s);
    {подпись оси Х}
    str(i,s);
    textout(x0+round(i*mx),y0+10,s);
    textout(x0-round(i*mx),y0+10,'-'+s);
  end;
{центр}
textout(x0+5,y0+10,'0');
{подписи концов осей}
textout(windowwidth-10,y0-10,'X');
textout(x0+10,10, 'Y');
dx:=0.001;
{график}
x:=a;
while x<=b do
 begin
  px:=x0+round(x*mx); {изменение по оси Х}
  y1:=F(-2,x); {график в реальных координатах}
  py:=y0-round(y1*my);{в экранных}
  setpixel(px,py,clRed);
  y2:=F(-0.5,x); {график в реальных координатах}
  py:=y0-round(y2*my);{в экранных}
  setpixel(px,py,clGreen);
  y3:=F(0.5,x); {график в реальных координатах}
  py:=y0-round(y3*my);{в экранных}
  setpixel(px,py,clBlue);
  y4:=F(2,x); {график в реальных координатах}
  py:=y0-round(y4*my);{в экранных}
  setpixel(px,py,clBrown);
  x:=x+dx; {наращиваем х}
 end;
setfontcolor(clBrown);
//settextstyle(0,0,1);
textout(95,10,'График функции Y=D*exp(-x*x)');
textout(145,30,'интервал [-5; 5]');
setfontcolor(clRed);
textout(185,50,'A=-2');
setfontcolor(clGreen);
textout(185,70,'A=-0.5');
setfontcolor(clBlue);
textout(185,90,'A=0.5');
setfontcolor(clYellow);
textout(185,110,'A=2');
readln
end.