%I
%S 0,1,2,3,4,5,6,7,8,9,11,22,33,44,101,111,121,131,141,151,161,171,191,
%T 202,212,222,232,242,252,272,303,313,323,333,353,404,414,434,515,787,
%U 868,949,1001,1111,1221,1331,2002,2112,2222,2992,3003,3113,3883,4004
%N Palindromes for which the sum of the digits is also a palindrome.
%e 4774 is in the sequence because (1) it is a palindrome and (2) the sum of its digits 4 + 7 + 7 + 4 = 22 is also a palindrome.
%p rev := proc (n) local nn; nn := convert(n, base, 10): add(nn[j]*10^(nops(nn)j), j = 1 .. nops(nn)) end proc: sd := proc (n) options operator, arrow: add(convert(n, base, 10)[j], j = 1 .. nops(convert(n, base, 10))) end proc: a := proc (n) if rev(n) = n and rev(sd(n)) = sd(n) then n else end if end proc: seq(a(n), n = 0 .. 5000); # _Emeric Deutsch_, Jan 18 2009
%t d[n_]:=IntegerDigits[n]; sod[n_]:=Total[d[n]]; palQ[n_]:=Reverse[x=d[n]]==x; t={}; Do[If[palQ[n] && palQ[sod[n]],AppendTo[t,n]],{n,0,4005}]; t (* _Jayanta Basu_, May 15 2013 *)
%Y Cf. A002113, A082207.
%K base,nonn
%O 1,3
%A _Amarnath Murthy_, Apr 10 2003
%E More terms from Luc Stevens (lms022(AT)yahoo.com), Apr 16 2006; corrected Apr 18 2006; corrected Apr 20 2006
%E Edited by _N. J. A. Sloane_, Jan 14 2009 at the suggestion of _Ian Kent_
