|
|
option base 0
DIM f(1000)
DIM g(1000)
LET f(0)=0
LET f(1)=6
INPUT n
FOR i=1 TO n
LET f(i+1)=(34*i^3+51*i^2+27*i+5)/(i+1)^3*f(i)-i^3/(i+1)^3*f(i-1)
NEXT i
PRINT "f(n)＝";f(n)
LET g(0)=1
LET g(1)=5
FOR i=1 TO n
LET g(i+1)=(34*i^3+51*i^2+27*i+5)/(i+1)^3*g(i)-i^3/(i+1)^3*g(i-1)
NEXT i
PRINT "g(n)＝";g(n)
PRINT "f(n)/g(n)＝";f(n)/g(n)
PRINT "比較：ζ(3)=1+1/2^3+1/3^3+1/4^3+･････"
PRINT " ＝1.2020569031595942853997381615114499907649･････"
END
１０００桁モードで見て下さい。
|
|