Problemas Método de bisección. Método de Newton. - CAA EII

Problemas Método de bisección. Método de Newton. - CAA EII Problemas Método de bisección. Método de Newton. - CAA EII

10.05.2013 Views

De hecho, cualquiera que sea la cantidad inicial A, el máximo de la concentración se alcanza en t = 3. Para que la cantidad máxima no supere 1 mg/ml, simplemente hay que poner un A menor del que resuelva 3e −1 A = 1, esto es, A ≤ e/3. b) Supuesta suministrada esta cantidad inicial, A = e/3, nos preguntan cuándo ocurrirá que c(t) = e 3 te−t/3 = 0 ′ 25. Tenemos que resolver pues la ecuación f(t) = e 3 te−t/3 − 0 ′ 25 = 0. Pero del análisis anterior sacamos que al máximo llega desde cero y que después tiende a cero, es decir, cualquier cantidad en (0, A), en particular 0 ′ 25, es alcanzada dos veces y nos interesa cuando lo alcanza por segunda vez. 1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 0 5 10 15 c(t) = e 3 te−t/3 . Usamos el método de Newton-Raphson con algún valor numérico que resulte adecuado (por lo dicho antes, más bien largo que corto, para que la convergencia sea hacia la solución por encima de las tres horas). Por la forma de campana, basta tomar cualquier valor inicial superior a tres horas, p. ej. x0 = 4 genera una sucesión que converge en 4 iteraciones (con una tolerancia de 10 −8 ; nos pedían precisición de un minuto, que pasado a horas es 0 ′ 01666, por tanto bastaba con 10 −3 ): x1 = 10 ′ 30879920381349 x2 = 11 ′ 02140006874606 x3 = 11 ′ 07757133060878 x4 = 11 ′ 07790357510362 horas= 11 horas, 4 minutos. 17. En este problema nos dan el modelo N(t) = N0e λt + v e λ λt − 1 , los datos N0 = 10 6 , v = 4 10 5 , N(1) = 1506 10 3 , y nos piden simplemente que despejemos el valor de λ : Ingeniería Técnica Forestal 10 Fundamentos Matemáticos Curso 2004/05

1506 10 3 = 10 6 e λ 4 105 + e λ λ − 1 . Simplificamos antes de resolver el cero de cierta función asociada: f(λ) = 1000e λ + 400 λ e λ − 1 − 1506 = 0. Con un dato inicial bajo x1 = 1 (lo esperable para la natalidad de una población que en un año sólo ha pasado de 1400000 a algo más de millón y medio de habitantes) la respuesta se obtiene en 4 iteraciones (con precisión 10 −8 ): x2 = 0 ′ 390820116865 x3 = 0 ′ 12534036450639 x4 = 0 ′ 08560613200322 x5 = 0 ′ 08483943705533 x6 = 0 ′ 08483915873218 ∼ λ. b) Consiste en sustituir t = 3 en la expresión (ahora totalmente conocida) de N(t) : N(3) = 2656373 ′ 589004676, que significa una población en torno a los 2.656373 habitantes. Ingeniería Técnica Forestal 11 Fundamentos Matemáticos Curso 2004/05

De hecho, cualquiera que sea la cantidad inicial A, el máximo <strong>de</strong> la concentración se<br />

alcanza en t = 3. Para que la cantidad máxima no supere 1 mg/ml, simplemente hay que<br />

poner un A menor <strong>de</strong>l que resuelva 3e −1 A = 1, esto es, A ≤ e/3.<br />

b) Supuesta suministrada esta cantidad inicial, A = e/3, nos preguntan cuándo ocurrirá<br />

que<br />

c(t) = e<br />

3 te−t/3 = 0 ′ 25.<br />

Tenemos que resolver pues la ecuación<br />

f(t) = e<br />

3 te−t/3 − 0 ′ 25 = 0.<br />

Pero <strong>de</strong>l análisis anterior sacamos que al máximo llega <strong>de</strong>s<strong>de</strong> cero y que <strong>de</strong>spués tien<strong>de</strong> a<br />

cero, es <strong>de</strong>cir, cualquier cantidad en (0, A), en particular 0 ′ 25, es alcanzada dos veces y nos<br />

interesa cuando lo alcanza por segunda vez.<br />

1<br />

0.9<br />

0.8<br />

0.7<br />

0.6<br />

0.5<br />

0.4<br />

0.3<br />

0.2<br />

0.1<br />

0<br />

0 5 10 15<br />

c(t) = e<br />

3 te−t/3 .<br />

Usamos el método <strong>de</strong> <strong>Newton</strong>-Raphson con algún valor numérico que resulte a<strong>de</strong>cuado<br />

(por lo dicho antes, más bien largo que corto, para que la convergencia sea hacia la solución<br />

por encima <strong>de</strong> las tres horas). Por la forma <strong>de</strong> campana, basta tomar cualquier valor inicial<br />

superior a tres horas, p. ej. x0 = 4 genera una sucesión que converge en 4 iteraciones (con<br />

una tolerancia <strong>de</strong> 10 −8 ; nos pedían precisición <strong>de</strong> un minuto, que pasado a horas es 0 ′ 01666,<br />

por tanto bastaba con 10 −3 ):<br />

x1 = 10 ′ 30879920381349<br />

x2 = 11 ′ 02140006874606<br />

x3 = 11 ′ 07757133060878<br />

x4 = 11 ′ 07790357510362 horas= 11 horas, 4 minutos.<br />

17. En este problema nos dan el mo<strong>de</strong>lo<br />

N(t) = N0e λt + v<br />

<br />

e<br />

λ<br />

λt <br />

− 1 ,<br />

los datos N0 = 10 6 , v = 4 10 5 , N(1) = 1506 10 3 , y nos pi<strong>de</strong>n simplemente que <strong>de</strong>spejemos el<br />

valor <strong>de</strong> λ :<br />

Ingeniería Técnica<br />

Forestal<br />

10 Fundamentos Matemáticos<br />

Curso 2004/05

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

Saved successfully!

Ooh no, something went wrong!