# Quadrica della nota 26.1.
f = function (x,y)
9*x^2-12*x*y+9*y^2

# Radice di f(z).
rfz = function (z)
{x=Re(z); y=Im(z)
sqrt(f(x,y))}

Gr.postscript('27-lem-2.ps',5,5)
larg=7
intervallo=c(-larg/2,larg/2)
Gr(intervallo,intervallo)
alfa=seq(0,2*pi,length=60)
z=exp(1i*alfa)
lines(z)
lines(z/rfz(z))
lines(z*rfz(z))
raggi=Mm(c(z,z*rfz(z),
  rep(NA,length(z))),righe=3)
lines(raggi)
dev.off()