1
2subroutine grexp(n, x, a)
3integer i,n
4double precision x(n),a
5call fseedi()
6do i = 1,n
7   call frexp(x(i), a)
8call fseedo()
9return
10end
11
12
13