Return to math.c CVS log

Up to [CSRG BSD Unix] / 43BSD / contrib / icon / pifuncs

Annotation of 43BSD/contrib/icon/pifuncs/math.c, revision 1.1.1.1

1.1       root        1: /*
                      2: #      MATH(3.icon)
                      3: #
                      4: #      Miscellaneous math functions
                      5: #
                      6: #      Ralph E. Griswold
                      7: #
                      8: #      Last modified 8/19/84
                      9: #
                     10: */
                     11: 
                     12: #include "../h/rt.h"
                     13: #include <errno.h>
                     14: 
                     15: int errno;
                     16: /*
                     17:  * exp(x), x in radians
                     18:  */
                     19: Xexp(nargs, arg1, arg0)
                     20: int nargs;
                     21: struct descrip arg1, arg0;
                     22:    {
                     23:    int t;
                     24:    double y;
                     25:    union numeric r;
                     26:    double exp();
                     27:    
                     28:    if ((t = cvreal(&arg1, &r)) == NULL) runerr(102, &arg1);
                     29:    y = exp(r.real);
                     30:    if (errno == ERANGE) runerr(252, NULL);
                     31:    mkreal(y,&arg0);
                     32:    }
                     33: Procblock(exp,1)
                     34: 
                     35: /*
                     36:  * log(x), x in radians
                     37:  */
                     38: Xlog(nargs, arg1, arg0)
                     39: int nargs;
                     40: struct descrip arg1, arg0;
                     41:    {
                     42:    int t;
                     43:    double y;
                     44:    union numeric r;
                     45:    double log();
                     46:    
                     47:    if ((t = cvreal(&arg1, &r)) == NULL) runerr(102, &arg1);
                     48:    y = log(r.real);
                     49:    if (errno == EDOM) runerr(251, NULL);
                     50:    mkreal(y,&arg0);
                     51:    }
                     52: Procblock(log,1)
                     53: 
                     54: /*
                     55:  * log10(x), x in radians
                     56:  */
                     57: Xlog10(nargs, arg1, arg0)
                     58: int nargs;
                     59: struct descrip arg1, arg0;
                     60:    {
                     61:    int t;
                     62:    double y;
                     63:    union numeric r;
                     64:    double log10();
                     65:    
                     66:    if ((t = cvreal(&arg1, &r)) == NULL) runerr(102, &arg1);
                     67:    y = log10(r.real);
                     68:    if (errno == EDOM) runerr(251, NULL);
                     69:    mkreal(y,&arg0);
                     70:    }
                     71: Procblock(log10,1)
                     72: 
                     73: /*
                     74:  * sqrt(x), x in radians
                     75:  */
                     76: Xsqrt(nargs, arg1, arg0)
                     77: int nargs;
                     78: struct descrip arg1, arg0;
                     79:    {
                     80:    int t;
                     81:    double y;
                     82:    union numeric r;
                     83:    double sqrt();
                     84:    
                     85:    if ((t = cvreal(&arg1, &r)) == NULL) runerr(102, &arg1);
                     86:    y = sqrt(r.real);
                     87:    if (errno == EDOM) runerr(251, NULL);
                     88:    mkreal(y,&arg0);
                     89:    }
                     90: Procblock(sqrt,1)

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