Lecture2 <- function()
{
 Part1()

}

# =============================================================================================

Part1 <- function()
{
 Nyear <-20
 rs <- c(0.1,0.2,0.3)
 Ks <- c(1000,2000,3000)
 H <- c(0.05,0.1,0.15)
 TheTab <- matrix(0,nrow=length(H),ncol=1+length(rs)*length(Ks))
 Probs <- rep(1/(length(rs)*length(Ks)),length=length(rs)*length(Ks))
 B <- rep(0,length=Nyear+1)
 Ipnt <- 0
 for (r in rs)
  for (K in Ks)
   {
    Ipnt <- Ipnt + 1
    for (IH in 1:length(H))
     {
      B[1] <- 0.2*K
      for (IY in 1:Nyear)
       B[IY+1] <- B[IY] + r*B[IY]*(1-B[IY]/K) - H[IH]*B[IY]
      TheTab[IH,Ipnt] <- B[IY+1]/K*100  
     }
   }
  for (IH in 1:length(H))
   TheTab[IH,Ipnt+1] <- sum(TheTab[IH,1:Ipnt]*Probs)
  print(round(TheTab,1))




}

# =============================================================================================

Lecture2()