『実例で学ぶ精度保証付き計算』サポートページ

『実例で学ぶ精度保証付き計算ー理論と実装ー』
臨時別冊・数理科学2011年10月, SCGライブラリ-85, サイエンス社,
(B5判; 216頁 発行日: 2011-10-25)
のプログラムコード群，正誤表，補遺 などのサポートページです．

• プログラムコード群
  • 2011年10月19日版 (zipファイル; "FOR_WEB_2011-10-19.zip"; 49KB)
  • フォルダ(ディレクトリ)構成; 本文に記述したコードに対応しています．
    1. intro: 第1章
    2. Emden: 第2章(2.10節)
    3. FS-Int: 第5章
    4. ode: 第6章
    5. pde: 第7章

• 正誤表
  1. p.83; 8行目; 埋め込みの向きが反対です。
  2. 局所一意性を検証するINTLABコードに誤りがありました。 以下の通り、修正をお願いします。 御指摘いただいた、早稲田大学 関根 晃太 様に深く感謝いたします。
     1. p.138 IN_Linz.m (4/4)

        【誤】

         s = alpha;  icon = 0;  it = 0;
         while icon == 0;
           s = s*1.01;
           s3 = (s+C*t1)*(3/two*(s+C*t1)+sqrt(two)*t3);
           Cg = Cp/epsilon^2*(s3+mag(delta)*C^2);  it = it + 1;
           if sup(Cg*M) >= 1, icon = 1; end
         end
         if it ~= 1, alphaU = s; end
        

        【正】

         si = alpha;  icon = 0;  it = 0;
         while icon == 0;
           s = si*1.01;
           s3 = (s+C*t1)*(3/two*(s+C*t1)+sqrt(two)*t3);
           Cg = Cp/epsilon^2*(s3+mag(delta)*C^2);  it = it + 1;
           if sup(Cg*M) >= 1
             icon = 1; 
           else
             si = s;
           end
         end
         if it ~= 1, alphaU = si; end
        

     2. p.165 IS_Res.m

        【誤】

              t = alpha; icon = 0;  it = 0;
              while icon == 0
                t = t*1.1;
                s = 2*t3*(Er(4)*t+Ei(4)*C*t1) + ...
                    3*a*lambda*(Er(8)*t+Ei(8)*C*t1)^2;
                Cg = Cp*t2+Er(4)*s;  it = it + 1;
                if sup(Cg*Cp) >= 1, icon = 1; end
              end
              if it ~= 1, alphaU = t; end
        

        【正】

              ti = alpha; icon = 0;  it = 0;
              while icon == 0
                t = ti*1.1;
                s = 2*t3*(Er(4)*t+Ei(4)*C*t1) + ...
                    3*a*lambda*(Er(8)*t+Ei(8)*C*t1)^2;
                Cg = Cp*t2+Er(4)*s;  it = it + 1;
                if sup(Cg*Cp) >= 1
                  icon = 1; 
                else 
                  ti = t;
                end
              end
              if it ~= 1, alphaU = ti; end
        

     3. p.168 IN_Linz.m (3/3)

        【誤】

              t = alpha; icon = 0;  it = 0;
              while icon == 0
                t = t*1.1;
                s = 3*a*lambda*(Er(8)*t+Ei(8)*C*t1)^2 + ...
                    2*t3*(Er(4)*t+Ei(4)*C*t1);
                Cg = Er(4)*s;  it = it + 1;
                if sup(Cg*M) >= 1, icon = 1; end
              end
              if it ~= 1, alphaU = t; end
        

        【正】

              ti = alpha; icon = 0;  it = 0;
              while icon == 0
                t = ti*1.1;
                s = 3*a*lambda*(Er(8)*t+Ei(8)*C*t1)^2 + ...
                    2*t3*(Er(4)*t+Ei(4)*C*t1);
                Cg = Er(4)*s;  it = it + 1;
                if sup(Cg*M) >= 1
                  icon = 1;
                else
                  ti = t;
                end
              end
              if it ~= 1, alphaU = ti; end