1 <?xml version="1.0" encoding="UTF-8"?>
4 * Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
5 * Copyright (C) 2008 - INRIA
8 * Copyright (C) 2012 - 2016 - Scilab Enterprises
10 * This file is hereby licensed under the terms of the GNU GPL v2.0,
11 * pursuant to article 5.3.4 of the CeCILL v.2.1.
12 * This file was originally licensed under the terms of the CeCILL v2.1,
13 * and continues to be available under such terms.
14 * For more information, see the COPYING file which you should have received
15 * along with this program.
19 <refentry xmlns="http://docbook.org/ns/docbook" xmlns:xlink="http://www.w3.org/1999/xlink"
20 xmlns:svg="http://www.w3.org/2000/svg" xmlns:ns5="http://www.w3.org/1999/xhtml"
21 xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:db="http://docbook.org/ns/docbook"
22 xmlns:scilab="http://www.scilab.org" xml:id="int2d" xml:lang="ja">
25 <refname>int2d</refname>
27 <refpurpose>求積法および立体求積法により2次元積分を定義</refpurpose>
36 [I, err] = int2d(X, Y, f)
37 [I, err] = int2d(X, Y, f, params)
56 N個の三角形の頂点の横座標軸の値を有する3 行 N 列の配列.
72 N個の三角形の頂点の縦座標軸の値を有する3 行 N 列の配列.
88 被積分関数 <literal>f(u,v)</literal>を定義する
106 実数ベクトル <literal>[tol, iclose, maxtri, mevals,
112 .デフォルト値は <literal>[1.d-10, 1, 50, 4000,
134 <literal>iflag=0</literal>の場合, <literal>tol</literal> は
138 <literal>iflag=</literal>1の場合,
156 LQM0 または LQM1の選択を定義する積分パラメータ.
158 <literal>iclose=1</literal>の場合, LQM1 が使用されます.
160 <literal>iclose</literal>がその他の値の場合,
164 LQM0 は三角形の内部の点における値のみ関数の値を使用します.
166 LQM1 は通常 LQM0 より正確ですが, 三角形の境界点を含む
170 通常,被積分関数が三角形の境界において特異とならない限り,
188 領域の最後の三角形分割における三角形の最大数
208 LQM1の場合には94*<literal>maxtri</literal>未満,
210 LQM0の場合には56*<literal>maxtri</literal>未満
212 とした場合にのみ計算負荷を制限する効果があります.
272 <literal>int2d</literal> は,
274 <literal>n</literal>個の三角形からなる領域において
276 関数<literal>f</literal>の2次元積分を計算します.
278 全体の推定誤差が得られ,このサブルーチンへの入力として指定された
280 許容誤差 <literal>tol</literal>と比較されます.
282 許容誤差は<literal>iflag</literal>の入力値に基づき
288 各入力三角形に適用され,積分の合計値が推定され,
292 ローカルな求積モジュールは,サブルーチン LQM0 または
296 入力変数<literal>iclose</literal>の値によりどちらを
306 絶対誤差が最大となる三角形は最も長い辺の中点で
310 この後,ローカルな求積モジュールが分割された各三角形に
312 適用され積分値と誤差の新たな推定値が得られます.
314 この処理は,(1)許容誤差の制約が満たされるか,
316 (2)生成される三角形の数が入力パラメータ<literal>maxtri</literal>
318 を超えるか,(3)被積分関数を評価した数が入力パラメータ
320 <literal>mevals</literal>を超えるか,もしくは(4)
322 丸め誤差が結果を汚染し始めたとこの関数が検知した場合,
334 <programlisting role="example"><![CDATA[
337 deff('z=f(x,y)','z=cos(x+y)')
338 [I,e] = int2d(X, Y, f)
339 // 正方形の区間 [0 1]x[0 1] で被積分関数を計算します
344 <refsection role="see also">
348 <simplelist type="inline">
350 <link linkend="mesh2d">mesh2d</link>
353 <link linkend="int3d">int3d</link>
356 <link linkend="intg">intg</link>
359 <link linkend="intl">intl</link>
362 <link linkend="intc">intc</link>