# b) Determine a função de fiabilidade do servidor "vsrv8" nos dias 1 de janeiro e 1 de
  # fevereiro 2017. Apresente e analise os respetivos gráficos.

# 1 janeiro 2017
fiab = function(t) {
  up = nrow(subset(vpnsessions, as.Date(data) == (as.Date("2017-01-01", format="%Y-%m-%d")) & duracao > 1 & servidor == 'vsrv8'))
  down = nrow(subset(vpnsessions, as.Date(data) == (as.Date("2017-01-01", format="%Y-%m-%d")) & duracao < 2 & servidor == 'vsrv8'))

  deltaT = (60 * 24) / (up + down)

  tempoOperacao = up * deltaT

  lambda = down / tempoOperacao

  # F(t) = e ^ - (0.001574t)

  return (exp(-lambda * t))
}
plot(fiab(1:1440), type='l', main="Função de fiabilidade do servidor vsrv8 no dia 1 de janeiro de 2017", ylab = "Fiabilidade (%)", xlab="Minutos do dia 1 de janeiro de 2017")


# 1 fevereiro 2017
fiab = function(t) {
  up = nrow(subset(vpnsessions, as.Date(data) == (as.Date("2017-02-01", format="%Y-%m-%d")) & duracao > 1 & servidor == 'vsrv8'))
  down = nrow(subset(vpnsessions, as.Date(data) == (as.Date("2017-02-01", format="%Y-%m-%d")) & duracao < 2 & servidor == 'vsrv8'))

  deltaT = (60 * 24) / (up + down)

  tempoOperacao = up * deltaT

  lambda = down / tempoOperacao

  # F(t) = e ^ -0.0016975t
  print(lambda)

  return (exp(-lambda * t))
}
plot(fiab(1:1440), type='l', main="Função de fiabilidade do servidor vsrv8 no dia 1 de fevereiro de 2017", ylab = "Fiabilidade (%)", xlab="Minutos do dia 1 de fevereiro de 2017")