Page 258 - riesgo2012

248 8. Teor´ ıa de la ruina: tiempo continuo

f<- array(1:R, dim=c(R))
f[1] <- 0.1
f[2] <- 0.1
f[3] <- 0.2
f[4] <- 0.3
f[5] <- 0.3
for (i in 5:R) { f[i] <- 0 }
#................................
#C\’alculodela densidaddeS
#................................
g0 <- p0
g<- array(1:R, dim=c(R))
g[1] <- (a+b)*f[1]*g0
for (r in 2: R) {
aux <- 0
for (i in 1:{r-1}) {
aux <- aux + (a+b*i/r)*f[i]*g[r-i]
}
aux <- aux + (a+b)*f[r]*g0
g[r] <- aux
}
#...............................
#Graficaci\’on
#...............................
#Nota:Se omiteen lagr\’aficaelvalorde ladensidaden cero "g0".
barplot(g,main="Funcin de densidad de S",xlab="r", ylab="g(r)")
#
#################################################################
#Fin dec\’odigo
#################################################################