%I #27 Sep 08 2022 08:45:15
%S 0,3,28,117,336,775,1548,2793,4672,7371,11100,16093,22608,30927,41356,
%T 54225,69888,88723,111132,137541,168400,204183,245388,292537,346176,
%U 406875,475228,551853,637392,732511,837900,954273,1082368,1222947
%N a(n) = n^4 + n^3 + n^2.
%C a(n) are the numbers m such that: j^2 = j + m + sqrt(j*m) with corresponding numbers j given by A002061(n+1), and with sqrt(j*m) = A027444(n) = n* A002061(n+1). - _Richard R. Forberg_, Sep 03 2013.
%H Vincenzo Librandi, <a href="/A100019/b100019.txt">Table of n, a(n) for n = 0..1000</a>
%F From _Indranil Ghosh_, Apr 15 2017: (Start)
%F G.f.: -x(3 + 13x + 7x^2 + x^3)/(x - 1)^5
%F E.g.f.: exp(x)*x*(3 + 11x + 7x^2 + x^3)
%F (End)
%p A100019:=n->n^4+n^3+n^2: seq(A100019(n), n=0..50); # _Wesley Ivan Hurt_, Apr 14 2017
%t f[n_]:=n^4+n^3+n^2;Array[f,50,0] (* _Vladimir Joseph Stephan Orlovsky_, Apr 12 2011 *)
%o (Magma) [n^4+n^3+n^2: n in [0..50]]; // _Vincenzo Librandi_, Jun 09 2011
%o (PARI) a(n) = n^4 + n^3 + n^2 \\ _Indranil Ghosh_, Apr 15 2017
%Y Cf. A000583, A002523, A091940, A071253, A059826, A027445, A053699.
%K nonn,easy
%O 0,2
%A Douglas Winston (douglas.winston(AT)srupc.com), Nov 19 2004