13.05.2013 Views

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 ...

SHOW MORE
SHOW LESS

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].

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!