****** check_fit in accord with the 'StarLight Mannual P40 Figure 3.' ********     
	real z(10),fra(10),err(10),fra_group(10),err_group(10)
	real fra_group_p(10),err_group_p(10)
	integer all_agn(10),group_agn(10),field_agn(10),group_gal(10),field_gal(10)
	integer i

      open(1,file='../data1/field_bin1.tbl',status='old')
   
      do i=1,10
         read(1,101)z(i),all_agn(i),group_agn(i),field_agn(i),
     +group_gal(i),field_gal(i),fra(i),err(i),fra_group(i),err_group(i)
     +,fra_group_p(i),err_group_p(i)

	fra(i)=fra(i)*100.
	err(i)=err(i)*100.

	fra_group(i)=fra_group(i)*100.
	err_group(i)=err_group(i)*100.

	 if(i.eq.1)print*,z(i),fra(i),err(i),fra_group(i),err_group(i)
         
         enddo     

101	format(f4.2,1x,i2,1x,i1,1x,i2,1x,i4,1x,i5,1x,f8.6,1x,f8.6
     +         ,1x,f8.6,1x,f8.6,2e11.4)

	call fit(z,fra,10,err,0,a,b,siga,sigb,chi2,q)
	print*, 'a=',a,'b=',b,'siga=',siga,'sigb=',sigb

	call pearsn(z,fra,10,r,prob,z00)
	print*,'r=',r,'prob=',prob

	call fit(z,fra_group,10,err_group,0,a1,b1,siga1,sigb1,chi22,q1)
      print*, 'a1=',a1,'b1=',b1,'siga1=',siga1,'sigb1=',sigb1

	call pearsn(z,fra_group,10,r1,prob1,z002)
	print*,'r1=',r1,'prob1=',prob1

     
************** begin to plot *******************

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

************ plot the AGNfraction_z **********

      call pgenv(0.,1.,-0.5,2.,0,0)
      call pglabel('\fiz\d\frmean\u','\fif\d\frAGN\u(%)','')

c        call pgtext(0.25,0.1,'r\ds\u=0.938;')
c	call pgtext(0.52,0.1,'slope=0.005\(2233)0.001')
	

      do i=1,10  
	call pgsci(2)        
	call pgpt1(z(i),fra(i),18)
	call pgerry(1,z(i),fra(i)+0.5*err(i),fra(i)-0.5*err(i),1.2)
	end do
	
	do i=1,10

	call pgsci(1)
	call pgpt1(z(i),fra_group(i),23)
	call pgerry(1,z(i),fra_group(i)+0.5*err_group(i),
     +              fra_group(i)-0.5*err_group(i),1.2)

	end do

     	call pgsci(2)
	call pgmove(0.0,a+b*0.0)
	call pgdraw(1.2,a+b*1.2)

	call pgsci(1)
	call pgmove(0.0,a1+b1*0.0)
	call pgdraw(1.2,a1+b1*1.2)

        call pgsch(1.3)
        call pgsci(2)
        call pgpt1(0.1,1.7,18)
        call pgtext(0.1,1.67,'  \fif\d\frAGN\u in field')
        call pgtext(0.2,-0.2,'\fir\d\frs\u=0.945')
        call pgtext(0.2,-0.35,'slope=0.414')
        call pgsci(1)
        call pgpt1(0.1,1.45,23)
        call pgtext(0.1,1.42,'  \fif\d\frAGN\u in \fiP\d\frmem\u>0 group')
        call pgtext(0.6,-0.2,'\fir\d\frs\u=0.461')
        call pgtext(0.6,-0.35,'slope=0.342')

*************************************
        call pgend
        

        close(1)
        close(2)

        end