1 3940506b 2005-01-13 devnull scale = 50
2 3940506b 2005-01-13 devnull define e(x) {
3 3940506b 2005-01-13 devnull auto a, b, c, d, e, g, w, y, t, r
9 3940506b 2005-01-13 devnull scale = t + .434*x + 1
12 3940506b 2005-01-13 devnull if(x<0) {
17 3940506b 2005-01-13 devnull while(x>2) {
27 3940506b 2005-01-13 devnull for(a=1; 1; a++) {
29 3940506b 2005-01-13 devnull c = c*a+b
32 3940506b 2005-01-13 devnull if(g == e) {
34 3940506b 2005-01-13 devnull while(y--) {
37 3940506b 2005-01-13 devnull scale = t
38 3940506b 2005-01-13 devnull if(w==1) {
39 3940506b 2005-01-13 devnull ibase = r
40 3940506b 2005-01-13 devnull return 1/g
42 3940506b 2005-01-13 devnull ibase = r
43 3940506b 2005-01-13 devnull return g/1
49 3940506b 2005-01-13 devnull define l(x) {
50 3940506b 2005-01-13 devnull auto a, b, c, d, e, f, g, u, s, t, r, z
52 3940506b 2005-01-13 devnull r = ibase
53 3940506b 2005-01-13 devnull ibase = A
54 3940506b 2005-01-13 devnull if(x <= 0) {
55 3940506b 2005-01-13 devnull z = 1-10^scale
56 3940506b 2005-01-13 devnull ibase = r
59 3940506b 2005-01-13 devnull t = scale
62 3940506b 2005-01-13 devnull scale += scale(x) - length(x) + 1
63 3940506b 2005-01-13 devnull s = scale
64 3940506b 2005-01-13 devnull while(x > 2) {
65 3940506b 2005-01-13 devnull s += (length(x)-scale(x))/2 + 1
66 3940506b 2005-01-13 devnull if(s>0) {
67 3940506b 2005-01-13 devnull scale = s
69 3940506b 2005-01-13 devnull x = sqrt(x)
72 3940506b 2005-01-13 devnull while(x < .5) {
73 3940506b 2005-01-13 devnull s += (length(x)-scale(x))/2 + 1
74 3940506b 2005-01-13 devnull if(s>0) {
75 3940506b 2005-01-13 devnull scale = s
77 3940506b 2005-01-13 devnull x = sqrt(x)
81 3940506b 2005-01-13 devnull scale = t + length(f) - scale(f) + 1
82 3940506b 2005-01-13 devnull u = (x-1)/(x+1)
84 3940506b 2005-01-13 devnull scale += 1.1*length(t) - 1.1*scale(t)
90 3940506b 2005-01-13 devnull for(a=3; 1; a=a+2){
92 3940506b 2005-01-13 devnull c = c*a + d*b
95 3940506b 2005-01-13 devnull if(g==e) {
96 3940506b 2005-01-13 devnull scale = t
97 3940506b 2005-01-13 devnull ibase = r
98 3940506b 2005-01-13 devnull return u*c/d
104 3940506b 2005-01-13 devnull define s(x) {
105 3940506b 2005-01-13 devnull auto a, b, c, s, t, y, p, n, i, r
107 3940506b 2005-01-13 devnull r = ibase
108 3940506b 2005-01-13 devnull ibase = A
109 3940506b 2005-01-13 devnull t = scale
110 3940506b 2005-01-13 devnull y = x/.7853
111 3940506b 2005-01-13 devnull s = t + length(y) - scale(y)
112 3940506b 2005-01-13 devnull if(s<t) {
115 3940506b 2005-01-13 devnull scale = s
116 3940506b 2005-01-13 devnull p = a(1)
118 3940506b 2005-01-13 devnull scale = 0
119 3940506b 2005-01-13 devnull if(x>=0) {
120 3940506b 2005-01-13 devnull n = (x/(2*p)+1)/2
122 3940506b 2005-01-13 devnull if(x<0) {
123 3940506b 2005-01-13 devnull n = (x/(2*p)-1)/2
125 3940506b 2005-01-13 devnull x -= 4*n*p
126 3940506b 2005-01-13 devnull if(n%2 != 0) {
130 3940506b 2005-01-13 devnull scale = t + length(1.2*t) - scale(1.2*t)
131 3940506b 2005-01-13 devnull y = -x*x
135 3940506b 2005-01-13 devnull for(i=3; 1; i+=2) {
137 3940506b 2005-01-13 devnull b *= i*(i-1)
139 3940506b 2005-01-13 devnull if(c==0){
140 3940506b 2005-01-13 devnull scale = t
141 3940506b 2005-01-13 devnull ibase = r
142 3940506b 2005-01-13 devnull return s/1
148 3940506b 2005-01-13 devnull define c(x) {
149 3940506b 2005-01-13 devnull auto t, r
151 3940506b 2005-01-13 devnull r = ibase
152 3940506b 2005-01-13 devnull ibase = A
153 3940506b 2005-01-13 devnull t = scale
154 3940506b 2005-01-13 devnull scale = scale+1
155 3940506b 2005-01-13 devnull x = s(x + 2*a(1))
156 3940506b 2005-01-13 devnull scale = t
157 3940506b 2005-01-13 devnull ibase = r
158 3940506b 2005-01-13 devnull return x/1
161 3940506b 2005-01-13 devnull define a(x) {
162 3940506b 2005-01-13 devnull auto a, b, c, d, e, f, g, s, t, r, z
164 3940506b 2005-01-13 devnull r = ibase
165 3940506b 2005-01-13 devnull ibase = A
166 3940506b 2005-01-13 devnull if(x==0) {
167 3940506b 2005-01-13 devnull return 0
169 3940506b 2005-01-13 devnull if(x==1) {
170 3940506b 2005-01-13 devnull z = .7853981633974483096156608458198757210492923498437764/1
171 3940506b 2005-01-13 devnull ibase = r
172 3940506b 2005-01-13 devnull if(scale<52) {
173 3940506b 2005-01-13 devnull return z
176 3940506b 2005-01-13 devnull t = scale
178 3940506b 2005-01-13 devnull while(x > .5) {
180 3940506b 2005-01-13 devnull x = -(1 - sqrt(1.+x*x))/x
183 3940506b 2005-01-13 devnull while(x < -.5) {
185 3940506b 2005-01-13 devnull x = -(1 - sqrt(1.+x*x))/x
188 3940506b 2005-01-13 devnull s = -x*x
193 3940506b 2005-01-13 devnull for(a=3; 1; a+=2) {
195 3940506b 2005-01-13 devnull c = c*a + d*b
198 3940506b 2005-01-13 devnull if(g==e) {
199 3940506b 2005-01-13 devnull scale = t
200 3940506b 2005-01-13 devnull ibase = r
201 3940506b 2005-01-13 devnull return x*c/d
207 3940506b 2005-01-13 devnull define j(n,x) {
208 3940506b 2005-01-13 devnull auto a,b,c,d,e,g,i,s,k,t,r
210 3940506b 2005-01-13 devnull r = ibase
211 3940506b 2005-01-13 devnull ibase = A
213 3940506b 2005-01-13 devnull t = scale
214 3940506b 2005-01-13 devnull k = 1.36*x + 1.16*t - n
215 3940506b 2005-01-13 devnull k = length(k) - scale(k)
216 3940506b 2005-01-13 devnull if(k>0) {
217 3940506b 2005-01-13 devnull scale += k
220 3940506b 2005-01-13 devnull s = -x*x/4
221 3940506b 2005-01-13 devnull if(n<0) {
227 3940506b 2005-01-13 devnull for(i=1; i<=n; i++) {
229 3940506b 2005-01-13 devnull c *= 2*i
234 3940506b 2005-01-13 devnull for(i=1; 1; i++) {
236 3940506b 2005-01-13 devnull b = b*i*(n+i) + a
237 3940506b 2005-01-13 devnull c *= i*(n+i)
239 3940506b 2005-01-13 devnull if(g==e) {
240 3940506b 2005-01-13 devnull scale = t
241 3940506b 2005-01-13 devnull ibase = r
242 3940506b 2005-01-13 devnull return g/1