function radiu2 (x, y, x0, y0)
c
c return the square of the radial distance of (x,y) from (x0,y0)
c
	real x, y, x0, y0
c
	radiu2 = ( (x-x0)**2 + (y-y0)**2 )
c
	return
	end