write(*,*) ' Input seed = '
	read(*,*) nseed
	x=rand(nseed)
	sum_x = 0.
	sum_xsq= 0.
	n = 1000

	do i=1, n
	  x=rand()
	  sum_x = sum_x + x
 	  sum_xsq=sum_xsq+x**2.
	end do

	emu=sum_x/n
	y =  (sum_xsq-n*emu**2)/(n-1)
	sigma = sqrt ( y)

	write(*,100) n, emu, sigma
100     format('Uniformly Distributed Random Number 0<-->1'/
     *         'Number of samples = ', i6/
     *         'Mean              = ', f10.5/
     *         'Std Deviation     = ', f10.5/)

	stop
	end