integer ic_gal(4,10),ic_agn(10)
	real ratio(4,10),d_ratio(4,10)
	real complete(4),purity(4)
c	data ic/4*0/
	data complete/0.97,0.94,0.92,0.80/
	data purity/0.63,0.67,0.69,0.78/

	open(1,file='../data/gal_count_zbin.tbl',status='old')
	read(1,*)
	do 100 i=1,4
	read(1,*) (ic_gal(i,j),j=1,10)
100	continue
	read(1,*) ic_agn
	close(1)

	  do 200 ip=1,4
	     do 400 iz=1,10
             ratio(ip,iz)=1.0*ic_agn(iz)/ic_gal(ip,iz)
	      if(ic_agn(iz).eq.0) then
		   d_ratio(ip,iz)=0.0 	
		else
	     d_ratio(ip,iz)=ratio(ip,iz)
     +          *sqrt(1.0/ic_agn(iz)+1.0/ic_gal(ip,iz))
     +          *sqrt(complete(ip)/purity(ip))
		end if
400	   continue
200	  continue	  
          
	do ip=1,4
	    write(*,'(i3,10e11.4/3x,10e11.4)') ip, (ratio(ip,iz),iz=1,10),
     +	          (d_ratio(ip,iz),iz=1,10)
	end do

	end