/* u.c - unit step function */

int u(x)
double x;
{
       if (x >= 0)
              return 1;
       else
              return 0;
}