N=100; nf <- rep(1,2); for(k in 3:N) nf[k] <- nf[k-1] + nf[k-2]  
plot(nf, log='y')
grid(lwd=4)
points(2^(1:N), col='blue') 
points(10^(1:N), col='blue') 
points(10^(1:N), col='red') 
points(exp(1:N), col='green') 
plot(nf, log='y', ylim=c(.9,1200) )
grid(lwd=4)
points(2^(1:N), col='blue') 
points(10^(1:N), col='blue') 
points(10^(1:N), col='red') 
points(exp(1:N), col='green') 
N=20; nf <- rep(1,2); for(k in 3:N) nf[k] <- nf[k-1] + nf[k-2]  
plot(nf, log='y', ylim=c(.9,1200) )
grid(lwd=4)
points(2^(1:N), col='blue') 
points(10^(1:N), col='blue') 
points(10^(1:N), col='red') 
points(exp(1:N), col='green') 
plot(log10(nf), )
grig(lwd=1.5)
grid(lwd=1.5)
N=100; nf <- rep(1,2); for(k in 3:N) nf[k] <- nf[k-1] + nf[k-2]  
plot(nf, log='y')
grid(lwd=4)
points(2^(1:N), col='blue') 
points(10^(1:N), col='blue') 
points(10^(1:N), col='red') 
points(exp(1:N), col='green') 
points(10^( 0.5* (1:N) ), col='blue', ty='l') 
points(10^( 0.1* (1:N) ), col='blue', ty='l') 
points(10^( 0.3* (1:N) ), col='blue', ty='l') 
points(10^( log10(2) * (1:N) ), col='blue', ty='l') 
points(10^( log10(exp(1)) * (1:N) ), col='blue', ty='l') 
points(exp( log(2) * (1:N) ), col='green', ty='l') 
points(exp( log(2) * (1:N) ), col='green', ty='l', lwd=2) 
N=100; nf <- rep(1,2); for(k in 3:N) nf[k] <- nf[k-1] + nf[k-2]  
plot(nf, log='y')
grid(lwd=4)
points(2^(1:N), col='blue') 
points(10^(1:N), col='blue') 
points(10^(1:N), col='red') 
points(exp(1:N), col='green') 
points(exp(1:N), col='green', ty='l', lwd=2) 
points(exp( log(10) * (1:N) ), col='green', ty='l', lwd=2) 
points(exp( log(2) * (1:N) ), col='green', ty='l', lwd=2) 
phi = (1+sqrt(5))/2
phi
points(exp( log(phi) * (1:N) ), col='green', ty='l', lwd=2) 
1/sqrt(5)
1/sqrt(5) -> alpha
points(alpha * exp( log(phi) * (1:N) ), col='green', ty='l', lwd=2) 
points(100 * alpha * exp( log(phi) * (1:N) ), col='green', ty='l', lwd=2) 
points(100^-1 * alpha * exp( log(phi) * (1:N) ), col='green', ty='l', lwd=2) 
points(1000^-1 * alpha * exp( log(phi) * (1:N) ), col='green', ty='l', lwd=2) 
points(1000000^-1 * alpha * exp( log(phi) * (1:N) ), col='green', ty='l', lwd=2) 
points(1000000 * alpha * exp( log(phi) * (1:N) ), col='green', ty='l', lwd=2) 
points(1000000 * alpha * exp( - log(phi) * (1:N) ), col='green', ty='l', lwd=2) 
points(exp( (1:N) ) , col='red', ty='l', lwd=2) 
points(exp( 0.001*(1:N) ) , col='red', ty='l', lwd=2) 
points( 0.001* exp((1:N) ) , col='red', ty='l', lwd=2) 
points( exp( -0.5 *  (1:N) ) , col='red', ty='l', lwd=2) 
points( 10^5* exp( -0.5 *  (1:N) ) , col='red', ty='l', lwd=2) 
points( 10^15* exp( -0.5 *  (1:N) ) , col='red', ty='l', lwd=2) 
x <- 1:15
y1 <- c(16,27,45,74,122,201,331,546,900,1484,2447,4034,6651,10966,18080)
y2 <- c(6.77e+00, 9.16e-01, 1.24e-01, 1.68e-02, 2.27e-03, 3.07e-04,
        4.16e-05, 5.63e-06, 7.61e-07, 1.03e-07, 1.39e-08, 1.89e-09,
2.55e-10, 3.46e-11, 4.68e-12)
plot(x, y1)
plot(x, y1, log='y')
plot(x, y2, log='y')
plot(x, y2)
plot(x, y2, log='y')
locator(2)
source("prova.R")
par
 readline(prompt="da a, b, c separati da virgola \n")  -> linea 
linea
strsplit(linea, ',')  
 unlist( strsplit(linea, ',')  )
 as.numeric(  unlist( strsplit(linea, ',')  ) )
q()