r - Rで製品関数をプロットする方法は?

Question

私は問題を試していますが、最初にプロットする必要がありますr

(x+1)(x+2)...(x+n),

n固定整数です。

このルーチンを作成する方法はありますか?

Answer 1

が-1xより大きい場合、これは関係を利用することによって最も効率的に計算される可能性があります。

(x + 1)*(x + 2)* ... *(x + n) = Gamma(x+n+1) / Gamma(x+1).

ガンマは対数で内部的に計算されるため、これらのログを次の形式で使用しますlgamma。

f <- function(x, n) exp(lgamma(x+n+1) - lgamma(x+1))

次に、たとえば次のように を介して プロットを取得できます。curve

curve(f(x,3), 0, pi)

Answer 2

このようなものが欲しいですか？

f <- function(x, n) {
  return(prod(1/(x+(1:n))))
}