## Task 1.4
for (k in 1:N) {
  if (is.na(app[k])) {
    dif[k] <- NA
  } else {
    dif[k] <- round(abs(app[k]-real_sqrt2), 5)
  }
}

for (k in 1:N) {
  if (is.na(app[k])) {
    acc[k] <- NA
  } else {
    d <- 0
    for (x in 0:15) {
      if (round(app[k],x) == round(real_sqrt2,x)) {
        d <- x
      } else {
        break
      }
    }
    acc[k] <- d
  }
}

show <- c(2,3,5,10,50,100,199)
data.frame(n=n[show], approx=round(app[show],10), diff_5dp=dif[show], accuracy=acc[show])