subroutine gen_random_seed()
integer :: s, di, clock
integer, dimension(:), allocatable :: seed
call random_seed(size = di)
allocate(seed(di))
call system_clock(count=clock)
seed = clock + 37 * (/ (s - 1, s = 1, di) /)
call random_seed(put=seed)
deallocate(seed)
end subroutine gen_random_seed
!Questo è un procedimento per generare numeri casuali tra 0 e 1 che non si ripetano
!viene utilizzato il valore del clock del computer per renderli realmente casuali.
!Le subroutine usate qui sono già presenti di default in fortran.


program eszero
    real,dimension(:),allocatable :: x,cont,ec        
    real,dimension(2) :: e
    integer :: ni,c,n,d
    real::a
    
    print*, "Inserire il numero di valori desiderati:"
    read*, n
    print*, "Inserire gli estremi dell'istogramma:"
    read*, e
    print*, "Inserire il numero di intervalli:"
    read*, ni
!Definizione di n,e,ni
    
   allocate(x(n),cont(ni),ec(ni+1))
!Prendo spazio in memoria per gli array
    cont=0.0
    a=(e(2)-e(1))/ni
!ampiezza dei bin
    ec(1)=e(1)

	call gen_random_seed()

    call random_number(x)

    do c=1,n
        x(c)=e(1)+(e(2)-e(1))*x(c)
!adatto all'intervallo di valori desiderato
        do d=1,ni
            ec(d+1)=e(1)+a*d
            if(ec(d)<x(c).and.(x(c)<ec(d+1))) cont(d)=1.0+cont(d)
        end do
!Conto quanti valori ci sono in ogni bin
    end do
    
    open(unit=001,file='distunif.x')
    do c=1,ni
        write(unit=001,fmt=*) (ec(c+1)+ec(c))/2.0, cont(c), sqrt(cont(c))
    end do
    close(unit=001)
!Salvo in un file il valore medio, il numero di conteggi e la radice quadrata di esso (usata come incertezza)
    
    open(unit=002,file='istounif.x')
        write(unit=002,fmt=*) "set title 'Istogramma della distribuzione uniforme'"
        write(unit=002,fmt=*) "set xlabel 'Numeri casuali'"
        write(unit=002,fmt=*) "set ylabel 'Numero di conteggi'"
        write(unit=002,fmt=*) "set xrange [",e(1),":",e(2),"]"
        write(unit=002,fmt=*) "set yrange [0:]"
        write(unit=002,fmt=*) "plot 'distunif.x' using 1:2:3 with error, 'distunif.x' using 1:2 with boxes"  
        write(unit=002,fmt=*) "replot 50*(1.0+0.1*cos(3.14*x-6.28))" 
    close(unit=002)
!Salvo in un secondo file le istruzioni da dare a gnuplot per graficare quanto desiderato.
!'1:2:3 with error' significa che uso la prima colonna del file come x, la seconda come y e la terza come incertezza su y.
!'with boxes' serve ad indicare che è un istogramma.
!Sovrappongo poi la curva desiderata.
end program eszero

!Nel terminale devo poi inserire i seguenti comandi:
!gfortran propedprova.f90 -o propro.x (per compilare il programma 'propedprova' nominando l'eseguibile 'propro')
!./propro.x (per eseguirlo)
!(inserisco '1000','2 4','20' quando chiede i valori)
!gnuplot
!load 'istounif.x' (per aprire l'istogramma)

!Ho anche stampato i grafici con i seguenti comandi in seguito a quelli scritti sopra:
!set term png (imposta il terminale nel formato .png)
!set output "istounif_24.png" (nome dell'immagine in output)
!replot (salva l'immagine nel file creato)
!set term x11 (riporta il terminal alle condizioni iniziali)