1
1
<?xml version="1.0" encoding="ISO-8859-1"?>
2
2
<refentry xmlns="http://docbook.org/ns/docbook" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:svg="http://www.w3.org/2000/svg" xmlns:ns3="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="int3d" xml:lang="en">
4
<refname>int3d</refname>
5
<refpurpose>integral definida 3d pelo m�todo da quadratura e
10
<title>Seq��ncia de Chamamento</title>
11
<synopsis>[result,err]=int3d(X,Y,Z,f [,nf[,params]])</synopsis>
14
<title>Par�metros</title>
20
um array 4 por <literal>NUMTET</literal> contendo as abscissas
21
dos v�rtices dos <literal>NUMTET</literal> tetraedros.
29
um array 4 por <literal>NUMTET</literal> contendo as ordenadas
30
dos v�rtices dos <literal>NUMTET</literal> tetraedros.
38
um array 4 por <literal>NUMTET</literal> contendo as terceiras
39
coordenadas dos v�rtices dos <literal>NUMTET</literal>
47
<para>fun��o externa (function, string ou lista) definindo o
48
integrando <literal>f(xyz,nf)</literal>, onde <literal>xyz</literal>
49
� o vetor das coordenadas de um ponto e nf os n�meros das
57
<para>o n�mero da fun��o a ser integrada (o padr�o � 1)</para>
64
vetor de reais <literal>[minpts, maxpts, epsabs,
67
.O valor padr�o � <literal>[0, 1000, 0.0,
76
<para>limite do erro absoluto desejado</para>
82
<para>limite do erro relativo desejado</para>
88
<para>n�mero m�nimo de avalia��es da fun��o</para>
94
<para>n�mero m�ximo de avalia��es da fun��o. O n�mero de
95
avalia��es da fun��o sobre cada subregi�o � 43
105
<para>o valor da integral, ou vetor de valores da integral.</para>
111
<para>estimativas do erro absoluto</para>
117
<title>Descri��o</title>
118
<para>A fun��o calcula uma aproxima��o a um dado vetor de integrais
121
<programlisting role="no-scilab-exec"><![CDATA[
4
<refname>int3d</refname>
5
<refpurpose>integral definida 3d pelo m�todo da quadratura e
10
<title>Seq��ncia de Chamamento</title>
11
<synopsis>[result,err]=int3d(X,Y,Z,f [,nf[,params]])</synopsis>
14
<title>Par�metros</title>
20
um array 4 por <literal>NUMTET</literal> contendo as abscissas
21
dos v�rtices dos <literal>NUMTET</literal> tetraedros.
29
um array 4 por <literal>NUMTET</literal> contendo as ordenadas
30
dos v�rtices dos <literal>NUMTET</literal> tetraedros.
38
um array 4 por <literal>NUMTET</literal> contendo as terceiras
39
coordenadas dos v�rtices dos <literal>NUMTET</literal>
47
<para>fun��o externa (function, string ou lista) definindo o
48
integrando <literal>f(xyz,nf)</literal>, onde <literal>xyz</literal>
49
� o vetor das coordenadas de um ponto e nf os n�meros das
57
<para>o n�mero da fun��o a ser integrada (o padr�o � 1)</para>
64
vetor de reais <literal>[minpts, maxpts, epsabs,
67
.O valor padr�o � <literal>[0, 1000, 0.0,
76
<para>limite do erro absoluto desejado</para>
82
<para>limite do erro relativo desejado</para>
88
<para>n�mero m�nimo de avalia��es da fun��o</para>
94
<para>n�mero m�ximo de avalia��es da fun��o. O n�mero de
95
avalia��es da fun��o sobre cada subregi�o � 43
105
<para>o valor da integral, ou vetor de valores da integral.</para>
111
<para>estimativas do erro absoluto</para>
117
<title>Descri��o</title>
118
<para>A fun��o calcula uma aproxima��o a um dado vetor de integrais
121
<programlisting role="no-scilab-exec"><![CDATA[
122
122
I I I (F ,F ,...,F ) dx(3)dx(2)dx(1),
124
124
]]></programlisting>
125
<para>onde a regi�o de integra��o s�o os NUMTET tetraedros e onde</para>
126
<programlisting role="no-scilab-exec"><![CDATA[
125
<para>onde a regi�o de integra��o s�o os NUMTET tetraedros e onde</para>
126
<programlisting role="no-scilab-exec"><![CDATA[
127
127
F = F (X(1),X(2),X(3)), J = 1,2,...,NUMFUN.
129
129
]]></programlisting>
130
<para>uma estrat�gia globalmente adaptativa � aplicada para se computar
131
aproxima��es <literal>result(k)</literal> esperando-se que satisfa�a, para
132
cada componente de I, � seguinte precis�o:
133
<literal>ABS(I(K)-RESULT(K))<=MAX(EPSABS,EPSREL*ABS(I(K)))</literal>
136
<literal>int3d</literal> repetidamente subdivide os tetraedros com
137
maiores erros estimados e estima as integrais e os erros sobre os novos
138
tetraedros at� que a exig�ncia de erro seja encontrada ou
139
<literal>MAXPTS</literal> avalia��es da fun��o tenham sido feitas.
141
<para>Uma regra de 43 pontos de integra��o com todos os pontos de
142
avalia��o dentro dos tetraedros � aplicada. A regra tem grau polinomial
146
Se os valores dos par�metros de entrada <literal>EPSABS</literal> ou
147
<literal>EPSREL</literal> s�o suficientemente grandes, uma regra de
148
integra��o � aplicada sobre cada tetraedro e os valores s�o aproximados
149
para se fornecer as aproxima��es <literal>RESULT(K)</literal>. Nenhuma
150
subdivis�o posterior dos tetraedros ser� feita.
153
Quando <literal>int3d</literal> computa estimativas a um vetor de
154
integrais, � dado tratamento igual a todos os componentes do vetor. Isto
155
�, <literal>I(Fj)</literal> e <literal>I(Fk)</literal> para<literal>
158
diferente de <literal>k</literal>, s�o estimados com a mesma
159
subdivis�o da regi�o de integra��o. Para integrais suficientemente
160
semelhantes, podemos economizar tempo aplicando <literal>int3d</literal> a
161
todos os integrando em uma chamada. Para integrais que variam
162
continuamente em fun��o de um par�metro, as estimativas produzidas por
163
<literal>int3d</literal> tamb�m ir�o variar continuamente quando a
164
subdivis�o � aplicada a todos os componentes. Este geralmente n�o ser� o
165
caso quando componentes diferentes s�o tratados separadamente.
167
<para>Por outro lado, este recurso deve ser utilizado com cautela quando
168
os componentes diferentes da integral requerem claramente subdivis�es
173
<title>Refer�ncias</title>
174
<para>Rotina FORTRAN dcutet.f</para>
177
<title>Exemplos</title>
178
<programlisting role="example"><![CDATA[
130
<para>uma estrat�gia globalmente adaptativa � aplicada para se computar
131
aproxima��es <literal>result(k)</literal> esperando-se que satisfa�a, para
132
cada componente de I, � seguinte precis�o:
133
<literal>ABS(I(K)-RESULT(K))<=MAX(EPSABS,EPSREL*ABS(I(K)))</literal>
136
<literal>int3d</literal> repetidamente subdivide os tetraedros com
137
maiores erros estimados e estima as integrais e os erros sobre os novos
138
tetraedros at� que a exig�ncia de erro seja encontrada ou
139
<literal>MAXPTS</literal> avalia��es da fun��o tenham sido feitas.
141
<para>Uma regra de 43 pontos de integra��o com todos os pontos de
142
avalia��o dentro dos tetraedros � aplicada. A regra tem grau polinomial
146
Se os valores dos par�metros de entrada <literal>EPSABS</literal> ou
147
<literal>EPSREL</literal> s�o suficientemente grandes, uma regra de
148
integra��o � aplicada sobre cada tetraedro e os valores s�o aproximados
149
para se fornecer as aproxima��es <literal>RESULT(K)</literal>. Nenhuma
150
subdivis�o posterior dos tetraedros ser� feita.
153
Quando <literal>int3d</literal> computa estimativas a um vetor de
154
integrais, � dado tratamento igual a todos os componentes do vetor. Isto
155
�, <literal>I(Fj)</literal> e <literal>I(Fk)</literal> para<literal>
158
diferente de <literal>k</literal>, s�o estimados com a mesma
159
subdivis�o da regi�o de integra��o. Para integrais suficientemente
160
semelhantes, podemos economizar tempo aplicando <literal>int3d</literal> a
161
todos os integrando em uma chamada. Para integrais que variam
162
continuamente em fun��o de um par�metro, as estimativas produzidas por
163
<literal>int3d</literal> tamb�m ir�o variar continuamente quando a
164
subdivis�o � aplicada a todos os componentes. Este geralmente n�o ser� o
165
caso quando componentes diferentes s�o tratados separadamente.
167
<para>Por outro lado, este recurso deve ser utilizado com cautela quando
168
os componentes diferentes da integral requerem claramente subdivis�es
173
<title>Refer�ncias</title>
174
<para>Rotina FORTRAN dcutet.f</para>
177
<title>Exemplos</title>
178
<programlisting role="example"><![CDATA[