Return to lib.b CVS log

Up to [MW Coherent from dump] / coherent / a / usr / lib

Annotation of coherent/a/usr/lib/lib.b, revision 1.1

1.1     ! root        1: /*
        !             2:  * Definition of pi.
        !             3:  * Used for large signs and funny
        !             4:  * arctangents.
        !             5:  */
        !             6: pi = 3.14159265358979323846264338327950288419716939937510582097494459230781640628620899862803482534211706798
        !             7: 
        !             8: /*
        !             9:  * exp(x) -- exponential function
        !            10:  * to the scale that is in effect
        !            11:  * at the time of the call.
        !            12:  */
        !            13: define exp(x) {
        !            14:        auto i, num, den, term, result;
        !            15: 
        !            16:        num = den = result = 1;
        !            17:        for (i=1; ; i++) {
        !            18:                num *= x;
        !            19:                den *= i;
        !            20:                if ((term = num/den) == 0)      return (result);
        !            21:                result += term;
        !            22:        }
        !            23: }
        !            24: 
        !            25: /*
        !            26:  * ln(x) -- natural logarithm
        !            27:  */
        !            28: define ln(x) {
        !            29:        auto i, xp1, xm1, num, den, term, result;
        !            30: 
        !            31:        if (x <= 0)     return (0);
        !            32:        xp1 = x+1;
        !            33:        xm1 = x-1;
        !            34:        num = 2*xm1;
        !            35:        den = xp1;
        !            36:        result = num/den;
        !            37:        for (i=3; ; i+=2) {
        !            38:                num *= xm1*xm1;
        !            39:                den *= xp1*xp1;
        !            40:                if ((term = num/(i*den)) == 0)  return (result);
        !            41:                result += term;
        !            42:        }
        !            43: }
        !            44: 
        !            45: /*
        !            46:  * sin(x) -- sine function.
        !            47:  */
        !            48: define sin(x) {
        !            49:        auto i, num, den, term, result, sign, pi2, minusx2;
        !            50: 
        !            51:        sign = den = 1;
        !            52:        pi2 = 2*pi;
        !            53:        if (x < 0) {
        !            54:                sign = -1;
        !            55:                x = -x;
        !            56:        }
        !            57:        while (x > pi2)         x -= pi2;
        !            58:        num = result = x;
        !            59:        minusx2 = -x*x;
        !            60:        for (i=3; ; i+=2) {
        !            61:                num *= minusx2;
        !            62:                den *= (i-1)*i;
        !            63:                if ((term = num/den) == 0)      return (sign*result);
        !            64:                result += term;
        !            65:        }
        !            66: }
        !            67: 
        !            68: /*
        !            69:  * cos(x) -- cosine function.
        !            70:  */
        !            71: define cos(x) {
        !            72:        auto i, num, den, term, result, pi2, minusx2;
        !            73: 
        !            74:        pi2 = 2*pi;
        !            75:        if (x < 0)      x = -x;
        !            76:        while (x > pi2)         x -= pi2;
        !            77:        num = den = result = 1;
        !            78:        minusx2 = -x*x;
        !            79:        for (i=2; ; i+=2) {
        !            80:                num *= minusx2;
        !            81:                den *= i*(i-1);
        !            82:                if ((term = num/den) == 0)      return (result);
        !            83:                result += term;
        !            84:        }
        !            85: }
        !            86: 
        !            87: /*
        !            88:  * atan(x) -- arctangent function of x
        !            89:  */
        !            90: define atan(x) {
        !            91:        auto i, sign, inverse, num, den, term, x2, x2p1, result;
        !            92: 
        !            93:        sign = 1;
        !            94:        inverse = 0;
        !            95:        if (x < 0) {
        !            96:                sign = -1;
        !            97:                x = -x;
        !            98:        }
        !            99:        if (x > 1) {
        !           100:                inverse = 1;
        !           101:                x = 1/x;
        !           102:        }
        !           103:        x2 = x*x;
        !           104:        x2p1 = x2+1;
        !           105:        result = num = den = 1;
        !           106:        for (i=2; ; i+=2) {
        !           107:                num *= x2*i;
        !           108:                den *= x2p1*(i+1);
        !           109:                if ((term = num/den) == 0) {
        !           110:                        result = x*result/x2p1;
        !           111:                        if (inverse != 0)       result = pi/2 - result;
        !           112:                        return (sign * result);
        !           113:                }
        !           114:                result += term;
        !           115:        }
        !           116: }
        !           117: 
        !           118: /*
        !           119:  * Bessel functions of first kind of
        !           120:  * integer order `n' for `x'.
        !           121:  */
        !           122: define j(n, x)
        !           123: {
        !           124:        auto i, num, den, result;
        !           125: 
        !           126:        if (n%1 != 0)   return (0);
        !           127:        den = 2;
        !           128:        num = 1;
        !           129:        if (n > 0) {
        !           130:                for (i=0; i<n; i++)     num *= x;
        !           131:        } else if (n < 0) {
        !           132:                for (i=0; i<-n; i++)    den *= x;
        !           133:        }
        !           134:        return (0);
        !           135: }

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.