## Task 1.2
for (k in 3:N) {
  Fn[k] <- Fn[k-1] + Fn[k-2]
  Pn[k] <- 2*Pn[k-1] + Pn[k-2]
}

data.frame(n=n[1:10], Fn=Fn[1:10], Pn=Pn[1:10])