Write a Python code using the below formula

to compute 1001 points randomly from a uniform distribution on a unit square. Fix one of these points and compute the square of the distance of all 1000 points in the interval of 0 D² 2. Compute the histogram of D² on [0,2].

Pairwise distances in R2 - Consider Vc(2,1) : Generate 1001 independent, identically distributed in Vc(2,1). - Fix one of the point and call it x=(x1,x2)T. - Compute for each of the rest of 1000 points D2(x,y)=[(x1y1)2+(x2y2)2](y=x). - Clearly 0D2(x,y)2 for all y=x since x1y11 and x2y21. - Histogram of D2(x,y) is fully supported on [0,2]. 1