"pete" <pfil...@mindspring.com> ha scritto nel messaggio
news:4F8A5BF2.17C5@mindspring.com...

#include <stdio.h>
#include <math.h>

#define  u32   unsigned
#define  u64   unsigned __int64
#define  i32   int
#define  u16   unsigned short
#define  u8    unsigned char
#define  i8    char
#define  i16   short
#define  sdC  __stdcall

#define  F   for
#define  R   return
#define  W   while
#define  G   goto
#define  P   printf

double ke=2.7182818284590452353602874713527;
int    con;

long double fs_expl(long double x)
{
    long unsigned n, square;
    long double b, e, old;

    con=0;
    for (square = 0; x > 1; x /= 2) {++con;
        ++square;
    }
    while (-1 > x) {++con;
        ++square;
        x /= 2;
    }
    e = b = n = 1;
    do {++con;
        b /= n++;
        b *= x;
        e += b;
        b /= n++;
        b *= x;
        old = e;
        e += b;
    } while (e > old);
    while (square-- != 0) {++con;
        e *= e;
    }
    return e;

 i=0;
 if(r==0)         R   0;
 *r=0.0;
 if(espon==0.0) {*r=1.0;
                  R   1;}
 else if(espon<0) x=-espon;
 else             x= espon;
 fp=modf(x, &ip);  t=ke;
 if(ip>0xFFFFFFF) R   0;
 F(w=1, v=ip; w<v; ++w)
         {++i; t=t*ke;}
 ip=t;
 v=2; t=1+fp; x=fp;
 F(;;){++i;
       x= x*(fp/v);
       if(x==0.0)  break;
       t+=x; ++v;
      }
 t*=ip;
 if(espon>=0) *r=    t;
 else         *r=1.0/t;
 R   i;