14
14
<refentry xmlns="http://docbook.org/ns/docbook" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:svg="http://www.w3.org/2000/svg" xmlns:ns5="http://www.w3.org/1999/xhtml" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:db="http://docbook.org/ns/docbook" version="5.0-subset Scilab" xml:id="int2d" xml:lang="en">
16
<refname>int2d</refname>
17
<refpurpose>integral definida 2d por quadratura e cubatura</refpurpose>
20
<title>Seq��ncia de Chamamento</title>
21
<synopsis>[I,err]=int2d(X,Y,f [,params])</synopsis>
24
<title>Par�metros</title>
29
<para>um array 3 por N contendo as abscissas dos vertices dos N
37
<para>um array 3 por N contendo as ordenadas dos vertices dos N
45
<para>fun��o externa (fun��o, string ou lista) definindo o
46
integrando <literal>f(u,v)</literal>;
54
vetor de reais <literal>[tol, iclose, maxtri, mevals,
57
.O valor padr�o � <literal>[1.d-10, 1, 50, 4000,
66
<para>o limite desejado do erro. Se
67
<literal>iflag=0</literal>, <literal>tol</literal> �
68
interpretado como um limite de erro relativo; se
69
<literal>iflag=</literal>1, o limite � de erro
77
<para>um inteiro que determina a sele��o dos m�todos LQM0 ou
78
LQM. Se <literal>iclose=1</literal> , ent�o LQM1 � utilizado.
79
Qualquer outro valor de <literal>iclose</literal> faz com que
80
LQM0 seja usado. LQM0 utiliza valores da fun��o apenas em
81
pontos interiores ao tri�ngulo. LQM1 geralmente � mais preciso
82
que LQM0 mas envolve a avalia��o do integrando em mais pontos,
83
incluindo em alguns pontos da fronteira do tri�ngulo.
84
Geralmente � melhor utilizar LQM1 a n�o ser que o integrando
85
possuia singularidades nas bordas do tri�ngulo.
92
<para>o n�mero m�ximo de tri�ngulos na triangulariza��o final
100
<para>o n�mero m�ximo de avalia��es da fun��o permitido. Este
101
n�mero ter� efeito na limita��o da computa��o se for menor que
102
94*<literal>maxtri</literal> quando LQM1 � especificado ou
103
56*<literal>maxtri</literal> quando LQM0 �
120
<para>o valor da integral</para>
126
<para>o erro estimado</para>
132
<title>Descri��o</title>
134
<literal>int2d</literal> computa a integral bidimensional de uma
135
fun��o <literal>f</literal> sobre uma regi�o que consiste de
136
<literal>n</literal> tri�ngulos. Um estimativa de erro total � obtida e
137
comparada a - <literal>tol</literal> - que � fornecida como entrada para a
138
subrotina. A toler�ncia de erro � tratada como relativa ou absoluta
139
dependendo do valor de entrada de <literal>iflag</literal>. Um "m�dulo de
140
quadratura local" ("Local Quadrature Module") � aplicado para cada
141
tri�ngulo de entrada e estimativas da integral total e do erro total s�o
142
computadas. O m�dulo de quadratura local � a subrotina LQM0 ou a subrotina
143
LQM1 e a escolha entre elas � determinada pelo valor da vari�vel
144
<literal>iclose</literal>.
146
<para>Se a estimativa de erro total excede a toler�ncia, o tri�ngulo com
147
maior erro absoluto � dividio em dois outro tri�ngulos tra�ando-se uma
148
mediana por seu maior lado. O m�dulo de quadratura local � ent�o aplicado
149
a cada um dos subtri�ngulos para se obter novas estimativas da integral e
150
do erro. Este processo � repetido at� que um dos seguintes (1) a
151
toler�ncia � satisfeita, (2) o n�mero de tri�ngulos gerados excede o
152
par�metro <literal>maxtri</literal>, (3) o n�mero de avalia��es do
153
integrando excede o par�metro <literal>mevals</literal>, ou (4) a fun��o
154
sente que um erro de arredondamento est� come�ando a contaminar o
159
<title>Exemplos</title>
160
<programlisting role="example"><![CDATA[
16
<refname>int2d</refname>
17
<refpurpose>integral definida 2d por quadratura e cubatura</refpurpose>
20
<title>Seq��ncia de Chamamento</title>
21
<synopsis>[I,err]=int2d(X,Y,f [,params])</synopsis>
24
<title>Par�metros</title>
29
<para>um array 3 por N contendo as abscissas dos vertices dos N
37
<para>um array 3 por N contendo as ordenadas dos vertices dos N
45
<para>fun��o externa (fun��o, string ou lista) definindo o
46
integrando <literal>f(u,v)</literal>;
54
vetor de reais <literal>[tol, iclose, maxtri, mevals,
57
.O valor padr�o � <literal>[1.d-10, 1, 50, 4000,
66
<para>o limite desejado do erro. Se
67
<literal>iflag=0</literal>, <literal>tol</literal> �
68
interpretado como um limite de erro relativo; se
69
<literal>iflag=</literal>1, o limite � de erro
77
<para>um inteiro que determina a sele��o dos m�todos LQM0 ou
78
LQM. Se <literal>iclose=1</literal> , ent�o LQM1 � utilizado.
79
Qualquer outro valor de <literal>iclose</literal> faz com que
80
LQM0 seja usado. LQM0 utiliza valores da fun��o apenas em
81
pontos interiores ao tri�ngulo. LQM1 geralmente � mais preciso
82
que LQM0 mas envolve a avalia��o do integrando em mais pontos,
83
incluindo em alguns pontos da fronteira do tri�ngulo.
84
Geralmente � melhor utilizar LQM1 a n�o ser que o integrando
85
possuia singularidades nas bordas do tri�ngulo.
92
<para>o n�mero m�ximo de tri�ngulos na triangulariza��o final
100
<para>o n�mero m�ximo de avalia��es da fun��o permitido. Este
101
n�mero ter� efeito na limita��o da computa��o se for menor que
102
94*<literal>maxtri</literal> quando LQM1 � especificado ou
103
56*<literal>maxtri</literal> quando LQM0 �
120
<para>o valor da integral</para>
126
<para>o erro estimado</para>
132
<title>Descri��o</title>
134
<literal>int2d</literal> computa a integral bidimensional de uma
135
fun��o <literal>f</literal> sobre uma regi�o que consiste de
136
<literal>n</literal> tri�ngulos. Um estimativa de erro total � obtida e
137
comparada a - <literal>tol</literal> - que � fornecida como entrada para a
138
subrotina. A toler�ncia de erro � tratada como relativa ou absoluta
139
dependendo do valor de entrada de <literal>iflag</literal>. Um "m�dulo de
140
quadratura local" ("Local Quadrature Module") � aplicado para cada
141
tri�ngulo de entrada e estimativas da integral total e do erro total s�o
142
computadas. O m�dulo de quadratura local � a subrotina LQM0 ou a subrotina
143
LQM1 e a escolha entre elas � determinada pelo valor da vari�vel
144
<literal>iclose</literal>.
146
<para>Se a estimativa de erro total excede a toler�ncia, o tri�ngulo com
147
maior erro absoluto � dividio em dois outro tri�ngulos tra�ando-se uma
148
mediana por seu maior lado. O m�dulo de quadratura local � ent�o aplicado
149
a cada um dos subtri�ngulos para se obter novas estimativas da integral e
150
do erro. Este processo � repetido at� que um dos seguintes (1) a
151
toler�ncia � satisfeita, (2) o n�mero de tri�ngulos gerados excede o
152
par�metro <literal>maxtri</literal>, (3) o n�mero de avalia��es do
153
integrando excede o par�metro <literal>mevals</literal>, ou (4) a fun��o
154
sente que um erro de arredondamento est� come�ando a contaminar o
159
<title>Exemplos</title>
160
<programlisting role="example"><![CDATA[
163
163
deff('z=f(x,y)','z=cos(x+y)')
164
164
[I,e]=int2d(X,Y,f)
165
165
// computa o integrando sobre o quadrado [0 1]x[0 1]
166
166
]]></programlisting>
168
<refsection role="see also">
169
<title>Ver Tamb�m</title>
170
<simplelist type="inline">
172
<link linkend="intc">intc</link>
175
<link linkend="intl">intl</link>
178
<link linkend="int3d">int3d</link>
181
<link linkend="intg">intg</link>
184
<link linkend="mesh">mesh</link>
168
<refsection role="see also">
169
<title>Ver Tamb�m</title>
170
<simplelist type="inline">
172
<link linkend="intc">intc</link>
175
<link linkend="intl">intl</link>
178
<link linkend="int3d">int3d</link>
181
<link linkend="intg">intg</link>
184
<link linkend="mesh">mesh</link>