问题描述
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 (将#修改为@)