|
%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
|
