Mais uma dúvida simples: O problema é no final do código ( em vermelho) : datas <- read.table("ddd.txt", header=F) datas <- as.Date(datas[,1], format="%m/%d/%y") i <- read.table("iii.txt", header=F) i <- i[,1] cupons <- read.table("cc1.txt", header=F) cupons <- as.Date(cupons[,1], format="%m/%d/%Y") fluxo <- read.table("ff1.txt", header=F) fluxo <- fluxo[,1] falta <- outer (cupons, datas, "-") *.69 falta <- matrix(falta,nrow(falta),ncol(falta)) falta <- ifelse(falta < 0, NA, falta) VP <- cbind((fluxo)/((1+i)^(falta))) plot(colSums(VP, na.rm=T), type="l") i tem 297 elementos, datas tem 297 elementos cupons tem 14 elementos, e fluxo tem 14 elementos; Logo, falta e VP tem 297 colunas e 14 linhas quanto ploto VP, aparecem 14 vales no gráfico. Acho que algo de errado em (1+i)^(falta). O que quero fazer é,(1+i[1])^(falta[,1]), (1+i[2])^(falta[,2]), ... , (1+i[n])^(falta[,n]), ou seja, elevar (1+ i), para cada i["n"] em todos os elementos da coluna "n" em "falta". Assim, acabaria também com uma matriz 14 por 297. E quero dividir os 14 elementos de "fluxo" pelos 14 elementos de cada coluna na matriz gerada. Deu pra entender? Obrigado, Luciano.