HalfFarey := proc (n) local a, b, c, d, k, s; if n<2 then RETURN([1]) fi; a:=0; b:=1; c:=1; d:=n; s:=NULL; do k := iquo(n+b, d); a, b, c, d := c, d, k*c-a, k*d-b; if b < 2*a then break fi; s := s, a/b od; [s] end:
A025555 := proc(n) local r; HalfFarey(n+1); subsop(nops(%) = NULL, %); mul(2*sin(Pi*r), r = %)^2 end: seq(round(evalf(A025555(i))), i=1..27); # Peter Luschny, Jun 09 2011
MATHEMATICA
nn=30; With[{trnos=Accumulate[Range[nn]]}, Table[LCM@@Take[trnos, n], {n, nn}]] (* Harvey P. Dale, Oct 21 2011 *)
f[x_] := x + 1; a[1] = f[1]; a[n_] := LCM[f[n], a[n - 1]]; Array[a, 30]/2 (* Robert G. Wilson v, Jan 04 2013 *)