Mathematics · Livro 2 · Grades 10–12

Matemática do ensino médio

Matemática do ensino médio · Grades 10–12

26Equações diferenciais

Uma equação diferencial relaciona uma função com suas derivadas. A física, a química, a biologia e a economia exprimem suas leis nessa forma: a taxa de resfriamento de um corpo, o decaimento de um núcleo radioativo, o crescimento de uma população são todos afirmações sobre yy'. Este capítulo resolve por completo as equações lineares de primeira ordem com coeficientes constantes.

26.1 A equação y=ayy' = ay

Definição 26.1 (Equação diferencial)

Uma equação diferencial é uma equação cuja incógnita é uma função yy e que envolve yy e suas derivadas. Uma solução em um intervalo II é uma função derivável que satisfaz a equação em todo ponto de II.

Teorema 26.2 (Soluções de y=ayy' = ay)

Seja aRa \in \R. As soluções em R\R da equação y=ayy' = ay são exatamente as funções

y(x)=Ceax,CR.y(x) = C\,\eu^{ax}, \qquad C \in \R .

Para cada par (x0,y0)(x_0, y_0) existe uma única solução com y(x0)=y0y(x_0) = y_0.

Demonstração. Cada y=Ceaxy = C\eu^{ax} satisfaz y=Caeax=ayy' = Ca\,\eu^{ax} = ay. Reciprocamente, seja yy uma solução qualquer e ponha z(x)=y(x)eaxz(x) = y(x)\,\eu^{-ax}. Então

z(x)=y(x)eaxay(x)eax=(y(x)ay(x))eax=0,z'(x) = y'(x)\,\eu^{-ax} - a\,y(x)\,\eu^{-ax} = \bigl(y'(x) - a y(x)\bigr)\eu^{-ax} = 0,

de modo que zz é constante, digamos CC, e y(x)=Ceaxy(x) = C\eu^{ax}. A condição inicial y(x0)=y0y(x_0) = y_0 obriga C=y0eax0C = y_0 \eu^{-a x_0}, de modo único.

Observação 26.3

A exponencial é, assim, caracterizada pela mais simples de todas as equações diferenciais: crescimento proporcional ao tamanho. É por isso que ela aparece em toda parte na natureza.

A equação y' = y prescreve um coeficiente angular em cada ponto do plano (segmentos cinzentos). As soluções C x — aqui C = 1 (azul), C = 0.3 (vermelho) e C = -0.5 (laranja) — são exatamente as curvas que seguem esse campo de inclinações.
A equação y=yy' = y prescreve um coeficiente angular em cada ponto do plano (segmentos cinzentos). As soluções CexC\eu^{x} — aqui C=1C = 1 (azul), C=0.3C = 0.3 (vermelho) e C=0.5C = -0.5 (laranja) — são exatamente as curvas que seguem esse campo de inclinações.

Exemplo 26.4 (Decaimento radioativo)

Uma quantidade radioativa satisfaz N=λNN' = -\lambda N com λ>0\lambda > 0, logo N(t)=N0eλtN(t) = N_0\,\eu^{-\lambda t}; veja o Exemplo 23.9 para a meia-vida.

26.2 A equação y=ay+by' = ay + b

Teorema 26.5 (Soluções de y=ay+by' = ay + b)

Sejam a0a \neq 0 e bRb \in \R. As soluções em R\R de y=ay+by' = ay + b são exatamente as funções

y(x)=Ceaxba,CR.y(x) = C\,\eu^{ax} - \frac{b}{a}, \qquad C \in \R .

A função constante yp=bay_p = -\frac ba é a solução de equilíbrio. Para cada (x0,y0)(x_0, y_0) existe uma única solução com y(x0)=y0y(x_0) = y_0; se a<0a < 0, toda solução tende ao equilíbrio ba-\frac ba quando x+x \to +\infty.

Demonstração. A constante yp=bay_p = -\frac{b}{a} satisfaz yp=0=ayp+by_p' = 0 = a y_p + b. Agora yy é solução se, e somente se,

(yyp)=y=ay+b=a(yyp)+ayp+b=0=a(yyp),(y - y_p)' = y' = ay + b = a(y - y_p) + \underbrace{a y_p + b}_{=\,0} = a (y - y_p),

isto é, se e somente se z=yypz = y - y_p resolve z=azz' = az. Pelo Teorema 26.2, z=Ceaxz = C\eu^{ax}, donde a fórmula, a existência e a unicidade. Se a<0a < 0, eax0\eu^{ax} \to 0 quando x+x \to +\infty, de modo que y(x)bay(x) \to -\frac ba.

Método 26.6 (Resolver y=ay+by' = ay + b com condição inicial)

  1. Encontre a solução de equilíbrio yp=bay_p = -\frac{b}{a} (resolva y=0y' = 0).
  2. Escreva a solução geral y=Ceax+ypy = C\eu^{ax} + y_p.
  3. Determine CC pela condição inicial.
  4. Confronte o comportamento a longo prazo com a intuição física (a solução converge para o equilíbrio?).

A mesma estratégia — solução particular mais solução geral da equação homogênea — estende-se a y=ay+f(x)y' = ay + f(x): veja a Proposição 26.8.

Exemplo 26.7 (Lei do resfriamento de Newton)

Uma xícara de café a 8080\,^\circC está numa sala a 2020\,^\circC. A lei de Newton afirma que a temperatura TT satisfaz T=k(T20)T' = -k(T - 20) para algum k>0k > 0, isto é, T=kT+20kT' = -kT + 20k. O equilíbrio é 2020, e T(t)=20+60ektT(t) = 20 + 60\,\eu^{-kt}: o café esfria exponencialmente depressa até a temperatura ambiente.

Qualquer que seja sua temperatura inicial, todas as soluções de T' = -k(T - 20) convergem exponencialmente para o equilíbrio T = 20.
Qualquer que seja sua temperatura inicial, todas as soluções de T=k(T20)T' = -k(T - 20) convergem exponencialmente para o equilíbrio T=20T = 20.

26.3 A equação y=ay+f(x)y' = ay + f(x)

Proposição 26.8 (Estrutura do conjunto de soluções)

Seja ff contínua em um intervalo II e seja ypy_p uma solução particular de

y=ay+f(x)y' = ay + f(x)

em II. Então as soluções em II são exatamente as funções y=Ceax+ypy = C\eu^{ax} + y_p, CRC \in \R.

Demonstração. Como no Teorema 26.5: yy é solução se, e somente se, z=yypz = y - y_p satisfaz z=(ay+f)(ayp+f)=azz' = (ay + f) - (ay_p + f) = az, se e somente se z=Ceaxz = C\eu^{ax}.

Método 26.9 (Adivinhar uma solução particular)

Procure uma solução particular da mesma forma que ff:

  • ff polinomial de grau nn: tente um polinômio de grau nn;
  • f(x)=αekxf(x) = \alpha\,\eu^{kx} com kak \neq a: tente yp=βekxy_p = \beta\,\eu^{kx};
  • f(x)=αeaxf(x) = \alpha\,\eu^{ax} (caso ressonante): tente yp=βxeaxy_p = \beta x\,\eu^{ax}.

Substitua na equação e identifique os coeficientes.

Exemplo 26.10

Resolva y=2y+4xy' = 2y + 4x. Tente yp=αx+βy_p = \alpha x + \beta: α=2(αx+β)+4x\alpha = 2(\alpha x + \beta) + 4x para todo xx obriga 2α+4=02\alpha + 4 = 0 e α=2β\alpha = 2\beta, logo α=2\alpha = -2 e β=1\beta = -1. Solução geral: y=Ce2x2x1y = C\eu^{2x} - 2x - 1.

26.4 Exercícios

Exercício 26.1

Resolva em R\R: (a) y=3yy' = 3y com y(0)=2y(0) = 2; (b) 2y+y=02y' + y = 0 com y(0)=1y(0) = -1; (c) y=y+5y' = -y + 5 com y(0)=0y(0) = 0.

Solução

Solução de Exercício 26.1.

(a) y=Ce3xy = C\eu^{3x}; y(0)=2y(0) = 2y=2e3xy = 2\eu^{3x}.

(b) y=12yy' = -\frac12 y, logo y=Cex/2y = C\eu^{-x/2}; y(0)=1y(0) = -1y=ex/2y = -\eu^{-x/2}.

(c) Equilíbrio yp=5y_p = 5; y=Cex+5y = C\eu^{-x} + 5; y(0)=0y(0) = 0C=5C = -5: y=5(1ex)y = 5\left(1 - \eu^{-x}\right).

Exercício 26.2

Uma população de bactérias cresce a uma taxa proporcional ao seu tamanho, dobrando a cada 33 horas. Escreva a equação diferencial satisfeita pela população N(t)N(t) e determine a constante de proporcionalidade.

Solução

Solução de Exercício 26.2.

N=kNN' = kN, logo N(t)=N0ektN(t) = N_0 \eu^{kt}. Dobrar em 33 horas significa e3k=2\eu^{3k} = 2, isto é,

k=ln230.231 h1.k = \frac{\ln 2}{3} \approx 0.231\ \text{h}^{-1}.

Exercício 26.3

Verifique que yp(x)=xexy_p(x) = x\,\eu^{x} é solução de y=y+exy' = y + \eu^x e dê todas as soluções em R\R.

Solução

Solução de Exercício 26.3.

yp(x)=ex+xex=yp(x)+exy_p'(x) = \eu^x + x\eu^x = y_p(x) + \eu^x: ypy_p é uma solução particular (este é o caso ressonante do Método 26.9). Pela Proposição 26.8, as soluções são y=Cex+xex=(C+x)exy = C\eu^{x} + x\eu^{x} = (C + x)\,\eu^x, CRC \in \R.

Exercício 26.4 ★★

Resolva y=2y+exy' = -2y + \eu^{x} com y(0)=1y(0) = 1. (Sugestão: procure uma solução particular da forma βex\beta\,\eu^{x}.)

Solução

Solução de Exercício 26.4.

Tente yp=βexy_p = \beta\eu^x: βex=2βex+ex\beta\eu^x = -2\beta\eu^x + \eu^x3β=13\beta = 1, logo yp=13exy_p = \frac13\eu^x. Solução geral y=Ce2x+13exy = C\eu^{-2x} + \frac13\eu^{x}; a condição y(0)=1y(0) = 1C=23C = \frac23:

y(x)=23e2x+13ex.y(x) = \frac{2}{3}\,\eu^{-2x} + \frac{1}{3}\,\eu^{x}.

Exercício 26.5 ★★

O carbono-14 decai com meia-vida de 57305730 anos. Uma amostra arqueológica contém 60%60\% do carbono-14 de um organismo vivo. Estime sua idade.

Solução

Solução de Exercício 26.5.

N(t)=N0eλtN(t) = N_0\,\eu^{-\lambda t} com λ=ln25730\lambda = \frac{\ln 2}{5730} (Exemplo 23.9). Resolvemos eλt=0.6\eu^{-\lambda t} = 0.6:

t=ln(1/0.6)λ=5730ln(5/3)ln25730×0.51080.69314220 anos.t = \frac{\ln(1/0.6)}{\lambda} = 5730\,\frac{\ln(5/3)}{\ln 2} \approx 5730 \times \frac{0.5108}{0.6931} \approx 4220 \text{ anos}.

Exercício 26.6 ★★

Um tanque contém 100100 L de água pura. Uma salmoura com 0.20.2 kg de sal por litro entra a 55 L/min, e a mistura (mantida uniforme) sai à mesma vazão. Seja m(t)m(t) a massa de sal no tanque no instante tt (em minutos).

  1. Justifique que m=1m20m' = 1 - \dfrac{m}{20}.
  2. Resolva e determine o limite de m(t)m(t) quando t+t \to +\infty. Interprete.
Solução

Solução de Exercício 26.6.

1. O sal entra a 0.2×5=10.2 \times 5 = 1 kg/min. A saída leva concentração m100\frac{m}{100} kg/L a 55 L/min, isto é, m20\frac{m}{20} kg/min. Logo m=1m20m' = 1 - \frac{m}{20}.

2. Equilíbrio mp=20m_p = 20; m(t)=20+Cet/20m(t) = 20 + C\eu^{-t/20}, e m(0)=0m(0) = 0C=20C = -20:

m(t)=20(1et/20)t+20 kg.m(t) = 20\left(1 - \eu^{-t/20}\right) \xrightarrow[t\to+\infty]{} 20 \text{ kg}.

A longo prazo, a concentração do tanque iguala a da salmoura que entra: 0.20.2 kg/L ×\times 100100 L =20= 20 kg.

Exercício 26.7 ★★

Uma paraquedista de 8080 kg cai sujeita à gravidade (g=9.8 m/s2g = 9.8\ \text{m/s}^2) e a uma resistência do ar proporcional à velocidade, de modo que sua velocidade satisfaz v=gkmvv' = g - \frac{k}{m}v com k=16k = 16 kg/s.

  1. Resolva a equação com v(0)=0v(0) = 0.
  2. Calcule a velocidade terminal limt+v(t)\lim_{t\to+\infty} v(t) e o tempo necessário para atingir 95%95\% dela.
Solução

Solução de Exercício 26.7.

1. v=gkmvv' = g - \frac{k}{m} v com km=1680=0.2\frac km = \frac{16}{80} = 0.2. Equilíbrio v=mgk=9.80.2=49v_\infty = \frac{mg}{k} = \frac{9.8}{0.2} = 49 m/s; v(t)=49+Ce0.2tv(t) = 49 + C\eu^{-0.2t}, e v(0)=0v(0) = 0

v(t)=49(1e0.2t).v(t) = 49\left(1 - \eu^{-0.2 t}\right).

2. Velocidade terminal 4949 m/s (176\approx 176 km/h). Queremos 1e0.2t=0.951 - \eu^{-0.2t} = 0.95, isto é, e0.2t=0.05\eu^{-0.2t} = 0.05:

t=ln200.23.000.215 s.t = \frac{\ln 20}{0.2} \approx \frac{3.00}{0.2} \approx 15 \text{ s}.

Exercício 26.8 ★★★

(Equação logística.) Uma população y(t)(0,1)y(t) \in \intoo{0}{1} (como fração da população máxima) satisfaz

y=y(1y).y' = y(1 - y).
  1. Seja z=1yz = \dfrac{1}{y}. Mostre que zz satisfaz a equação linear z=z+1z' = -z + 1.
  2. Resolva em zz e depois em yy, com y(0)=110y(0) = \frac{1}{10}.
  3. Mostre que y(t)1y(t) \to 1 quando t+t \to +\infty e esboce a forma da curva solução.
Solução

Solução de Exercício 26.8.

1. z=1yz = \frac1yz=yy2=y(1y)y2=1yy=1y+1=z+1z' = -\frac{y'}{y^2} = -\frac{y(1-y)}{y^2} = -\frac{1-y}{y} = -\frac1y + 1 = -z + 1.

2. Equilíbrio zp=1z_p = 1, logo z(t)=1+Cetz(t) = 1 + C\eu^{-t}. De y(0)=110y(0) = \frac{1}{10} vem z(0)=10z(0) = 10, logo C=9C = 9 e

y(t)=11+9et.y(t) = \frac{1}{1 + 9\,\eu^{-t}} .

3. Quando t+t \to +\infty, 9et09\eu^{-t} \to 0 e y(t)1y(t) \to 1. A curva é a clássica curva logística em S (sigmoide): crescimento lento no início (yy pequeno, yyy' \approx y), crescimento máximo quando y=12y = \frac12 (onde y=y(1y)y' = y(1-y) é máximo) e depois saturação rumo à capacidade de suporte 11.

Exercício 26.9 ★★★

Seja yy uma solução de y=ay+by' = ay + b e seja un=y(n)u_n = y(n) para nNn \in \N. Mostre que (un)(u_n) satisfaz uma recorrência aritmético-geométrica un+1=qun+ru_{n+1} = q u_n + r e exprima qq e rr em função de aa e bb. A que corresponde a condição q<1\abs{q} < 1 para a equação diferencial?

Solução

Solução de Exercício 26.9.

Pelo Teorema 26.5, y(x)=Ceaxbay(x) = C\eu^{ax} - \frac ba. Logo

un+1=Cea(n+1)ba=ea(Ceanba)+ba(ea1)=qun+ru_{n+1} = C\eu^{a(n+1)} - \frac ba = \eu^{a}\left(C\eu^{an} - \frac ba\right) + \frac ba\left(\eu^a - 1\right) = q\,u_n + r

com q=eaq = \eu^{a} e r=ba(ea1)r = \frac{b}{a}\left(\eu^{a} - 1\right). A condição q<1\abs q < 1 significa ea<1\eu^a < 1, isto é, a<0a < 0: exatamente a condição sob a qual as soluções da equação diferencial convergem para o equilíbrio ba-\frac ba — e, de fato, o ponto fixo da recorrência é r1q=ba\frac{r}{1 - q} = -\frac{b}{a}.

26.5 Problema: O relógio dentro das coisas

Problema 26.1

Problema de fim de semana — o carbono-14 data as cavernas, a lei do resfriamento de Newton cronometra um crime, e uma equaçãozinha veste quatro trajes

Uma equação diferencial é uma lei de mudança; resolvê-la transforma a lei em um relógio. A mesma equaçãozinha y=ay+by' = ay + b (Teorema 26.5) faz tique-taque dentro de carvão pré-histórico, de café esfriando, de gotas de chuva caindo e de soros hospitalares — e ler esses relógios é o ofício deste problema: ele data as pinturas de Lascaux, fixa uma hora da morte e audita as próprias hipóteses, como faz um bom cientista.

Parte I — Fluência.

  1. Resolva y=3yy' = 3y, y(0)=2y(0) = 2; e depois y=2y+6y' = -2y + 6, y(0)=0y(0) = 0 (Método 26.6).
  2. Rededuza o truque de unicidade por trás do Teorema 26.2: se y=ayy' = ay, calcule a derivada de y(t)eaty(t)\,\eu^{-at} e conclua que toda solução é CeatC\eu^{at}.
  3. Para y=2y+6y' = -2y + 6: encontre a solução de equilíbrio e descreva o destino de todas as demais soluções quando t+t \to +\infty. (Os pontos fixos das recorrências do Problema 13.1, agora contínuos.)
  4. Para a lei de decaimento y=yτy' = -\frac{y}{\tau}: mostre que a meia-vida é t1/2=τln2t_{1/2} = \tau \ln 2.
  5. Esboce (ou descreva) a família de soluções de y=2y+6y' = -2y + 6: o que fazem as soluções que começam acima de 33? E abaixo de 33? E em 33?

Parte II — Carbono-14. O tecido vivo mantém uma proporção constante de carbono-14 radioativo; na morte a entrada cessa e o estoque decai: N=λNN' = -\lambda N, com meia-vida de 57305\,730 anos.

  1. Calcule λ\lambda (por ano).
  2. Um osso retém 20%20\,\% da proporção do vivo: qual é a idade dele?
  3. Carvão das cavernas pintadas de Lascaux retém cerca de 15%15\,\%: date as pinturas.
  4. O homem do gelo, Ötzi, encontrado em uma geleira alpina, media cerca de 53%53\,\%: date-o (os arqueólogos falam em cerca de 53005\,300 anos — como foi o seu resultado?).
  5. Por que o carbono-14 não pode datar dinossauros? Calcule a fração restante após dez meias-vidas, enuncie a fração após 6565 milhões de anos como potência de 22 e conclua.
  6. Barras de erro: se os 53%53\,\% de Ötzi só são conhecidos com precisão de ±1%\pm 1\,\%, calcule a faixa de idade. Que precisão de laboratório compra que precisão de datação?

Parte III — O resfriamento de Newton e um crime. Um corpo à temperatura TT em um ambiente à temperatura constante TaT_a esfria segundo T=k(TTa)T' = -k\,(T - T_a).

  1. Resolva a equação (substitua z=TTaz = T - T_a):

    T(t)=Ta+(T0Ta)ekt.T(t) = T_a + (T_0 - T_a)\,\eu^{-kt} .
  2. Um café servido a 9090\,^\circC em uma sala a 2020\,^\circC marca 7070\,^\circC após 55 minutos. Determine kk e, em seguida, o tempo de espera até os 5555\,^\circC em que se pode beber.
  3. A questão do leite: para tomar o café mais quente possível daqui a dez minutos, o leite frio deve ser acrescentado agora ou no último momento? Responda com a equação (o que acrescentar leite faz com a diferença TTaT - T_a, e como essa diferença comanda a perda?).
  4. Perícia: um corpo é encontrado à meia-noite a 3030\,^\circC em uma sala a 2020\,^\circC; uma hora depois ele marca 2828\,^\circC. Supondo 3737\,^\circC no momento da morte, recupere kk com as duas medições e calcule a hora da morte.
  5. Audite o relógio do perito: cite três violações reais das hipóteses do modelo e o sentido em que cada uma enviesaria a hora da morte estimada.

Parte IV — Quatro trajes e uma curva em S.

  1. A equação logística do Exercício 26.8, y=y(1y)y' = y(1 - y) com y(0)=0.1y(0) = 0.1: usando a substituição do exercício, deduza y(t)=11+9ety(t) = \dfrac{1}{1 + 9\eu^{-t}} e encontre o instante do ponto de inflexão y=12y = \frac12 — a data do epidemiologista do Problema 22.1, agora calculável.
  2. Uma gota de chuva obedece a v=10v2v' = 10 - \frac v2 (gravidade menos arrasto), com v(0)=0v(0) = 0. Encontre a velocidade terminal, resolva em v(t)v(t) e calcule quando a gota atinge 95%95\,\% da velocidade terminal.
  3. Um soro hospitalar libera um medicamento a taxa constante enquanto o organismo o elimina proporcionalmente: c=40.5cc' = 4 - 0.5c, com c(0)=0c(0) = 0. Encontre a concentração de regime e o tempo para alcançar metade dela. Enuncie em uma frase o paralelo com as recorrências de empréstimo do Exercício 26.9.
  4. Final — uma equação, quatro trajes: decaimento, resfriamento, queda, dosagem, todos y=ay+by' = ay + b com sinais e nomes diferentes. Recite o ciclo de modelagem que este problema percorreu quatro vezes (lei \to equação \to solução \to calibração \to previsão \to auditoria) e nomeie o fenômeno que exige yy'' — e que problema já o encontrou.
Solução

Solução de Problema 26.1.

1. y=2e3ty = 2\eu^{3t}. Para a segunda: equilíbrio 33, logo y=3+Ce2ty = 3 + C\eu^{-2t} com y(0)=0y(0) = 0: C=3C = -3: y=3(1e2t)y = 3\left(1 - \eu^{-2t}\right).

2. (yeat)=yeatayeat=(yay)eat=0\left(y\eu^{-at}\right)' = y'\eu^{-at} - ay\eu^{-at} = (y' - ay)\eu^{-at} = 0: o produto é uma constante CC, logo y=Ceaty = C\eu^{at} — nenhuma solução escapa.

3. Equilíbrio: y3y \equiv 3. Toda outra solução é 3+Ce2t3 + C\eu^{-2t} com C0C \neq 0: a exponencial morre e a solução desliza até 33 — um ponto fixo estável, o gêmeo contínuo dos pontos fixos das recorrências do Problema 13.1.

4. y=y0et/τy = y_0\eu^{-t/\tau} cai à metade quando et/τ=12\eu^{-t/\tau} = \frac12: t=τln2t = \tau\ln 2.

5. Todas as soluções são translações verticais do decaimento rumo a 33: as que começam acima caem até 33, as que começam abaixo sobem até 33, e a solução constante 33 fica parada — um funil em torno do equilíbrio.

6. λ=ln257301.21×104\lambda = \frac{\ln 2}{5730} \approx 1.21 \times 10^{-4} por ano.

7. eλt=0.2\eu^{-\lambda t} = 0.2: t=ln5λ13300t = \frac{\ln 5}{\lambda} \approx 13\,300 anos.

8. t=ln(1/0.15)λ15700t = \frac{\ln(1/0.15)}{\lambda} \approx 15\,700 anos: os touros de Lascaux são do fim da era glacial.

9. t=ln(1/0.53)λ5250t = \frac{\ln(1/0.53)}{\lambda} \approx 5\,250 anos: dentro de uma vida humana do valor dos arqueólogos — o relógio funciona.

10. Dez meias-vidas deixam 2100.1%2^{-10} \approx 0.1\,\%: no limite da medição. Após 6565 milhões de anos a fração é 265000000/57302113442^{-65\,000\,000/5\,730} \approx 2^{-11\,344}: nenhum átomo do estoque original resta em qualquer fóssil — os dinossauros são datados por relógios mais lentos (potássio-argônio e parentes).

11. 52%52\,\%5405\approx 5\,405 anos e 54%54\,\%5095\approx 5\,095: o ±1%\pm 1\,\% da química vira cerca de ±155\pm 155 anos de história — precisão no laboratório é precisão na etiqueta do museu.

12. z=TTaz = T - T_a satisfaz z=kzz' = -kz: z=(T0Ta)ektz = (T_0 - T_a)\eu^{-kt}, e T=Ta+zT = T_a + z.

13. A diferença inicial é T0Ta=70T_0 - T_a = 70; após cinco minutos ela é 7020=5070 - 20 = 50, logo 50=70e5k50 = 70\eu^{-5k}: k=ln(70/50)50.067k = \frac{\ln(70/50)}{5} \approx 0.067 por minuto. Para beber: 5520=35=70ekt55 - 20 = 35 = 70\eu^{-kt}: t=ln2k10.3t = \frac{\ln 2}{k} \approx 10.3 minutos.

14. A taxa de perda é proporcional à diferença TTaT - T_a: café fervendo sangra calor mais depressa. Acrescentar o leite agora corta a diferença imediatamente, de modo que se perde menos calor ao longo dos dez minutos; deixar para o fim permite que o café esfrie a toda velocidade antes. Para a xícara mais quente: leite primeiro. (Anfitriões impacientes têm a termodinâmica ao contrário.)

15. Diferenças em relação ao ambiente: à meia-noite 1010, uma hora depois 88: ek=0.8\eu^{-k} = 0.8: k=ln1080.223k = \ln\frac{10}{8} \approx 0.223 por hora. Na morte a diferença era 1717; ela decaiu até 1010 à meia-noite: s=ln(17/10)k2.4s = \frac{\ln(17/10)}{k} \approx 2.4 horas. Hora da morte: por volta das 21:4021{:}40.

16. Uma sala aquecida ou com corrente de ar (TaT_a não constante) entorta o relógio para qualquer lado; roupas ou a massa corporal mudam kk (as tabelas do perito corrigem isso) — um kk errado reescala todo o intervalo; e um corpo movido de outro lugar reinicia TaT_a no meio do decaimento, forjando uma morte mais cedo ou mais tarde. A equação é honesta; o risco está nas hipóteses.

17. z=1yz = \frac1y obedece a z=1zz' = 1 - z: z=1+9etz = 1 + 9\eu^{-t} (de z(0)=10z(0) = 10), logo y=11+9ety = \frac{1}{1 + 9\eu^{-t}}. Inflexão em y=12y = \frac12: 9et=19\eu^{-t} = 1: t=ln92.2t = \ln 9 \approx 2.2 — a data de virada da epidemia, direto da fórmula.

18. Velocidade terminal: v=0v' = 0 em v=20v = 20 m/s. Solução: v(t)=20(1et/2)v(t) = 20\left(1 - \eu^{-t/2}\right). Então 0.950.95: et/2=0.05\eu^{-t/2} = 0.05: t=2ln206t = 2\ln 20 \approx 6 s — as gotas de chuva alcançam sua velocidade de cruzeiro em segundos, e é por isso que a chuva não mata.

19. Regime permanente: c=8c = 8. Metade disso: 4=8(1et/2)4 = 8\left(1 - \eu^{-t/2}\right): et/2=12\eu^{-t/2} = \frac12: t=2ln21.4t = 2\ln 2 \approx 1.4 hora. O soro é o gêmeo contínuo do empréstimo: entrada constante, saída proporcional — o Exercício 26.9 torna o dicionário exato.

20. Decaimento (a<0a < 0, b=0b = 0), resfriamento (y=TTay = T - T_a), queda (a<0a < 0, b>0b > 0: aproximação do equilíbrio por baixo), dosagem (a mesma coisa, numa veia): uma equação, quatro mundos. O ciclo: enunciar a lei da mudança; escrever a equação; resolver (Método 26.6); calibrar as constantes com medições; prever; e então auditar as hipóteses (questão 16). A oscilação exige y=ω2yy'' = -\omega^2 y — encontrada, resolvida e posta a balançar no Problema 24.1.