1 | test_product = ->
|
2 | run_test [
|
3 |
|
4 | # compute pi using Viete's formula ------------
|
5 |
|
6 | # note that this is not an array, this
|
7 | # defines a recursive function
|
8 | "a(n)=test(n=0,0,sqrt(2+a(n-1)))",
|
9 | "",
|
10 |
|
11 | # not very efficient because evaluation of
|
12 | # a(n) is not memoized, so there
|
13 | # is quadratic cost as n increases.
|
14 | "float(2*product(2/a(k),k,1,9))",
|
15 | "3.14159",
|
16 |
|
17 | "a = quote(a)",
|
18 | "",
|
19 |
|
20 | # Wallis' product
|
21 | "2*product(float(4*k^2/(4*k^2-1)),k,1,100)",
|
22 | "3.13379",
|
23 |
|
24 |
|
25 | ]
|