написал эту функцию

ifun <- function(m)  {
  o = c() 
  for (k in 1:m) {
    o[k] = prod(1:k) / prod(2 * (1:k) + 1)
    }
  o_sum = 2 * (1 + sum(o))  # Final result

  print(o_sum)
}

Эта функция приближает постояннуюpiОднако послеm > 50 аппроксимация застревает, то есть аппроксимация имеет то же значение и не улучшается. Как я могу это исправить? Благодарю.