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