Fortran Coder

关于算常数e的10000位数 大神来呀

查看数: 8805 | 评论数: 8 | 收藏 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2014-9-26 04:56

正文摘要:

最近看见这个代码,不清楚这个e是怎么算出来的,求解答或者共同探讨 代码见下: (数组的初附值都为0) 以及e的泰勒级数 e=\sum_{n=0}^{\infty}\frac{1}{n!}=\frac{1}{0!}+\frac{1}{1!}+\frac{1}{2!}+\frac{1}{3! ...

回复

wjy8888840 发表于 2014-9-26 22:40:37
楚香饭 发表于 2014-9-26 15:40
你给出的只是部分代码,我无法帮你查看是什么问题。

我也没有能力根据部分代码,推测出所有代码。 ...

sorrysorry  这部分的已经想出来了!谢谢!!
楚香饭 发表于 2014-9-26 15:40:55
你给出的只是部分代码,我无法帮你查看是什么问题。

我也没有能力根据部分代码,推测出所有代码。
wjy8888840 发表于 2014-9-26 15:25:05

  是少了什么,,求具体点,, 试了很多了,,
楚香饭 发表于 2014-9-26 15:19:34
需完整代码
wjy8888840 发表于 2014-9-26 15:11:19
楚香饭 发表于 2014-9-26 12:01
具体要看 Silverfrost Ftn95 的限定了,通常来说,数组的长度限定在 4294967296 左右 ...

compile只有write部分一直跳出来missing expression

什么情况 怎么解决呢
[Fortran] 纯文本查看 复制代码
open(60,file="e.txt")
write(60,'(1a)'),"e=2."
write(60,'(8x,77i1)')(ns(i),i=2,77)
write(60,'(80i1)')(ns(i),i=78,nk)
write(*,'(a)'),'The ten digital primes in e:'
write(60,'(a)'),'The ten digital primes in e:'
do i=2,nk-10
C i=2
tmp=0
if(ns(i).ne.0) then
do j=i,i+9
tmp=tmp*10+ns(j)
end do
nst=sqrt(tmp)
do it=2,nst
if (mod(tmp,it).eq.0) then 
goto 13
end if 
end do
13 if (it.gt.nst) then
write(*,'(i5,f15.1)'),i,tmp
write(60,'(i5,f15.1)'),i,tmp
end if
end if
end do
楚香饭 发表于 2014-9-26 12:01:18
具体要看 Silverfrost Ftn95 的限定了,通常来说,数组的长度限定在 4294967296 左右
wjy8888840 发表于 2014-9-26 11:22:37

非常有用!
如果是用silverfrost  plato的话,对array出来的长度有没有限制呢?
楚香饭 发表于 2014-9-26 09:38:00
这个应该是用
e=\sum_{n=0}^{\infty}\frac{1}{n!}=\frac{1}{0!}+\frac{1}{1!}+\frac{1}{2!}+\frac{1}{3!}+\frac{1}{4!}+...

公式的。

通过除法,竖式加法等基本运算获得。(回忆一下小学的除法和加法)

有篇文章你可以参考,《4种计算自然常数e的方法及精度比较》
http://www.doc88.com/p-990521706366.html

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-11-23 01:20

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表