real a(1000,1000)
	
	open(10,file="out.dat")
	r=800
	x1=-200
	y1=500
	x2=1200
	y2=500
	
	do i=1,1000
	do j=1,1000
	if( sqrt((i-x1)**2+(j-y1)**2).lt.R.
     *  and.sqrt((i-x2)**2+(j-y2)**2).lt.R) then
        a(i,j)=2.
        else
        a(i,j)=1.
	endif
	enddo
	enddo
	write(10,"(1000F15.5)") a
	stop
	end