dane=RandomReal[{-1,1},{10000,2}];
Show[ListPlot[dane,AspectRatio->1],ParametricPlot[{Cos@θ,Sin@θ},{θ,0,2π}]];
Show[ListPlot[{Select[#,Norm[#]>1&],Select[#,Norm[#]<1&]}&@dane,AspectRatio->1,PlotStyle->{Red,Blue}],ParametricPlot[{Cos@θ,Sin@θ},{θ,0,2π}]]
4 Length@Select[#,Norm[#]<1&]/Length@#&@dane//N

(((6^4)-(5^4))/(6^4))//N

(((36^24)-(35^24))/(36^24))//N

gora[x_]:=Chop@(1/14) (({{\[Piecewise], ({{6 (Sqrt[49-(x^2)]), x≥3||x<-3}, {21+12 (Sqrt[10])+7 x-6 (Sqrt[10]) (Sqrt[-(-1+x) (3+x)]), -3≤x<-1}, {21+12 (Sqrt[10])-7 x-6 (Sqrt[10]) (Sqrt[-(-3+x) (1+x)]), 1≤x<3}, {(7/4) (-11 Abs[3-4 x]+6 (12+Abs[1-2 x]+Abs[1+2 x])-11 Abs[3+4 x]), True}})}}));
dol[x_]:=Chop[(1/2)*(3*Sqrt[1-(x/7)^2]+Sqrt[1-(Abs[Abs[x]-2]-1)^2]+Abs[x/2]-((3*Sqrt[33]-7)/112)*x^2-3)*((x+4)/Abs[x+4]-(x-4)/Abs[x-4])-3*Sqrt[1-(x/7)^2]]

(955/48-2/7 (-3 Sqrt[10]+2 Sqrt[33]+7 π+3 Sqrt[10] π)+21 ArcCos[3/7]+21 ArcCos[4/7])
%//N

Plot[{gora[x],dol[x]},{x,-7,7},PlotStyle->Gray]

n=1000;
dane=Partition[Riffle[RandomReal[{-7,7},{n}],RandomReal[{-3,3},{n}]],2];
bat=Select[dane,If[#[[2]]>0,#[[2]]<gora[#[[1]]],#[[2]]>dol[#[[1]]]]&];
Length@bat/n*(14*6)//N
Show[ListPlot[{dane,bat},PlotStyle->{Black,Yellow}],Plot[{gora[x],dol[x]},{x,-7,7},PlotStyle->Gray]]

Show[Histogram[RandomReal[{0,1},{100}],{0.1},"CDF"],Plot[x,{x,0,1}]]