version 1.1, 2003/09/25 19:34:22 |
version 1.2, 2007/01/29 11:02:53 |
|
|
*/ |
*/ |
|
|
scale = 20 |
scale = 20 |
define e(x){ |
define e(x) { |
auto a, b, c, d, e, g, t, w, y |
auto a, b, c, d, e, g, t, w, y |
|
|
t = scale |
t = scale |
scale = t + .434*x + 1 |
scale = t + .434*x + 1 |
|
|
w = 0 |
w = 0 |
if(x<0){ |
if (x < 0) { |
x = -x |
x = -x |
w = 1 |
w = 1 |
} |
} |
y = 0 |
y = 0 |
while(x>2){ |
while (x > 2) { |
x = x/2 |
x = x/2 |
y = y + 1 |
y = y + 1 |
} |
} |
|
|
a=1 |
a = 1 |
b=1 |
b = 1 |
c=b |
c = b |
d=1 |
d = 1 |
e=1 |
e = 1 |
for(a=1;1==1;a++){ |
for (a = 1; 1 == 1; a++) { |
b=b*x |
b = b*x |
c=c*a+b |
c = c*a + b |
d=d*a |
d = d*a |
g = c/d |
g = c/d |
if(g == e){ |
if (g == e) { |
g = g/1 |
g = g/1 |
while(y--){ |
while (y--) { |
g = g*g |
g = g*g |
} |
} |
scale = t |
scale = t |
if(w==1) return(1/g) |
if (w == 1) return(1/g) |
return(g/1) |
return(g/1) |
} |
} |
e=g |
e = g |
} |
} |
} |
} |
|
|
define l(x){ |
define l(x) { |
auto a, b, c, d, e, f, g, u, s, t |
auto a, b, c, d, e, f, g, u, s, t |
if(x <=0) return(1-10^scale) |
if (x <= 0) return (1 - 10^scale) |
t = scale |
t = scale |
|
|
f=1 |
f = 1 |
scale = scale + scale(x) - length(x) + 1 |
scale = scale + scale(x) - length(x) + 1 |
s=scale |
s = scale |
while(x > 2){ |
while (x > 2) { |
s = s + (length(x)-scale(x))/2 + 1 |
s = s + (length(x) - scale(x))/2 + 1 |
if(s>0) scale = s |
if (s > 0) scale = s |
x = sqrt(x) |
x = sqrt(x) |
f=f*2 |
f = f*2 |
} |
} |
while(x < .5){ |
while (x < .5) { |
s = s + (length(x)-scale(x))/2 + 1 |
s = s + (length(x) - scale(x))/2 + 1 |
if(s>0) scale = s |
if (s > 0) scale = s |
x = sqrt(x) |
x = sqrt(x) |
f=f*2 |
f = f*2 |
} |
} |
|
|
scale = t + length(f) - scale(f) + 1 |
scale = t + length(f) - scale(f) + 1 |
u = (x-1)/(x+1) |
u = (x - 1)/(x + 1) |
|
|
scale = scale + 1.1*length(t) - 1.1*scale(t) |
scale = scale + 1.1*length(t) - 1.1*scale(t) |
s = u*u |
s = u*u |
|
|
c = b |
c = b |
d = 1 |
d = 1 |
e = 1 |
e = 1 |
for(a=3;1==1;a=a+2){ |
for (a = 3; 1 == 1 ; a = a + 2) { |
b=b*s |
b = b*s |
c=c*a+d*b |
c = c*a + d*b |
d=d*a |
d = d*a |
g=c/d |
g = c/d |
if(g==e){ |
if (g == e) { |
scale = t |
scale = t |
return(u*c/d) |
return (u*c/d) |
} |
} |
e=g |
e = g |
} |
} |
} |
} |
|
|
define s(x){ |
define s(x) { |
auto a, b, c, s, t, y, p, n, i |
auto a, b, c, s, t, y, p, n, i |
t = scale |
t = scale |
y = x/.7853 |
y = x/.7853 |
s = t + length(y) - scale(y) |
s = t + length(y) - scale(y) |
if(s<t) s=t |
if (s < t) s = t |
scale = s |
scale = s |
p = a(1) |
p = a(1) |
|
|
scale = 0 |
scale = 0 |
if(x>=0) n = (x/(2*p)+1)/2 |
if (x >= 0) n = (x/(2*p) + 1)/2 |
if(x<0) n = (x/(2*p)-1)/2 |
if (x < 0) n = (x/(2*p) - 1)/2 |
x = x - 4*n*p |
x = x - 4*n*p |
if(n%2!=0) x = -x |
if (n % 2 != 0) x = -x |
|
|
scale = t + length(1.2*t) - scale(1.2*t) |
scale = t + length(1.2*t) - scale(1.2*t) |
y = -x*x |
y = -x*x |
a = x |
a = x |
b = 1 |
b = 1 |
s = x |
s = x |
for(i=3; 1==1; i=i+2){ |
for (i =3 ; 1 == 1; i = i + 2) { |
a = a*y |
a = a*y |
b = b*i*(i-1) |
b = b*i*(i - 1) |
c = a/b |
c = a/b |
if(c==0){scale=t; return(s/1)} |
if (c == 0) { |
s = s+c |
scale = t |
|
return (s/1) |
|
} |
|
s = s + c |
} |
} |
} |
} |
|
|
define c(x){ |
define c(x) { |
auto t |
auto t |
t = scale |
t = scale |
scale = scale+1 |
scale = scale + 1 |
x = s(x+2*a(1)) |
x = s(x + 2*a(1)) |
scale = t |
scale = t |
return(x/1) |
return (x/1) |
} |
} |
|
|
define a(x){ |
define a(x) { |
auto a, b, c, d, e, f, g, s, t |
auto a, b, c, d, e, f, g, s, t |
if(x==0) return(0) |
if (x == 0) return(0) |
if(x==1) { |
if (x == 1) { |
if(scale<52) { |
if (scale < 52) { |
return(.7853981633974483096156608458198757210492923498437764/1) |
return(.7853981633974483096156608458198757210492923498437764/1) |
} |
} |
} |
} |
t = scale |
t = scale |
f=1 |
f = 1 |
while(x > .5){ |
while (x > .5) { |
scale = scale + 1 |
scale = scale + 1 |
x= -(1-sqrt(1.+x*x))/x |
x = -(1 - sqrt(1. + x*x))/x |
f=f*2 |
f = f*2 |
} |
} |
while(x < -.5){ |
while (x < -.5) { |
scale = scale + 1 |
scale = scale + 1 |
x = -(1-sqrt(1.+x*x))/x |
x = -(1 - sqrt(1. + x*x))/x |
f=f*2 |
f = f*2 |
} |
} |
s = -x*x |
s = -x*x |
b = f |
b = f |
c = f |
c = f |
d = 1 |
d = 1 |
e = 1 |
e = 1 |
for(a=3;1==1;a=a+2){ |
for (a = 3; 1 == 1; a = a + 2) { |
b=b*s |
b = b*s |
c=c*a+d*b |
c = c*a + d*b |
d=d*a |
d = d*a |
g=c/d |
g = c/d |
if(g==e){ |
if (g == e) { |
scale = t |
scale = t |
return(x*c/d) |
return (x*c/d) |
} |
} |
e=g |
e = g |
} |
} |
} |
} |
|
|
define j(n,x){ |
define j(n,x) { |
auto a,b,c,d,e,g,i,s,k,t |
auto a,b,c,d,e,g,i,s,k,t |
|
|
t = scale |
t = scale |
k = 1.36*x + 1.16*t - n |
k = 1.36*x + 1.16*t - n |
k = length(k) - scale(k) |
k = length(k) - scale(k) |
if(k>0) scale = scale + k |
if (k > 0) scale = scale + k |
|
|
s= -x*x/4 |
s = -x*x/4 |
if(n<0){ |
if (n < 0) { |
n= -n |
n = -n |
x= -x |
x = -x |
} |
} |
a=1 |
a = 1 |
c=1 |
c = 1 |
for(i=1;i<=n;i++){ |
for (i = 1; i <= n; i++) { |
a=a*x |
a = a*x |
c = c*2*i |
c = c*2*i |
} |
} |
b=a |
b = a |
d=1 |
d = 1 |
e=1 |
e = 1 |
for(i=1;1;i++){ |
for (i = 1; 1; i++) { |
a=a*s |
a = a*s |
b=b*i*(n+i) + a |
b = b*i*(n + i) + a |
c=c*i*(n+i) |
c = c*i*(n + i) |
g=b/c |
g = b/c |
if(g==e){ |
if (g == e) { |
scale = t |
scale = t |
return(g/1) |
return (g/1) |
} |
} |
e=g |
e = g |
} |
} |
} |
} |