# Diagramma polare dell'intensità di una lampada a incandescenza
# in funzione dell'angolo (intensità in "candele")
plot(c(-160,160),c(-160,160),type="n",xlab="",ylab="",axes=FALSE,asp=1)
rad <- function(x) x*pi/180; A <- seq(0,330,30)*pi/180
for (i in 1:6) symbols(0,0,circles=i*25, inches=FALSE, add=TRUE, fg="blue")
for(i in 0:6) abline(0,tan(rad(30*i)),lty=3,col="grey50")
int <- c(135,150,125,115,125,105,0,105,115,130,115,130)
polygon(cos(A)*int,sin(A)*int,angle=75,density=10,border="black",col="red")
text(cos(rad(170))*150,sin(rad(170))*150,"150")
text(cos(rad(192))*100,sin(rad(192))*100,"100")