让a1成为第一个词，r是每个词乘以得到下一个词并且n是术语数的常量，几何级数为： ai = a1*r**(i-1)，pn n条款和sn n条款总和的乘积。
我有计算公式，但Fortran 95（Plato2）不承认我需要的精度。 （例如：我不能得到-1.234E + 00567890）。
如何“充分”使用这个“巨大”数字的双精度？
答案 0 :(得分：1)
如此高的数字（您的示例-1.234E + 00567890）对于Fortran标准提供的任何内在数字类型来说太大了。它们也大于物理和工程应用中使用的数字。例如，我的gfortran支持以下类型：
huge(1.0_real32) 3.40282347E+38
huge(1.0_real64) 1.7976931348623157E+308
huge(1.0_real128) 1.18973149535723176508575932662800702E+4932
据我所知，没有可用的Fortran编译器具有更大的内部浮点类型。
出于特殊目的，与您的一样，需要专门的库。本网站不适用于软件推荐，因此我不会推荐任何特定的网站。在http://crd-legacy.lbl.gov/~dhbailey/mpdist/查看其中一些列表，当然还有更多（GNU科学图书馆将有一些我肯定的）。