Распределение Максвелла при исходном двумерном законе Гаусса определяется формулой

uses crt,graphABC;
const n=10;
      r:array[1..n] of real=(0,0.2,0.4,0.6,0.8,1,1.5,2,2.5,3);
var f:array[1..n] of real;
    i,x0,y0:integer;
    a,mx,my,max:real;
    s:string;
begin
a:=0.8;
max:=0;
for i:=1 to n do
 begin
  f[i]:=pi*r[i]*exp(-pi*sqr(r[i])/(4*a*a))/(2*a*a);
  if f[i]>max then max:=f[i];
  writeln('r=',r[i]:0:1,'  f=',f[i]:0:4);
 end;
readln;
clrscr;
hidecursor;
x0:=60; //начало координат
y0:=windowheight-60;
mx:=(windowwidth-2*x0)/r[n];  //масштаб по Х
my:=(y0-60)/max;  //по У
setpencolor(clBlack);  //оси
line(x0,y0,windowwidth-x0,y0);
line(x0,y0,X0,30 );
//засечкт и подписи на осях
for i:=1 to n do
  begin
    line(x0+round(i*mx),y0-3,x0+round(i*mx),Y0+3);
    str(i,s);
    textout(x0+round(i*mx)-15,y0+10,s);
    line(x0+3,y0-round(i*my/10),x0-3,y0-round(i*my/10));
    str(i/10:0:1,s);
    textout(x0-30,y0-round(i*my/10),s);
  end;
textout(x0+5,y0+10,'0');
textout(windowwidth-60,y0-15,'X');
textout(x0-10,10, 'Y');
//график
moveto(x0+round(r[1]*mx),y0-round(f[1]*my));
for i:=1 to n do
 begin
  setpencolor(clred);
  lineto(x0+round(r[i]*mx),y0-round(f[i]*my)); //линии
  setpencolor(clgreen);
  circle(x0+round(r[i]*mx),y0-round(f[i]*my),3); //точки
 end;
end.