Con le istruzioni allegate mediante R ottieni l'immagine qui a destra.   Aggiungi delle istruzioni in modo da ottenere l'immagine qui a sinistra (eventualmente senza colorarne l'interno).

source("http://macosa.dima.unige.it/r.R")
BF=3; HF=3; PLANE(-3,3, -3,3)
PLANE(-3,3, -3,3)
circl(0,0,1, "black"); g=degrees
P=c(cos(30*g),sin(30*g)); Q=c(cos((30+120)*g),sin((30+120)*g)); R=c(cos((30+240)*g),sin((30+240)*g))
text(P[1]+0.4,P[2]+0.4,"P"); text(Q[1]-0.4,Q[2]+0.4,"Q"); text(R[1],R[2]-0.5,"R")
d=point_point(P[1],P[2], Q[1],Q[2]); d
circl(R[1],R[2],d, "seagreen"); circl(Q[1],Q[2],d, "blue"); circl(P[1],P[2],d, "red")
POINT(P[1],P[2],"red"); POINT(Q[1],Q[2],"blue"); POINT(R[1],R[2],"seagreen")
A=circle_circle(R[1],R[2],d, P[1],P[2],d);A; B=circle_circle(R[1],R[2],d, Q[1],Q[2],d);B
C=circle_circle(P[1],P[2],d, Q[1],Q[2],d);C
text(A[3]+0.4,A[4]-0.4,"A"); text(B[1]-0.4,B[2]-0.4,"B"); text(C[3],C[4]+0.5,"C")
circl3p(c(C[3],Q[1],P[1]), c(C[4],Q[2],P[2]),"brown" ); C1=center3p
circl3p(c(A[3],P[1],R[1]), c(A[4],P[2],R[2]),"brown" ); A1=center3p
circl3p(c(B[1],R[1],Q[1]), c(B[2],R[2],Q[2]),"brown" ); B1=center3p
POINT(P[1],P[2],"red"); POINT(Q[1],Q[2],"blue"); POINT(R[1],R[2],"seagreen")
POINT(B[1],B[2],"red"); POINT(A[3],A[4],"blue"); POINT(C[3],C[4],"seagreen")
# Aggiungi l'immagine che vuoi ottenere
circle(0,0,1, "brown")
ARC(C1[1],C1[2], 1, 90-120,90+120, "brown"); ARC(A1[1],A1[2],1, -30-120,-30+120,"brown")
ARC(B1[1],B1[2],1, 210-120,210+120,"brown")
ARC(R[1],R[2],d, 0,180, "brown"); ARC(P[1],P[2], d, 120,120+180, "brown")
ARC(Q[1],Q[2], d, 240,240+180, "brown")
# Ottieni solo l'immagine
PLANEww(-2,2,-1.5,2.5)
circl(0,0,1, "brown")
arc(C1[1],C1[2], 1, 90-120,90+120, "brown"); arc(A1[1],A1[2],1, -30-120,-30+120,"brown")
arc(B1[1],B1[2],1, 210-120,210+120,"brown")
arc(R[1],R[2],d, 0,180, "brown"); arc(P[1],P[2], d, 120,120+180, "brown")
arc(Q[1],Q[2], d, 240,240+180, "brown")
#
# Per colorare la figura:
Q1 = function(x,y) (x^2+y^2 < 1 & x^2+(y+1)^2 > d^2)
P1 = function(x,y) (x^2+(y-1)^2 < 1 & ((x-P[1])^2+(y-P[2])^2 > d^2 | (x-Q[1])^2+(y-Q[2])^2 > d^2))
Q2 = function(x,y) (x^2+y^2 < 1 & (x-P[1])^2+(y-P[2])^2 > d^2)
P2 = function(x,y) ((x-B1[1])^2+(y-B1[2])^2 < 1 & ((x-Q[1])^2+(y-Q[2])^2 > d^2 | (x)^2+(y+1)^2 > d^2))
Q3 = function(x,y) (x^2+y^2 < 1 & (x-Q[1])^2+(y-Q[2])^2 > d^2)
P3 = function(x,y) ((x-A1[1])^2+(y-A1[2])^2 < 1 & ((x-P[1])^2+(y-P[2])^2 > d^2 | (x)^2+(y+1)^2 > d^2))
RR = function(x,y) Q1(x,y) | P1(x,y) | Q2(x,y) | P2(x,y) | Q3(x,y) | P3(x,y)
for(i in 1:100) FIGURE(RR, -2,2, -2,2, "red")
circl(0,0,1, "brown")
arc(C1[1],C1[2], 1, 90-120,90+120, "brown"); arc(A1[1],A1[2],1, -30-120,-30+120,"brown")
arc(B1[1],B1[2],1, 210-120,210+120,"brown")
arc(R[1],R[2],d, 0,180, "brown"); arc(P[1],P[2], d, 120,120+180, "brown")
arc(Q[1],Q[2], d, 240,240+180, "brown")