parameter(n=28)
        parameter(m=7984)
        real x(n),y(m)

	open(1,file='../data/agn_pmem.tbl',status='old')
	open(2,file='../data/gal_pmem.tbl',status='old')

	do 10 i=1,n
	 read(1,'(f7.5)')x(i)
10	continue

	do 20 j=1,m
	 read(2,'(e15.4)')y(j)
20	continue


c	call pgbegin(0,'N_Pmem1.eps/vcps',2,1)
	call pgbegin(0,'?',2,1)
c	call pgbegin(0,'N_Pmem1.eps/vcps',1,1)
c	call pgpaper(0.0,1.6)
c	call pgpaper(0.0,0.8)
	call pgpaper(0.0,0.4)
	call pgscf(2)
	call pgsls(1)
	call pgslw(2)
	call pgsch(1.5)

************ plot the N_AGN-P_mem **********
	call pgenv(0.,1.2,0.,25.,0,0)
	call pglabel('\fiP\fr\dmem\u','\fiN\fr\dAGN\u','')
	call pgtext(1.08,1.5,'(a)')

	call pgsci(4)
	call pgslw(4)
	call pghist(n,x,0.,1.3,13,5)
	call pgslw(2)
	
	call pgsci(1)
	call pgtext(0.1,21.,'same histogram for all samples')


************ plot the N_gal-P_mem **********
	call pgsch(1.8)
	call pgenv(0.,1.2,0.,3000.,0,0)
	call pglabel('\fiP\fr\dmem\u','\fiN\fr\dgal\u','')
       
	call pgtext(1.09,200.,'(b)')
        
	call pgsch(1.5)
	call pgsls(2)
	call pgmove(0.1,2600.)
	call pgdraw(0.1,0.0)
	call pgmove(0.3,2200.0)
	call pgdraw(0.3,0.0)
	call pgmove(0.5,1800.0)
	call pgdraw(0.5,0.0)
	call pgmove(0.7,1400.0)
	call pgdraw(0.7,0.0)
	call pgmove(1.0,2600.0)
	call pgdraw(1.0,0.0)

        call pgsls(1)
	

	call pgmove(0.5,1800.)
	call pgdraw(1.0,1800.)
	call pgmove(0.5,1850.)
	call pgdraw(0.5,1750.)
	call pgmove(1.0,1850.)
	call pgdraw(1.0,1750.)
	call pgtext(0.62,1900.,'sample3')

	call pgmove(0.3,2200.)
	call pgdraw(1.0,2200.)
	call pgmove(0.3,2250.)
	call pgdraw(0.3,2150.)
	call pgmove(1.0,2250.)
	call pgdraw(1.0,2150.)
	call pgtext(0.52,2300.,'sample2')

	call pgmove(0.1,2600.)
	call pgdraw(1.0,2600.)
	call pgmove(0.1,2650.)
	call pgdraw(0.1,2550.)
	call pgmove(1.0,2650.)
	call pgdraw(1.0,2550.)
	call pgtext(0.42,2700.,'sample1')

	call pgsci(4)
	call pgslw(4)
        call pghist(m,y,0.,1.3,13,5)
	call pgslw(2)

        call pgend

        close(1)
        end