为什么我的代码使用Ramanujan公式来计算pi的值,但为什么不能产生正确的值?

问题描述

package ValueOfPi;
import java.math.*;
import java.util.Scanner;

public class RamanujanAlgorithm {

    public static void main(String[] args) {
        
        double  constant = 0.0002885855652;
        
        Scanner sc = new Scanner(system.in);
        
        System.out.println("Enter limit of summation : ");
        
        int limit = sc.nextInt();
        
        BigDecimal pireci = new BigDecimal("0.00");
        
        for(int n = 0; n<=limit;n++) {
            
            BigDecimal summation = new BigDecimal((factorial(4*n)/Math.pow(factorial(n),4))* 
                    (26390*n + 1103)/Math.pow(396,4*n));
            
            
            pireci = pireci.add(summation);
        }
        
        pireci = pireci.multiply(BigDecimal.valueOf(constant)); 
        
    
        
        BigDecimal pi = BigDecimal.ONE.divide(pireci,pireci.scale(),RoundingMode.HALF_UP);
        
    
        
        System.out.println(pi);
    }
    
    public static long factorial(int num) {
        
        long factorial = 1;
        
        for(int i = 2; i<=num; i++) {
            
            factorial = factorial*i;
            
        }
        return factorial;
        
        
    }

}

公式来源-

我使用的拉马努詹公式是-

https://en.wikipedia.org/wiki/Pi#Rapidly_convergent_series 您看到的第一个公式是我使用的公式...

我的代码解释

在我的代码中,我首先将root(8)/ 9801的值存储在名为常量的变量中。 然后,我接受了用户输入的总和上限... 然后循环像Sigma一样工作,其中公式中sigma右边的所有内容都被计算并存储在BigDecimal pireci中,最后从循环中出来后,我将pireci与变量constant乘以

结果存储在pireci中。 因此,为了得到pi,我找到了pireci的倒数(因为公式给我的值是1 / pi)并打印结果 还有一种叫阶乘的方法,我已经计算出公式中的阶乘

产生的结果

在推杆限制= 1上我得到-

3.14159265383525216155408817653305805189369961082170142700979067898169315680169553

在推杆限制= 10上我得到-

3.141592653835251522017825683538367933441136808152187507749300560259859937277166891211043856677465587830817787411902252466648641215763228301502439410432665788310356378854402610761037687401200267944051297889143892295707008454487636735727641182639754208990826247667192290712119662965838581189023690921554268043065123435639930191440467906529034361442072184293368529887693375871414064191030918467799966026346716917502424050190813

而pi的真实值是

3.1415926535 8979323846 2643383279 5028841971 6939937510 5820974944 5923078164 0628620899 8628034825 3421170679 8214808651 3282306647 0938446095 ......

我的问题

您看到,在9位数字之后,我的代码获得的结果开始与实际值有所不同。为什么会这样??? 是由于公式还是由于我的代码??? 请帮我。 我还读到,精度随着求和极限的增加增加。我尝试使用更大的限制,但仍然出问题了。

在此先谢谢你

解决方法

暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!

如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@)