editing
approved
editing
approved
P:=proc(n) local i, j, k, w; for i from 1 by 1 to n do k:=i; w:=i-2; while w>0 do k:=k*w; w:=w-2; od; j:=i+1; w:=i-1; while w>0 do j:=j*w; w:=w-2; od; print(j-k); od; end: P(100); # Paolo P. Lava, Jun 14 2007
approved
editing
(MAGMAMagma) DoubleFactorial:=func< n | &*[n..2 by -2] >; [DoubleFactorial((n-1))-DoubleFactorial(n-2): n in [3..30]]; // Vincenzo Librandi, Aug 08 2017
reviewed
approved
proposed
reviewed
editing
proposed
(n-1)*a(n+2) = a(n+1) + n^2*a(n). - Robert Israel, Aug 08 2017
f:= gfun:-rectoproc({(-n+1)*a(2+n)+a(1+n)+n^2*a(n), a(2)=0, a(3)=1}, a(n), remember):
map(f, [$3..100]); # Robert Israel, Aug 08 2017
proposed
editing
editing
proposed
1, 1, 5, 7, 33, 57, 279, 561, 2895, 6555, 35685, 89055, 509985, 1381905, 8294895, 24137505, 151335135, 468934515, 3061162125, 10033419375, 68000295825, 234484536825, 1645756410375, 5943863027025, 43105900812975, 162446292283275, 1214871076343925, 4761954230608575
Table[(n - 1)!! - (n - 2)!!, {n, 3, 30}] (* Vincenzo Librandi, Aug 08 2017 *)
(MAGMA) DoubleFactorial:=func< n | &*[n..2 by -2] >; [DoubleFactorial((n-1))-DoubleFactorial(n-2): n in [3..30]]; // Vincenzo Librandi, Aug 08 2017
reviewed
editing
proposed
reviewed