A181184 - OEIS

%I #10 Mar 31 2012 10:24:25

%S 2,5,3,7,5,6,1,3,3,5,2,4,0,1,0,1,6,1,0,2,6,2,1,4,3,5,5,1,1,1,4,1,3,6,

%T 7,7,3,6,7,7,2,6,0,2,2,0,5,7,0,6,6,3,0,4,2,8,4,1,5,0,4,1,7,0,4,0,2,6,

%U 5,6,0,6,4,6,7,2,7,0,6,5,5,1,1,2,1,2,3,3,6,2,1,2,4,7,2,8,6,1,3,0,5,0,2,6,1,4

%N | (n-th digit of Pi) - (n-th digit of phi (golden ratio)) |.

%p std:=105:   Digits:=2*std:

%p num1:=Pi:   num2:=(1+sqrt(5))/2:

%p A181184:=proc(n) local a, b, c, d, i, am;

%p a:=evalf(num1-trunc(num1),2*std):

%p b:=evalf(num2-trunc(num2),2*std):

%p am:=[ abs(trunc(num1)-trunc(num2)) ]:

%p for i from 0 by 1 to n do

%p    c:=trunc(10*a);   d:=trunc(10*b);

%p    a:=frac(evalf(10*a, 2*std));   b:=frac(evalf(10*b, 2*std));

%p    am:=[op(am),abs(c-d)]:

%p od: end: A181184(std);

%t Abs[RealDigits[Pi, 10, 150][[1]] - RealDigits[GoldenRatio, 10, 150][[1]]]

%Y Cf. A000796, A001622, A073223.

%K nonn,base

%O 0,1

%A _Jani Melik_, Jan 26 2011