float gasdev(long *idum);