не могу разобраться с ошибкой

[Smirnov]

посмотрите пожалуйста мою программу...тут лажа какая то, не считает нифига
сама программа

uses crt,graph,mod1;
var
grDriver : Integer;
grMode : Integer;
Errcode,l : Integer;
r: real;
Lx1,Lx2,Ly1,Ly2:integer;
kx, {kofix}
ky, {kofiy}
min,
max:real;
otstx,
otsty,
v,u,o,R1:integer;
BEGIN
Rp[1]:=0.25; Rp [2]:=0.13; Rp [3]:=0.06; {proektivniy probeg}
dr [1]:=0.07; dr [2]:=0.05; dr [3]:=0.02; {standartnoe otklonenie}

ClrScr;
writeln ('vvedite Xmax');
Readln (Xmax);
writeln ('vvedite Xmin');
Readln (Xmin);
writeln ('vvedite m');
Readln (m);
zap;

grDriver:= Detect;
InitGraph(grDriver, grMode, 'D:\BP\bgi');
ErrCode:= GraphResult;
if ErrCode = grOk then
Lx1:=80; {Xo}
Lx2:=480; {Xkon}
Ly1:=40; {Yo}
Ly2:=440; {Ykon}
otstx:=80; {otstup po x}
otsty:=39; {otstyp po y}
v:=11;
r:=40; {shag setki}

R1:=Q*3;
min:=A[2,1];
max:=A[2,1];
for u:=2 to 4 do
for o:=1 to R1 do
begin



if A[u,o]<min then min:=A[u,o];
if A[u,o]>max then max:=A[u,o];
end;

kx:=(Xmax-Xmin)/(Lx2-Lx1);
ky:=(Ly2-Ly1)/(max-min);
writeln (kx);
writeln (ky);

for l:=1 to v do
begin
Line(80, 40+round(r*(l-1)), 480, 40+round(r*(l-1)));
Line(80+round(r*(l-1)), 40, 80+round(r*(l-1)) , 440);
end;
readln;

end.
это модуль
unit mod1;
interface
const
Pi=3.1415926536;
N=1e14; {doza}
var
xmax, {diaoapozon izmeneniy}
xmin,
{proektivniy probeg}
{ctandartnoe otklonenie}
m, {shag izmenenia x}
Cp: {pikovya concentr}
real;
i,
Q: {chislo povtorov}
integer;
A:array[1..4,1..100] of real;
Rp,dr:array[1..3] of real;

procedure zap;
function prov (Xmax,Xmin:real):boolean;
function rasch (x:real;i:integer):real;

implementation

function prov(Xmax,Xmin:real):boolean;
begin
prov:=True;
If ((Xmin < 0) or (Xmax > 1) or (Xmin>=Xmax)) then
begin
Write('vveden neverniy diapazon');
prov:=False;
Exit;
end;
end;

function rasch (x:real;i:integer):real;
begin
Cp:=N/(sqrt(2*Pi)*dr[i]);
rasch:=Cp*(exp(-sqr(x-Rp[i])/(2*sqr(dr[i]))));
end;

procedure zap;
var x:real;
i,j:integer;
begin
If prov(Xmax, Xmin) then
begin
x:=Xmin;
Q:=Round((Xmax-Xmin)/m);
Writeln(' x,mkm C(B) C(P) C(As) ');
for i:=1 to Q+1 do
begin
A[1,i]:=x;
for j:=1 to 3 do
A[j+1,i]:= rasch(x,j);
writeln (x:3,' ', A[2,i]:3,' ', A[3,i]:3,' ', A[4,i]:3);
x:=x+m;

end;
end;

readln
end;
end.