Домашняя работа по информатике "Вложенные циклы FOR, решение задач"

Program kvur;
var
x1,x2,d,e,f, a, b, c: real;
begin
writeln('Решение квадратного уравнения(ax^2+bx+c=0).');
write('Введите a: ');
readln(a);
write('Введите b: ');
readln(b);
write('Введите c: ');
readln(c);
d:=sqr(b)-4*a*c;
writeln('вы ввели:');
writeln('a= ',a:6:2);
writeln('b= ',b:6:2);
writeln('c= ',c:6:2);
writeln('D=b*b-4*a*c=',d:6:2);
if (d>=0)
then begin
e:= (-b/(2*a));
f:= (sqrt(d)/(2*a));
x1:= (e+f);
x2:= (e-f);
writeln('D>= 0, значит корни действительные: ');
writeln('x1= ',x1:6:2);
writeln('x2= ',x2:6:2);
end;
if(d<0) then begin
e:= (-b/(2*a));
f:= (sqrt(-d)/(2*a));
x1:= (e+f);
x2:= (e-f);
writeln(' D<0, значит корни комплексные: ');
writeln('x1= ', e:6:2, ' + ', f:6:2, 'i');
writeln('x2= ', e:6:2, ' - ', f:6:2, 'i');
end;
end.

1. Вот так сейчас решают подобные задачи:

// PascalABC.NET 3.1, сборка 1174 от 22.02.2016
begin
  Writeln(Range(1,100).Select(i->sin(i)*cos(i)).Where(x->x<>0).Average)
end.

Тестовое решение:
-0.00136006072493969

2. А вот так учат писать это же школьные учителя:

// PascalABC.NET 3.1, сборка 1174 от 22.02.2016
var
  m:array[1..100] of real;
  i,k:integer;
  s:real;

begin
  s:=0;
  k:=0;
  for i:=1 to 100 do begin
    m[i]:=sin(i)*cos(i);
    if m[i]<>0 then begin
      s:=s+m[i];
      k:=k+1
      end
    end;
  Writeln(s/k)
end.

Тестовое решение:
-0.00136006072493969