x=read.table(file="input_piano_ammortamento_27_02_09.txt",sep=";",header=F)
i <- x[1,1]
S <- x[1,2]
n <- x[1,3]


nu <- 1/(1+i)
a_n_i <- nu * ( 1 - nu^n )/(1 - nu)
R <- S / a_n_i


amm <- matrix(c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0),ncol=6,byrow=TRUE)
colnames(amm) <- c("  k","  R_k","  C_k","  I_k","  D_k","  E_k")


I <- 1
amm[I,1] <- I-1
amm[I,2] <- 0
amm[I,3] <- 0
amm[I,4] <- 0
amm[I,5] <- S
amm[I,6] <- 0


lunghezza_ciclo <- n+1

for (I in 2:lunghezza_ciclo)   {
amm[I,1] <- I - 1
amm[I,2] <- R
amm[I,3] <- R - amm[I-1,5] * i
amm[I,4] <- amm[I-1,5] * i
amm[I,5] <- amm[I-1,5] - amm[I,3]
amm[I,6] <- amm[I-1,6] + amm[I,3]
}


write.table(amm, file="output_piano_ammortamento_27_02_09.txt", row.names=FALSE, col.names=TRUE)