Решить квадратное уравнение

procedure kvur (a,b,c: real; var x1,x2: real);
var d: real;
begin
d:=sqr(b)-4*a*c;
if d>0 then 
begin
x1:=(-b+sqrt(d))/(2*a);
x2:=(-b-sqrt(d))/(2*a);
end
else if d=0 then 
begin
x1:=-b/(2*a);
x2:=x1;
end;
end;

program proc;
var
d:array[1..10,1..10] of integer;
x:array[1..8,1..8] of integer;
y:array[1..6,1..6] of integer;
a1,b1,c1,xp,xp2:real;
i,j,a,n,b,m,c,l:integer;
procedure kvur (a,b,c: real; var x1,x2: real);
var d: real;
begin
d:=sqr(b)-4*a*c;
if d>0 then
begin
x1:=(-b+sqrt(d))/(2*a);
x2:=(-b-sqrt(d))/(2*a);
end
else if d=0 then
begin
x1:=-b/(2*a);
x2:=x1;
end;
write ('x1=',x1:3:3,' x2=',x2:3:3);
end;
begin
randomize;
for i:=1 to 10 do
begin
writeln;
for j:=1 to 10 do
begin
d[i,j]:=random(21)-10;
write (d[i,j],' ');
if d[i,j]>0 then
begin
a:=a+d[i,j];
n:=n+1;
end;
end;
end;
writeln ('a=',a,' n=',n);
for i:=1 to 8 do
begin
writeln;
for j:=1 to 8 do
begin
x[i,j]:=random(21)-10;
write (x[i,j],' ');
if x[i,j]>0 then
begin
b:=b+x[i,j];
m:=m+1;
end;
end;
end;
writeln ('b=',b,' m=',m);
for i:=1 to 6 do
begin
writeln;
for j:=1 to 6 do
begin
y[i,j]:=random(21)-10;
write (y[i,j],' ');
if y[i,j]>0 then
begin
c:=c+y[i,j];
l:=l+1;
end;
end;
end;
writeln ('c=',c,' l=',l);
a1:=a/n;
b1:=b/m;
c1:=c/l;
kvur(a1,b1,c1,xp,xp2);
readln;
end.