r <- 0.05
term <- 40
period <- 1/12
payments <- 100
n <- floor(term/period)
pension <- 0
for (i in 1:n) {
pension[i+1] <- pension[i]*(1 + r*period)+ payments
}
time <- (0:n)*period
plot(time, pension,pch=19,cex=0.8,col=1)
plot(time, pension,type="l",col=2,lwd=2)

pop <- 1:1000000
codpop <- pop%%5
codic <- sample(0:5,1,replace=F)
syssample <- pop[codpop == codic]
length(syssample)
syssample[1:100]
rm(list=ls())

F <- c(1, 1)
n <- 2
while (F[n] <= 100) {
n <- n + 1
F[n] <- F[n-1] + F[n-2]
}
cat("Il primo numero di Fibonacci > 100 è F(", n, ") =", F[n], "\n")

F <- c(1, 1)
for (n in 3:100) {
F[n] <- F[n-1] + F[n-2]
}
cat("Il centesimo numero di Fibonacci è F(", 100, ") =", F[100], "\n")