Составить алгоритм и программу для замены в заданном массиве

program P; 
var IshMas, NulID, PolMas, OtrMas:array [1..N] of integer;   
i, PolCount, OtrCount: integer;   
SredKvad, SredGeom: real; 
begin   
PolCount := 0;   
OtrCount := 0;   
for i := 1 to N do     
  if IshMas[i] > 0 then     
    begin       
     PolCount := PolCount + 1;       
     PolMas[PolCount] := IshMas[i];     
   end;     
  if IshMas[i] < 0 then     
    begin       
      OtrCount := OtrCount + 1;       
      OtrMas[OtrCount] := IshMas[i];     
    end;   
   SredKvad := 0;   
   SredGeom := 0;   
   if OtrCount > 0 then   
     begin     
        for i := 1 to OtrCount do       
           SredKvad := SredKvad + OtrMas[i]*OtrMas[i];     
           SredKvad := SredKvad/OtrCount;     
           SredKvad := sqrt(SredKvad);   
     end;   
     if PolCount > 0 then   
        begin     
           for i := 1 to PolCount do       
             SredGeom := SredGeom * PolMas[i]; 
            {    SredKvad := SredKvad/OtrCount;     
            SredKvad := sqrt(SredKvad);}   
        end;   
     for i := 1 to N do     
        if IshMas[i] = 0 then       
          if i mod 2 = 0 then 
            IshMas[i] := SredGeom       
          else IshMas[i] := SredKvad; 
  end;

//n-ая степень числа a 
function Stepen(a: real; n: integer): real; 
var j: integer; 
begin   
  Result := 1;   
  for i := 1 to n do     
    Result := Result * a; 
end;   
//корень n-ой степени числа a 
function Koren(a: real; n: integer): real; 
var j: integer;     x_1, x_2: real; 
begin   
x_1 := 0;   
x_2 := a/2;   
repeat     
  x_1 := x_2;     
  x_2 := ((n - 1)*x_1 + a/Stepen(x_1, n-1))/n   
until abs(x_1 - x_2) < 0.001 
end;