MathAction SandBoxCS224

• FrontPage
  • SandBox
    • SandBoxCS224 <-- You are here.

(1) -> artanh:(Float) -> Float

artanh(x) == 
 y := x 
 k := 1
 for k in 1..precision() repeat
 z := x^(2*k+1) / (2 * k + 1)
 y := z + y
 return y

ln2:() -> Float

ln2() == artanh(0.5) + artanh(1.0 / 7)

gamma:(NonNegativeInteger) -> Float

gamma(e) ==
 precision(2 ^ e + 100)
 n := 2.0 ^ e
 eps := 1.0 / n
 A := - e * ln2()
 B := 1.0
 U := A 
 V := 1.0
 k := 1
 repeat
 B := B * n^2 / k^2
 A := (A * n^2 / k + B) / k
 if (A < eps) /\ (B < eps) then
 return U / V
 U := U + A
 V := V + B
 k := k + 1

fractions:(NonNegativeInteger, Integer) -> List(Integer)

fractions(e, n) ==
 y := gamma(e)
 l : List(Integer) := []
 for i in 0..n repeat
 f := floor(y)
 l := cons(f, l)
 y := y - f
 if y > 0 then 
 y := 1 / y
 reverse(l)

fractions(12, 1000)

Compiling function artanh with type Float -> Float

Compiling function ln2 with type () -> Float

Compiling function gamma with type NonNegativeInteger -> Float

Compiling function fractions with type (NonNegativeInteger, Integer)
 -> List(Integer)