# ------ funzione per disegnare un piano cartesiano vuoto
piano.cartesiano <- function(Z=10) {
   plot(NULL, xlim=c(-Z,Z), ylim=c(-Z,Z), xlab='x', ylab='y', asp=1)
   grid()
   abline(h=0)
   abline(v=0)
}

piano.cartesiano(10)

# valori di phi per tracciare il cerchio mediante
# le equazioni parametriche

n=101
phi <- seq(0, 2*pi, len=n)  

while(1) {   # loop infinito: attenzione!

  cat(sprintf("\n Clicca su due punti (coincidenti per uscire)\n"))

  p <- locator(2)

  # CALCOLARE r come distanza fra due punti nel piano
  # (per le coordinate dei punti vedi ossatura_cerchi_2p.interattivo.R)

  r <- .....
  
  if ( r < 0.2) break    # se è 'piccolo' vuol dire che vogliamo uscire

  # CALCOLARE i vettori x e y che contengono i punti sulla circonferenza
  x <- .....
  y <- .....

  colore <- sample(rainbow(7))[1]

  # centro
  points(p$x[1], p$y[1], pch=3, col=colore)
  # circonferenza
  points(x, y, ty='l', lwd=2, col= colore)

}