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