関数 Sn=LeftRiemannSum(f,left,right,N) を作成したとします。これは、N 個のサブドメインで左から右への間隔で左リーマン和を計算します。
すべてのサブドメイン i=0 から N-1 のSn=sum(f(xi) h)。f は私の関数であり、xi=left+i h であるため、入力引数は left=x0 および right=xN です。
f を無名関数 (例 f=@(x)(x.*log(1+x)) ) とします。
また、N を 10 から 100000 まで変化させた場合の Sn も推定しました。
ここで、N -> 無限大の場合の系列 Sn の値を計算したいだけです。関数内に for ループ [ for i=0:(N-1)] があるので、無限ループになります ..
関数を何らかの方法で「limit」コマンドに渡すことはできますか? 手がかりはありますか？
ありがとう！
PS: 私の関数 LeftRiemann Sum のコードの主要部分は次のとおりです。
for i=0:(N-1)
x=x0+i.*h;
y=f(x);
A=y.*h;
S=S+A
end
Sn=S