При вводе отрицательных чисел, программа глючит(пишет: Программа выполняеца), а при вводе любых других чисел работает нормально

const n=3;
type mas=array[1..n,1..n] of real;
function f(x,y:real;k:integer):real;
begin
 f:=1/exp(1/2)-exp(1/2)*sin((x-y)/k);
end;
function g(x,y:real;k:integer):real;
begin
 if x<y then g:=1/sqrt(2)-sin((x*x-y)/k*k)
        else g:=1/ln(3)-cos((x*x+y)/k*k)
end;
function p(a:mas):real;
var i,j:integer;
k:real;
begin
for i:=1 to n do
for j:=1 to n do
begin
 if (abs(a[i,j])<1)and(a[i,j]>k) then
        k:=a[i,j];
        p:=a[i,j];
end;
end;
function r(a:mas):real;
var i,j,k:integer;
     q:real;
begin
 for i:=1 to n do
  for j:=1 to n do
   if i>j then
      begin
        q:=q+a[i,j];
        inc(k);
      end;
      r:=q/k;
end;
procedure k12(a:mas; var k1,k2:integer);
var i,j:integer;
 
begin
for k1:=1 to n do
for k2:=1 to n do
for i:=1 to n do
for j:=1 to n do
if (abs(a[i,j])<a[k1,k2]) then
      begin
        k1:=i;
        k2:=j;
end;
end;
function w(k2:integer;p,r:real):real;
var k:integer;
    s:real;
begin
s:=0;
 for k:=1 to k2 do
  s:=s+f(p,r,k)-g(p,r,k)/2;
w:=s;
end;
function z(k1:integer;p,r:real):real;
var k:integer;
    d:real;
begin
d:=1;
 for k:=1 to k1 do
  d:=d*((f(r,p,k)-sqrt(2))/g(r,p,k));
z:=d;
end;
function a(z,w:real):real;
begin
 a:=z/sqrt(z*z+w*w);
end;
function b(w,z:real):real;
begin
 b:=w/sqrt(z*z+w*w);
end;
function h(k:integer;t,a,b:real):real;
begin
 h:=cos(exp(1/3*ln(a*a*t+b*b)))/sqr(1+k);
end;
var o:mas;
    i,j,k,k1,k2:integer;
    s,t,v:real;
begin
s:=0;
for i:=1 to n do
 for j:=1 to n do
  readln(o[i,j]);
  k12(o,k1,k2);
t:=0.2;
while t<0.5 do
begin
k:=0;
repeat
k:=k+1;
v:=h(k,t,a(z(k1,p(o),r(o)),w(k2,p(o),r(o))),b(z(k1,p(o),r(o)),w(k2,p(o),r(o))));
s:=s+v;
until v<10e-6;
writeln('T=',t:4:2,' S=',s:2:1);
t:=t+0.04;
end;
readln;
end.