codepad
[
create a new paste
]
login
|
about
Language:
C
C++
D
Haskell
Lua
OCaml
PHP
Perl
Plain Text
Python
Ruby
Scheme
Tcl
import std.cstream; import std.math; void main() { int N = 10; real[] a = new real[N], b = new real[N], t = new real[N], pi = new real[N]; a[0] = 1.0; b[0] = 1.0 / sqrt(2.0); t[0] = 1.0 / 4.0; for(int i=0; i<N-1; i++) { a[i+1] = (a[i] + b[i]) / 2.0; b[i+1] = sqrt(a[i] * b[i]); t[i+1] = t[i] - pow(2.0, cast(real)i) * pow(a[i] - a[i+1], 2.0); } for(int i=0; i<N; i++) { dout.writefln("%40.37f", pow(a[i] + b[i], 2.0) / (4.0 * t[i])); } dout.writefln("%40.37f", real.epsilon); }
Private
[
?
]
Run code
Submit