Cap´ıtulo 21 Introduç˜ao `a Integral: Cálculo de´Areas e Integrais ...
Cap´ıtulo 21 Introduç˜ao `a Integral: Cálculo de´Areas e Integrais ...
Cap´ıtulo 21 Introduç˜ao `a Integral: Cálculo de´Areas e Integrais ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
W.Bianchini, A.R.Santos 295<br />
> numeros:=[seq(rand(),k=1..31)]; k:=’k’:<br />
numeros := [97396414947, 780422731613, 987785640265, 674198272844,<br />
134050365811, 754869582636, 140810856859, 347877704841, 433599229456,<br />
898724880795, 485531802023, 255050614524, 952922474293, 642065329619,<br />
154912668026, 856069438450, 681407641506, 962917791070, 874166946435,<br />
905950292905, 549552888716, 84125842236, 67060541266, 6<strong>21</strong>757734462,<br />
223575905687, 273574099511, 410424381304, 659501247275, 887974857856,<br />
234450269247, 606386273485]<br />
Como queremos pontos pertencentes ao intervalo [0, 1], vamos converter os pontos gerados pelo comando acima<br />
para este intervalo, por uma mudança de escala:<br />
> pts:=map(x->evalf(x/10^12),numeros);<br />
pts := [.09739641495, .7804227316, .9877856403, .6741982728, .1340503658,<br />
.7548695826, .1408108569, .3478777048, .4335992295, .8987248808,<br />
.4855318020, .2550506145, .9529224743, .6420653296, .1549126680,<br />
.8560694385, .6814076415, .9629177911, .8741669464, .9059502929,<br />
.5495528887, .08412584224, .06706054127, .6<strong>21</strong>7577345, .2235759057,<br />
.2735740995, .4104243813, .6595012473, .8879748579, .2344502692,<br />
.6063862735]<br />
Para formar os pontos da partição, precisamos colocar esta última seqüência em ordem crescente. Isto é feito<br />
utilizando-se o comando sort:<br />
> part:=sort(pts);<br />
part := [.06706054127, .08412584224, .09739641495, .1340503658, .1408108569,<br />
.1549126680, .2235759057, .2344502692, .2550506145, .2735740995,<br />
.3478777048, .4104243813, .4335992295, .4855318020, .5495528887,<br />
.6063862735, .6<strong>21</strong>7577345, .6420653296, .6595012473, .6741982728,<br />
.6814076415, .7548695826, .7804227316, .8560694385, .8741669464,<br />
.8879748579, .8987248808, .9059502929, .9529224743, .9629177911,<br />
.9877856403]<br />
Podemos agora, calcular a soma de Riemann associada a esta partição do intervalo [0, 1], como se segue.<br />
> f:=x->x^3+3*x^2+2*x-5;<br />
f := x → x 3 + 3 x 2 + 2 x − 5<br />
> S:=sum(f(part[2*j])*(part[2*j+1]-part[2*j-1]),j=1..15);<br />
S := −2.403062293<br />
1. Repita este processo mais cinco vezes e guarde os resultados. Calcule a média das suas 6 tentativas e descreva<br />
como este processo forma uma soma de Riemann geral e como por meio dele se chega a uma aproximação do<br />
valor da integral da função no intervalo [0, 1]. Ilustre geometricamente.<br />
2. Explique como é possível melhorar a precisão do resultado e aplique as suas conclusões para melhorar o resultado<br />
obtido acima.<br />
3. Ache por este processo uma aproximação para a integral da função f(x) = x 3 + x + 2 no intervalo [0, 1].<br />
<strong>21</strong>.10.2 Somas de Riemann e funções monótonas<br />
O objetivo deste projeto é calcular integrais de funções monótonas por meio de somas de Riemann com um erro<br />
máximo prefixado.<br />
1. Considere a função f(x) = x 3 + x + 2.<br />
(a) Mostre que f é monótona no intervalo [0, 2].