Return to jn.c CVS log

Up to [MW Coherent from dump] / coherent / b / lib / libm

Annotation of coherent/b/lib/libm/jn.c, revision 1.1.1.1

1.1       root        1: /*
                      2:  * Evaluate the Bessel function of the first kind for the
                      3:  * n'th order.
                      4:  */
                      5: #include <math.h>
                      6: 
                      7: double
                      8: jn(n, x)
                      9: register int n;
                     10: double x;
                     11: {
                     12:        double r, ja, jb;
                     13:        register int m, i;
                     14: 
                     15:        if (x == 0.0)
                     16:                return (0.0);
                     17:        m = 0;
                     18:        if (x < 0.0) {
                     19:                x = -x;
                     20:                m = 1;
                     21:        }
                     22:        if (n < 0) {
                     23:                n = -n;
                     24:                if (n&1)
                     25:                        m ^= 1;
                     26:        }
                     27:        switch (n) {
                     28:        case 0:
                     29:                r = j0(x);
                     30:                break;
                     31:        case 1:
                     32:                r = j1(x);
                     33:                break;
                     34:        default:
                     35:                jb = j0(x);
                     36:                ja = j1(x);
                     37:                for (i=1; i<n; i++) {
                     38:                        r = (2.0*i/x)*ja - jb;
                     39:                        jb = ja;
                     40:                        ja = r;
                     41:                }
                     42:                break;
                     43:        }
                     44:        if (m)
                     45:                r = -r;
                     46:        return (r);
                     47: }

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