GeezerVisonThe World Wide Web, restored to 800 pixels of glory. |
| Home|About|Gaming cabinet|Archives|How it works|Starter sites |
GeezerVison is showing github.com — Java/Others/GuassLegendre.java at master · feiteng/Java Expand file tree/ Copy pathGuassLegendre.javaMore file actions 46 lines (36 loc) · 1017 Bytes / Copy pathGuassLegendre.javaFile metadata and controls46 lines (36 loc) · 1017 Bytes 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 package Others; import java.lang.Math; /** * Guass Legendre Algorithm * ref https://en.wikipedia.org/wiki/Gauss–Legendre_algorithm * * @author AKS1996 */ public class 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, b = Math.pow(2, -0.5), t = 0.25, 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; } } |
| You are visitor 01402520 — GeezerVison does not store the pages you look up. Best viewed with Netscape Navigator 4.7 or Internet Explorer 6 at 800×600. No cookies required. |