MVLib - Short Writeup - Version 4.0
Example Number 33 - Normalized histogram
The subroutines mvbookget and mvbookget allow a simple way to normalize an histogram providing the so called PDF (Probability distribution Function).
The fortran code is available here and the complete script running the code is available here.
Previous Example List of Examples Next Example
MVLib Home Page 
	real v(40) ; character title*60 
	call  mvbook(10,0.0,10.0,40,' ')
	do i=1,100000
	   call mvbookfill(10,10*gauss(0.5,0.2),1.0)
	enddo
	call mvbookplot(10)
	perc=22 ; call mvbookcentile(10,perc,xx)
	write(title,'(a,i2,a,f3.1)') 'The ',nint(perc),'th percentile is ',xx
	call displayexample('example33','Histogram of observed distribution',title)
	call mvbookget (10,v)
	tot=0 ; do i=1,20 ; tot=tot+v(i) ; enddo ; v=v/tot
	call mvbookset (10,v)
	call mvbookplot(10)
	call displayexample('example33','Probability distribution function',' ')
	end