-2

の収束率を計算して決定する MATLAB コードを記述します。

(exp(h)-(1+h+1/2*h^2))/h ‍‍‍‍‍‍ ‍‍‍‍‍‍ with h=1/2, 1/2^2,..., 1/2^10

私のコードは次のとおりです。

h0=(0.5)^i;
TOL=10^(-8);
N=10;
i=1;
flag=0;
table=zeros(30,1);
table(1)=h0

while i < N
    h=(exp(h0)-(1+h0+0.5*h0^2))/h0;
    table (i+1)=h;
    if abs(h-h0)< TOL
        flag=1;
        break;
    end
    i=i+1;
    h0=h;
end

if flag==1
    h
else
    error('failed');
end

私が受け取った答えはまったく意味がありません。助けてください。

4

1 に答える 1