BSD 3 development

[unix-history] / usr / lib / lib.b

diff --git a/usr/lib/lib.b b/usr/lib/lib.b
index 902162e..e09bf26 100644 (file)
--- a/usr/lib/lib.b
+++ b/usr/lib/lib.b
@@ -1,8 +1,23 @@
 scale = 20
 define e(x){
 scale = 20
 define e(x){

-       auto a, b, c, d, e, g
+       auto a, b, c, d, e, g, w, y
+
+       t = scale
+       scale = t + .434*x + 1
+
+       w = 0
+       if(x<0){
+               x = -x
+               w = 1
+       }
+       y = 0
+       while(x>2){
+               x = x/2
+               y = y + 1
+       }
+

        a=1
        a=1

-       b=10
+       b=1

        c=b
        d=1
        e=1
        c=b
        d=1
        e=1


@@ -11,69 +26,44 @@ define e(x){
                c=c*a+b
                d=d*a
                g = c/d
                c=c*a+b
                d=d*a
                g = c/d

-               if(g == e) return(g/10)
-               e=g
-       }
-}
-scale = 20
-define s(x){
-       auto a, b, c, d, e, g, y
-       y = -x*x
-       a=1
-       b=x
-       c=b
-       d=1
-       e=1
-       for(a=3;1==1;a=a+2){
-               b=b*y
-               c=c*a*(a-1) + b
-               d=d*a*(a-1)
-               g=c/d
-               if(g==e) return(g)
-               e=g
-       }
-}
-scale = 20
-define c(x){
-       auto a, b, c, d, e, g, y
-       y = -x*x
-       a=1
-       b=1
-       c=b
-       d=1
-       e=1
-       for(a=2;1==1;a=a+2){
-               b=b*y
-               c=c*a*(a-1) + b
-               d=d*a*(a-1)
-               g=c/d
-               if(g==e) return(g)
+               if(g == e){
+                       g = g/1
+                       while(y--){
+                               g = g*g
+                       }
+                       scale = t
+                       if(w==1) return(1/g)
+                       return(g/1)
+               }

                e=g
        }
 }
                e=g
        }
 }

-scale = 20
+

 define l(x){
        auto a, b, c, d, e, f, g, u, s, t
        if(x <=0) return(1-10^scale)
        t = scale
 define l(x){
        auto a, b, c, d, e, f, g, u, s, t
        if(x <=0) return(1-10^scale)
        t = scale

-       scale = 0
-       f = 1
-       s = x
-       while(s > 0){
-               s = s/10
-               f = f + 1
-       }
-       scale = t + f
+

        f=1
        f=1

+       scale = scale + scale(x) - length(x) + 1
+       s=scale

        while(x > 2){
        while(x > 2){

+               s = s + (length(x)-scale(x))/2 + 1
+               if(s>0) scale = s

                x = sqrt(x)
                f=f*2
        }
        while(x < .5){
                x = sqrt(x)
                f=f*2
        }
        while(x < .5){

+               s = s + (length(x)-scale(x))/2 + 1
+               if(s>0) scale = s

                x = sqrt(x)
                f=f*2
        }
                x = sqrt(x)
                f=f*2
        }

+
+       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)

        s = u*u
        b = 2*f
        c = b
        s = u*u
        b = 2*f
        c = b


@@ -91,10 +81,51 @@ define l(x){
                e=g
        }
 }
                e=g
        }
 }

-scale = 20
+
+define s(x){
+       auto a, b, c, s, t, y, p, n, i
+       t = scale
+       y = x/.7853
+       s = t + length(y) - scale(y)
+       if(s<t) s=t
+       scale = s
+       p = a(1)
+
+       scale = 0
+       if(x>=0) n = (x/(2*p)+1)/2
+       if(x<0) n = (x/(2*p)-1)/2
+       x = x - 4*n*p
+       if(n%2!=0) x = -x
+
+       scale = t + length(1.2*t) - scale(1.2*t)
+       y = -x*x
+       a = x
+       b = 1
+       s = x
+       for(i=3; 1==1; i=i+2){
+               a = a*y
+               b = b*i*(i-1)
+               c = a/b
+               if(c==0){scale=t; return(s/1)}
+               s = s+c
+       }
+}
+
+define c(x){
+       auto t
+       t = scale
+       scale = scale+1
+       x = s(x+2*a(1))
+       scale = t
+       return(x/1)
+}
+

 define a(x){
        auto a, b, c, d, e, f, g, s, t
        if(x==0) return(0)
 define a(x){
        auto a, b, c, d, e, f, g, s, t
        if(x==0) return(0)

+       if(x==1)
+               if(scale<52)
+return(.7853981633974483096156608458198757210492923498437764/1)

        t = scale
        f=1
        while(x > .5){
        t = scale
        f=1
        while(x > .5){


@@ -124,9 +155,15 @@ define a(x){
                e=g
        }
 }
                e=g
        }
 }

-scale = 20
+

 define j(n,x){
 define j(n,x){

-auto a,b,c,d,e,g,i,s
+auto a,b,c,d,e,g,i,s,k,t
+
+       t = scale
+       k = 1.36*x + 1.16*t - n
+       k = length(k) - scale(k)
+       if(k>0) scale = scale + k
+

 s= -x*x/4
 if(n<0){
        n= -n
 s= -x*x/4
 if(n<0){
        n= -n


@@ -147,7 +184,8 @@ for(i=1;1;i++){
        c=c*i*(n+i)
        g=b/c
        if(g==e){
        c=c*i*(n+i)
        g=b/c
        if(g==e){

-               return(g)
+               scale = t
+               return(g/1)

                }
        e=g
        }
                }
        e=g
        }