blob: b743de976e328db95c33aba82749eec0f83a67ec (
plain)
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
|
> ||this tests that floating point overflow is handled correctly
Date: Mon, 16 Apr 90 17:15:03 CST
From: mccrosky@ishmael.usask.ca
Sorry for the delay in sending my bignum problem, I've been out of town.
> keith n = (shownum n) ++ " " ++ (shownum ((log (sumterm n))/n)) where
> sumterm n = sum (map prodterm [0..entier (n/2)]) where
> prodterm p = prod (map term [0..p-1]) where
> term j = x*(x-1)/((p-j)^2) where
> x = n-(2*j)
> prod = foldr (*) 1
> lim = lay (from 1) where from n = (keith n):(from (n*2))
******** This is the execution:
(We believe the results up to n=256 are correct).
Miranda lim
1 0.0
2 0.549306144334
4 0.736109744792
8 0.876176116589
16 0.966470019952
32 1.021895160467
64 1.054884461071
128 1.07405223475
256 1.084981322415
512
should trap floating point overflow here
|
