mirror of
https://github.com/TheAlgorithms/Java.git
synced 2025-07-25 13:34:54 +08:00
49 lines
1.1 KiB
Java
49 lines
1.1 KiB
Java
package com.thealgorithms.others;
|
||
|
||
/**
|
||
* Guass Legendre Algorithm ref
|
||
* https://en.wikipedia.org/wiki/Gauss–Legendre_algorithm
|
||
*
|
||
* @author AKS1996
|
||
*/
|
||
public final class GuassLegendre {
|
||
private GuassLegendre() {
|
||
}
|
||
|
||
public static void main(String[] args) {
|
||
for (int i = 1; i <= 3; ++i) {
|
||
System.out.println(pi(i));
|
||
}
|
||
}
|
||
|
||
static double pi(int l) {
|
||
/*
|
||
* l: No of loops to run
|
||
*/
|
||
|
||
double a = 1;
|
||
double b = Math.pow(2, -0.5);
|
||
double t = 0.25;
|
||
double p = 1;
|
||
for (int i = 0; i < l; ++i) {
|
||
double[] temp = update(a, b, t, p);
|
||
a = temp[0];
|
||
b = temp[1];
|
||
t = temp[2];
|
||
p = temp[3];
|
||
}
|
||
|
||
return Math.pow(a + b, 2) / (4 * t);
|
||
}
|
||
|
||
static double[] update(double a, double b, double t, double p) {
|
||
double[] values = new double[4];
|
||
values[0] = (a + b) / 2;
|
||
values[1] = Math.sqrt(a * b);
|
||
values[2] = t - p * Math.pow(a - values[0], 2);
|
||
values[3] = 2 * p;
|
||
|
||
return values;
|
||
}
|
||
}
|