matlab - Newton-Raphson 反復法の実装

Question

Newton Raphson 反復を使用して後方オイラー スキームを実装しようとしています。各反復で、最初の推測を行い、残差を計算し、変化を解くことを理解しています。私の場合、変更は del w です。その後、値を w^m に追加し、次の m 反復で w の更新値を取得することを知っています。反復が進むにつれて解の収束を確認することを知っています。私が抱えている問題は、時間ステップ dt を t=0:Tmax/dt として実装する方法です。ここで、Tmax は 10 です。時間ステップがどのように行われるかについて混乱しています。その間、どんな助けも大歓迎です。ありがとうございました！

 while Rw(m)>10^-6      % Convergence condition
    drdw(m)=(1-2*dt+2*t(n+1)^2*w(m)*dt);
    Dw(m)=Rw(m)\drdw(m); %Inverse
    w(m+1)=w(m)+Dw(m);  %Newton method
    Rw(m+1)=(-(w(m)-v(1)-2*w(m)*dt+t(n+1)^2*w(m)^2*dt));   %New Residual value
    if Rw(m+1)>10^-6  %Check on the level of convergence
        m=m+1;
    else
        Rw=1;  % I was thinking I should make the Residual 1 for the next time step. 
        break

    end