Вычислить s=1/33+1/53+1/73+..+1/(2n+1)3 в паскале

Var
  S, tmp: real;
  n, i: integer;

begin
      write('Введите n:');
      readln(n);
      for i := 1 to n do
      begin
        tmp := 1 / ((2 * i + 1) * 3);
        S := S + tmp;
      end;
      writeln(S:8:10);
      readln;
end.

// PascalABC.NET 3.1, сборка 1218 от 12.04.2016
uses GraphABC;

procedure OutLightGreen(n:integer);
begin
  var OldColor:=FontColor;
  SetFontColor(clLightGreen);
  Write(n:5);
  SetFontColor(OldColor)
end;

begin
  var a:=MatrixRandom(10,14,-180,320);
  var m:=ReadInteger('n=');
  Writeln;
  SetFontColor(clWhite);
  Font.Style:=fsBold; Font.Name:='Courier';
  Brush.Color:=clBlack; Window.Clear(clBlack);
  var k:=0;
  for var i:=0 to 9 do begin
    for var j:=0 to 13 do
      if a[i,j]>m then begin
        Inc(k); OutLightGreen(a[i,j]) end
      else Write(a[i,j]:5);
    Writeln
    end;
  SetFontColor(clRed);
  Writeln('Количество элементов равно ',k)
end.

Так как язык не указан, приведу пример на SWI-Prolog.

Код:

read_int(Int) :- read(Int), integer(Int).split_int_by_numbers(0, []) :- !.split_int_by_numbers(N, [Number|Ints]) :- Number is mod(N, 10), RestN is div(N, 10), split_int_by_numbers(RestN, Ints).test_to_div(_, []).test_to_div(N, [Number|Ints]) :- mod(N, Number) =:= 0, test_to_div(N, Ints). test(Int) :- split_int_by_numbers(Int, Numbers), test_to_div(Int, Numbers), write(Int), write(" - Yes!"), nl.test(Int) :- write(Int), write(" - No!"), nl.?- read_int(Int), test(Int).