Решение уравнения

const m=3;
      n=5;
      p=4;
var x:array[1..m,1..n] of real;
    y:array[1..p,1..p] of real;
    a,a1,b,b1,x1,x2:real;
    i,j,k:integer;
begin
randomize;
a:=0;
writeln('Матрица Х');
for i:=1 to m do
 begin
  a1:=1;
  for j:=1 to n do
   begin
    x[i,j]:=-9+18*random;
    a1:=a1*x[i,j];
    write(x[i,j]:6:2);
   end;
  a:=a+a1;
  writeln;
 end;
writeln('Матрица Y');
b:=0;
for k:=1 to p do
 begin
  b1:=1;
  for j:=1 to p do
   begin
    y[i,j]:=-9+18*random;
    b1:=b1*y[k,j];
    write(y[i,j]:6:2);
   end;
  b:=b+b1;
  writeln;
 end;
writeln('a=',a:0:2,'  b=',b:0:2);
if -b/a<0 then write('Уравнение не имеет решений')
else
 begin
  x1:=-sqrt(-b/a);
  x2:=sqrt(-b/a);
  write('x1=',x1:0:2,'  x2=',x2:0:2);
 end;
 
end.

type mas=array[1..100,1..100] of real;
var x,y:mas;
    a,a1,b,b1,x1:real;
    i,j,k,n,m,t,q:integer;
    
procedure mas1(var x:mas; n,m:integer);
var i,j:integer;
begin
for i:=1 to n do
 begin
  for j:=1 to m do
readln(x[i,j]);
end;
for i:=1 to n do 
    begin
      for j:=1 to m do 
        write(x[i,j]:4);
        writeln();
    end;
    writeln;
end;
 
 
    
function ai( x:mas; n, m:integer):real;
var a,a1:real;
    i,j:integer;
begin
a:=0;
for i:=1 to n do
 begin
  a1:=1;
  for j:=1 to m do
   begin
    a1:=a1*x[i,j];
   end;
  a:=a+a1;
 end;
 ai:=a;
end;
 
 
    begin
q:=4;
t:=4;
n:=3;
m:=5;
writeln('Введите  массивы');
mas1(x,n,m);
mas1(y,q,t);
x1:= -ai(y,q,t)/ai(x,n,m);
write('x^2=', x1:2:2);
end.