# carichiamo il file seguente, per battere meno comandi
source("http://macosa.dima.unige.it/R/curve.txt")
# quindi ...
boxm(-5,5, -5,5)
# un po' di punti del cerchio
t <- seq(0,2*pi, len=5000)
x1 <- cos(t); y1 <- sin(t)
punt(x1,y1,"red")
# i loro trasformati, usando il primo valore di q
q <- 1/2
x2 <- x1/(y1-x1-q)
y2 <- (x1+y1)/(y1-x1-q)
punt(x2,y2, "blue")
# usando il secondo valore di q
q <- sqrt(2)
x3 <- x1/(y1-x1-q)
y3 <- (x1+y1)/(y1-x1-q)
punt(x3,y3, "green4")
#
# Se invece di sqrt(2) prendo ...
boxm(-5,20, -10,15)
punt(x1,y1, "red")
punt(x3,y3, "green4")
q <- 1.45
x4 <- x1/(y1-x1-q)
y4 <- (x1+y1)/(y1-x1-q)
punt(x4,y4,"brown")