k <- 1:10
k
1/2^k
sum(1/2^k)
cumsum(1/2^k)
nf <- rep(1,)
nf
nf <- rep(1,2)
nf
N=100; nf <- rep(1,2); for(k in 3:N) nf[k] <- nf[k-1] + nf[k-2]  
N=10; nf <- rep(1,2); for(k in 3:N) nf[k] <- nf[k-1] + nf[k-2]  
nf
N=100; nf <- rep(1,2); for(k in 3:N) nf[k] <- nf[k-1] + nf[k-2]  
N=10; nf <- rep(1,2); for(k in 3:N) nf[k] <- nf[k-1] + nf[k-2]  
plot(nf)
N=100; nf <- rep(1,2); for(k in 3:N) nf[k] <- nf[k-1] + nf[k-2]  
plot(nf)
plot(nf, log='y')
N=20; nf <- rep(1,2); for(k in 3:N) nf[k] <- nf[k-1] + nf[k-2]  
plot(nf, log='y')
grid()
grid(lwd=2)
grid(lwd=4)
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') 
nf[51]/nf[50]
source("int_rette.R") 
m1
m2
q1
q2
source("int_rette.R") 
source("int_rette.R") 
source("int_rette.R") 
source("int_rette.R") 
( p <- locator(2) )
p$x
p$y
q()