Math ClubMath Club
v1 · padrão canônico

Lição 91 — Introdução às Equações Diferenciais Ordinárias

EDO: equação relacionando função e suas derivadas. Classificação, solução geral vs. particular, modelagem em ciência e engenharia.

Used in: Ano 3 EM — arco cálculo aplicado · Equiv. Spécialité Maths francesa (Terminale) · Equiv. Math III japonês avançado · Equiv. Leistungskurs DE (Klasse 12)

y=f(x,y)y' = f(x,\, y)

Forma normal de uma EDO de 1ª ordem: a derivada de y é dada como função de x e do próprio y. Resolver a EDO é encontrar y(x) que, ao ser derivada, satisfaz essa relação em todo intervalo de validade. Uma condição inicial y(x0)=y0y(x_0) = y_0 determina a solução particular.

Choose your door

Rigorous notation, full derivation, hypotheses

Definição e classificação

Equação Diferencial Ordinária

"Uma equação diferencial é uma equação que contém uma ou mais funções de uma variável independente e as derivadas dessas funções." — OpenStax Calculus Volume 2, §4.1

Classificação

Solução geral e particular

"A solução geral de y=f(x)y' = f(x) é y=F(x)+Cy = F(x) + C, onde FF é uma antiderivada de ff e CC é uma constante arbitrária. Para determinar um valor único para CC, uma condição inicial é necessária." — OpenStax Calculus Volume 2, §4.1

Existência e unicidade (Picard-Lindelöf)

Consequência prática: verifica-se Picard-Lindelöf para toda EDO antes de afirmar unicidade. A EDO y=yy' = \sqrt{|y|}, y(0)=0y(0) = 0 viola a hipótese (derivada parcial descontinua em y=0y = 0) e tem infinitas soluções.

A EDO fundamental: crescimento/decaimento exponencial

y=ky    y(x)=y0ekxy' = ky \implies y(x) = y_0\, e^{kx}
what this means · A derivada de y é proporcional a y. Solução: exponencial. Aparece em juros contínuos, decaimento radioativo, resfriamento de Newton, farmacocinética, crescimento bacteriano.
xyk > 0k < 0y₀

Família de soluções de y' = ky. Crescimento (k maior que 0, curva azul) e decaimento (k menor que 0, curva laranja). Todas partem de y₀ em x = 0.

Exemplos resolvidos

Exercise list

40 exercises · 10 with worked solution (25%)

Application 23Understanding 4Modeling 11Proof 2
  1. Ex. 91.1Application

    Classifique y+yy=0y''' + y'' - y = 0: determine a ordem, se é linear e se é homogênea.

    Show solution
    Maior derivada presente: yy'''. Ordem: 3. As incógnitas yy, yy', yy''' aparecem linearmente (grau 1, sem produtos). É linear. O lado direito é zero, portanto é homogênea.
  2. Ex. 91.2ApplicationAnswer key

    Classifique y+xy=0y' + xy = 0: ordem, linear/não-linear, homogênea.

    Show solution
    Ordem 1. yy' e yy aparecem linearmente (coeficiente de yy' é 1, coeficiente de yy é xx). É linear. Lado direito zero: homogênea.
  3. Ex. 91.3ApplicationAnswer key

    Classifique y+y2=xy' + y^2 = x: ordem e linear/não-linear.

    Show solution
    Ordem 1. O termo y2y^2 eleva yy ao quadrado — não é linear. Classificação: não-linear.
    Show step-by-step (with the why)
    1. Identifique a maior derivada: yy'. Ordem = 1.
    2. Para testar linearidade, escreva na forma padrão: y=xy2y' = x - y^2. O lado direito contém y2y^2 — função não-linear de yy.
    3. Uma EDO linear de 1ª ordem tem a forma a1(x)y+a0(x)y=g(x)a_1(x)y' + a_0(x)y = g(x) — os coeficientes dependem só de xx, e yy aparece com grau 1. Aqui y2y^2 viola essa exigência.
    4. Conclusão: não-linear. Linearidade é sobre os graus de yy e suas derivadas, não sobre xx.
    Macete: substitua yy por cycy (constante cc). Se a equação fica multiplicada por cc em todos os termos, ela é linear. Aqui (cy)2=c2y2(cy)^2 = c^2 y^2 — fator c2c^2, não cc. Não-linear confirmado.
  4. Ex. 91.4Application

    Classifique y+y=sinxy'' + y = \sin x: ordem, linear/não-linear, homogênea.

    Show solution
    Ordem 2. Linear: sim (coeficientes constantes, yy e yy'' com grau 1). Lado direito sinx0\sin x \neq 0: não-homogênea.
  5. Ex. 91.5Application

    Verifique que y=e2xy = e^{2x} é solução de y=2yy' = 2y.

    Show solution
    Calcule yy': derivando y=e2xy = e^{2x}, obtemos y=2e2xy' = 2e^{2x}. Verifique na EDO: y=2e2x=2e2x=2yy' = 2e^{2x} = 2 \cdot e^{2x} = 2y. Igualdade satisfeita para todo xx. Portanto y=e2xy = e^{2x} é solução de y=2yy' = 2y.
  6. Ex. 91.6Application

    Verifique que y=sinxy = \sin x é solução de y+y=0y'' + y = 0.

    Show solution
    Calcule y=cosxy' = \cos x e y=sinxy'' = -\sin x. Substitua: y+y=sinx+sinx=0y'' + y = -\sin x + \sin x = 0. Igualdade satisfeita para todo xx.
  7. Ex. 91.7Application

    Verifique que y=x2+3y = x^2 + 3 é solução do PVI y=2xy' = 2x, y(0)=3y(0) = 3.

    Show solution
    Derive y=x2+3y = x^2 + 3: y=2xy' = 2x. A EDO pede y=2xy' = 2x. Igualdade confirmada. Verifique CI: y(0)=0+3=3y(0) = 0 + 3 = 3. Correto.
  8. Ex. 91.8Application

    Verifique que y=etcosty = e^{-t}\cos t é solução de y+2y+2y=0y'' + 2y' + 2y = 0.

    Show solution
    Calcule y=etcostetsinty' = -e^{-t}\cos t - e^{-t}\sin t e y=2etsinty'' = 2e^{-t}\sin t. Substitua em y+2y+2yy'' + 2y' + 2y: 2etsint+2(etcostetsint)+2etcost=02e^{-t}\sin t + 2(-e^{-t}\cos t - e^{-t}\sin t) + 2e^{-t}\cos t = 0. Verificado.
  9. Ex. 91.9Application

    Mostre que y(t)=Cekty(t) = Ce^{-kt} é a solução geral de y=kyy' = -ky.

    Show solution
    Seja y(t)=Cekty(t) = Ce^{-kt}. Calcule y=kCekt=kCekt=kyy' = -kCe^{-kt} = -k \cdot Ce^{-kt} = -ky. A EDO é satisfeita para qualquer CRC \in \mathbb{R}. Como a EDO é de 1ª ordem, esperamos 1 constante arbitrária — confirmado.
    Show step-by-step (with the why)
    1. Derive y=Cekty = Ce^{-kt}: y=C(k)ekt=kCekty' = C \cdot (-k)e^{-kt} = -kCe^{-kt}.
    2. Note que kCekt=k(Cekt)=ky-kCe^{-kt} = -k \cdot (Ce^{-kt}) = -ky.
    3. Portanto y=kyy' = -ky — a EDO é satisfeita para qualquer valor de CC.
    4. A família y=Cekty = Ce^{-kt} com CC livre é a solução geral. Uma CI y(0)=y0y(0) = y_0 determina C=y0C = y_0.
    Observação: A EDO y=kyy' = -ky é linear de 1ª ordem, portanto a solução geral tem exatamente 1 constante arbitrária. Confirmar o número de constantes é sempre o primeiro teste de completude da solução.
  10. Ex. 91.10ApplicationAnswer key

    A queda livre é modelada por y=gy'' = -g (aceleração gravitacional constante). Encontre a solução geral por dupla integração.

    Show solution
    Integre uma vez: y=gt+C1y' = -gt + C_1. Integre novamente: y=12gt2+C1t+C2y = -\tfrac{1}{2}gt^2 + C_1 t + C_2. EDO de 2ª ordem: solução geral com 2 constantes. Interpretação: C1=v0C_1 = v_0 (velocidade inicial), C2=h0C_2 = h_0 (posição inicial).
  11. Ex. 91.11Application

    Qual é a solução geral de yy=0y'' - y = 0?

    Select the correct option
    Select an option first
    Show solution
    A EDO yy=0y'' - y = 0 tem a forma y=yy'' = y. Buscamos funções cuja segunda derivada seja igual a si mesmas. As exponenciais exe^x e exe^{-x} satisfazem: (ex)=ex(e^x)'' = e^x e (ex)=ex(e^{-x})'' = e^{-x}. Solução geral: y=C1ex+C2exy = C_1 e^x + C_2 e^{-x}. As outras opções: (B) seria solução de y=xyy' = xy; (C) de y+y=0y'' + y = 0 (sinal trocado); (D) de y=0y'' = 0.
  12. Ex. 91.12Understanding

    O que uma condição inicial determina na solução de uma EDO?

    Select the correct option
    Select an option first
    Show solution
    Uma condição inicial especifica o valor de yy (e possivelmente suas derivadas) num ponto x0x_0. Isso seleciona, dentro da família de soluções gerais (parameterizada por constantes C1,C2,C_1, C_2, \ldots), a única curva que passa pelo ponto (x0,y0)(x_0, y_0). Não determina velocidade de convergência numérica nem periodicidade — esses são propriedades da estrutura da EDO.
  13. Ex. 91.13ApplicationAnswer key

    Resolva y=3x2y' = 3x^2, y(0)=5y(0) = 5.

    Show solution
    Integre: y=3x2dx=x3+Cy = \int 3x^2\,dx = x^3 + C. Aplique CI: y(0)=0+C=5    C=5y(0) = 0 + C = 5 \implies C = 5. Solução: y=x3+5y = x^3 + 5.
  14. Ex. 91.14Application

    Resolva y=sinxy' = \sin x, y(0)=1y(0) = 1.

    Show solution
    Integre: y=sinxdx=cosx+Cy = \int \sin x\,dx = -\cos x + C. CI: y(0)=cos0+C=1+C=1    C=2y(0) = -\cos 0 + C = -1 + C = 1 \implies C = 2. Solução: y=cosx+2y = -\cos x + 2.
  15. Ex. 91.15ApplicationAnswer key

    Resolva y=6xy'' = 6x, y(0)=1y(0) = 1, y(0)=2y'(0) = 2.

    Show solution
    Primeira integração: y=6xdx=3x2+C1y' = \int 6x\,dx = 3x^2 + C_1. CI y(0)=2y'(0) = 2: C1=2C_1 = 2. Segunda integração: y=x3+2x+C2y = x^3 + 2x + C_2. CI y(0)=1y(0) = 1: C2=1C_2 = 1. Solução: y=x3+2x+1y = x^3 + 2x + 1.
  16. Ex. 91.16Application

    Resolva y=kyy' = ky, y(0)=y0y(0) = y_0. Expresse em termos de kk e y0y_0.

    Show solution
    Solução geral de y=kyy' = ky: y=Cekty = Ce^{kt}. CI y(0)=y0y(0) = y_0: C=y0C = y_0. Solução particular: y(t)=y0ekty(t) = y_0 e^{kt}. Para k>0k > 0: crescimento exponencial. Para k<0k < 0: decaimento exponencial.
  17. Ex. 91.17Application

    Resolva y=2yy' = 2y, y(0)=5y(0) = 5. Calcule y(3)y(3).

    Show solution
    PVI y=2yy' = 2y, y(0)=5y(0) = 5: solução y(t)=5e2ty(t) = 5e^{2t}. Em t=3t = 3: y(3)=5e65×403,42017y(3) = 5e^6 \approx 5 \times 403{,}4 \approx 2017.
    Show step-by-step (with the why)
    1. A EDO é y=2yy' = 2y — forma y=kyy' = ky com k=2k = 2.
    2. Solução geral: y=Ce2ty = Ce^{2t}.
    3. Aplique CI y(0)=5y(0) = 5: 5=Ce0=C    C=55 = Ce^0 = C \implies C = 5.
    4. Solução particular: y(t)=5e2ty(t) = 5e^{2t}.
    5. Calcule y(3)=5e6y(3) = 5e^6. Como e6403,4e^6 \approx 403{,}4, temos y(3)2017y(3) \approx 2017.
    Macete: para qualquer y=kyy' = ky, a solução é sempre y0ekty_0 e^{kt}. Memorize essa forma — ela aparece em todo modelo de crescimento/decaimento.
  18. Ex. 91.18Application

    Resolva y=0,1yy' = -0{,}1\,y, y(0)=100y(0) = 100. Calcule y(20)y(20).

    Show solution
    Solução: y(t)=100e0,1ty(t) = 100 e^{-0{,}1 t}. Em t=20t = 20: y(20)=100e2100×0,135=13,5y(20) = 100 e^{-2} \approx 100 \times 0{,}135 = 13{,}5. A quantidade decaiu a 13,5% do valor inicial.
  19. Ex. 91.19ApplicationAnswer key

    Resolva y=exy' = e^x, y(0)=0y(0) = 0.

    Show solution
    Integre: y=exdx=ex+Cy = \int e^x\,dx = e^x + C. CI y(0)=0y(0) = 0: 0=e0+C=1+C    C=10 = e^0 + C = 1 + C \implies C = -1. Solução: y=ex1y = e^x - 1.
  20. Ex. 91.20Application

    Resolva y=0y'' = 0, y(0)=1y(0) = 1, y(0)=3y'(0) = -3.

    Show solution
    y=0y'' = 0 integrado duas vezes: y=C1y' = C_1 e y=C1x+C2y = C_1 x + C_2. CI y(0)=1y(0) = 1: C2=1C_2 = 1. CI y(0)=3y'(0) = -3: C1=3C_1 = -3. Solução: y=3x+1y = -3x + 1.
  21. Ex. 91.21Application

    Resolva y=1/xy' = 1/x, y(1)=0y(1) = 0 (domínio x>0x > 0).

    Show solution
    Integre para x>0x > 0: y=1/xdx=lnx+Cy = \int 1/x\,dx = \ln x + C. CI y(1)=0y(1) = 0: 0=ln1+C=0+C    C=00 = \ln 1 + C = 0 + C \implies C = 0. Solução: y=lnxy = \ln x.
  22. Ex. 91.22ApplicationAnswer key

    Resolva y=yy'' = -y, y(0)=1y(0) = 1, y(0)=0y'(0) = 0.

    Show solution
    Solução geral de y=yy'' = -y: y=C1cosx+C2sinxy = C_1 \cos x + C_2 \sin x. CI y(0)=1y(0) = 1: C1=1C_1 = 1. CI y(0)=0y'(0) = 0: y(0)=C1sin0+C2cos0=C2=0y'(0) = -C_1 \sin 0 + C_2 \cos 0 = C_2 = 0. Solução: y=cosxy = \cos x.
    Show step-by-step (with the why)
    1. Reconheça que y=yy'' = -y é o oscilador harmônico. Solução geral: y=C1cosx+C2sinxy = C_1 \cos x + C_2 \sin x (verifique derivando duas vezes).
    2. Aplique y(0)=1y(0) = 1: C1cos0+C2sin0=C1=1C_1 \cos 0 + C_2 \sin 0 = C_1 = 1.
    3. Calcule y=C1sinx+C2cosxy' = -C_1 \sin x + C_2 \cos x.
    4. Aplique y(0)=0y'(0) = 0: C10+C21=C2=0-C_1 \cdot 0 + C_2 \cdot 1 = C_2 = 0.
    5. Solução particular: y=cosxy = \cos x.
    Observação: y+y=0y'' + y = 0 é a equação do movimento harmônico simples (mola sem atrito). As soluções são sempre combinações de senos e cossenos com mesma frequência.
  23. Ex. 91.23Application

    Resolva y+2y+2y=0y'' + 2y' + 2y = 0, y(0)=0y(0) = 0, y(0)=1y'(0) = 1.

    Show solution
    Solução geral de y+2y+2y=0y'' + 2y' + 2y = 0: y=ex(C1cosx+C2sinx)y = e^{-x}(C_1 \cos x + C_2 \sin x). CI y(0)=0y(0) = 0: C1=0C_1 = 0. CI y(0)=1y'(0) = 1: y=ex(C1cosxC1sinxC2cosx+C2sinx+)y' = e^{-x}(-C_1 \cos x - C_1 \sin x - C_2 \cos x + C_2 \sin x + \ldots); com C1=0C_1 = 0, avalia-se y(0)=C2=1y'(0) = C_2 = 1. Solução: y=exsinxy = e^{-x} \sin x.
  24. Ex. 91.24Application

    Capacitor descarrega: V=V/(RC)V' = -V/(RC). Para V0=12V_0 = 12 V e RC=1RC = 1 s, calcule V(2)V(2).

    Show solution
    EDO V=V/(RC)V' = -V/(RC) com RC=1RC = 1: V=VV' = -V. Solução: V(t)=V0etV(t) = V_0 e^{-t}. Para V0=12V_0 = 12, V(2)=12e21,62V(2) = 12 e^{-2} \approx 1{,}62 V. O capacitor descarregou a 13,5% da tensão original em 2 s.
  25. Ex. 91.25Modeling

    Colônia de bactérias dobra a cada hora. População inicial: 100. Escreva a EDO e calcule N(5)N(5).

    Show solution
    EDO: N=rNN' = rN com r=ln2r = \ln 2 h1^{-1} (dobra a cada hora). Solução: N(t)=100e(ln2)t=1002tN(t) = 100 \cdot e^{(\ln 2) t} = 100 \cdot 2^t. Em t=5t = 5: N(5)=10025=10032=3200N(5) = 100 \cdot 2^5 = 100 \cdot 32 = 3200 bactérias.
    Show step-by-step (with the why)
    1. Identifique o tipo: "dobra a cada hora" significa crescimento exponencial.
    2. Escreva a lei: taxa de crescimento proporcional à população atual: N=rNN' = rN.
    3. Determine rr: se N(1)=2N(0)N(1) = 2 N(0) então N0er=2N0    r=ln2N_0 e^r = 2N_0 \implies r = \ln 2.
    4. Solução com N(0)=100N(0) = 100: N(t)=1002tN(t) = 100 \cdot 2^t.
    5. Calcule N(5)=10032=3200N(5) = 100 \cdot 32 = 3200.
    Atalho mental: "dobra a cada período τ\tau" sempre leva a N(t)=N02t/τN(t) = N_0 \cdot 2^{t/\tau}. Generalização direta do decaimento radioativo com sinal positivo.
  26. Ex. 91.26Modeling

    Investimento de R$ 1.000 a 5% a.a. com juros contínuos. Escreva a EDO e calcule o montante em 10 anos.

    Show solution
    EDO: S=0,05SS' = 0{,}05\,S, S(0)=1000S(0) = 1000. Solução: S(t)=1000e0,05tS(t) = 1000\,e^{0{,}05 t} (com tt em anos). Para t=10t = 10: S(10)=1000e0,51648S(10) = 1000\,e^{0{,}5} \approx 1648 reais.
  27. Ex. 91.27Modeling

    Café a 90 °C, sala 25 °C, k=0,04k = 0{,}04 min1^{-1}. Escreva a EDO, resolva, e determine em quanto tempo a temperatura chega a 50 °C.

    Show solution
    EDO de resfriamento: T=0,04(T25)T' = -0{,}04(T - 25), T(0)=90T(0) = 90. Solução: T(t)=25+65e0,04tT(t) = 25 + 65\,e^{-0{,}04 t}. Para T=50T = 50: t=ln(30/65)/0,0419,3t = -\ln(30/65)/0{,}04 \approx 19{,}3 min.
  28. Ex. 91.28Modeling

    Carbono-14 (τ1/2=5730\tau_{1/2} = 5730 anos). Um fóssil tem 25% do C-14 inicial. Calcule sua idade.

    Show solution
    EDO: N=λNN' = -\lambda N, λ=ln2/5730\lambda = \ln 2/5730. Para N(t)/N0=0,25N(t)/N_0 = 0{,}25: eλt=0,25    t=ln(0,25)/λ=2ln2/λ=2×5730=11460e^{-\lambda t} = 0{,}25 \implies t = -\ln(0{,}25)/\lambda = 2\ln 2/\lambda = 2 \times 5730 = 11460 anos (exato: duas meias-vidas para reduzir a 25%).
  29. Ex. 91.29Modeling

    Medicamento: meia-vida de 6 h, dose 200 mg. Escreva a EDO e calcule o quanto resta após 18 h.

    Show solution
    EDO: m=kmm' = -km com meia-vida de 6 h, portanto k=ln2/6k = \ln 2/6. Dose inicial 200 mg: m(t)=200e(ln2/6)t=200(1/2)t/6m(t) = 200\,e^{-(\ln 2/6)t} = 200 \cdot (1/2)^{t/6}. Em t=18t = 18 h: m(18)=200(1/2)3=25m(18) = 200 \cdot (1/2)^3 = 25 mg. Restam 12,5% da dose.
    Show step-by-step (with the why)
    1. Meia-vida de 6 h: λ=ln2/60,1155\lambda = \ln 2/6 \approx 0{,}1155 h1^{-1}.
    2. EDO de decaimento: m=λmm' = -\lambda m.
    3. Solução com m(0)=200m(0) = 200: m(t)=2002t/6m(t) = 200 \cdot 2^{-t/6}.
    4. Em t=18t = 18 h (três meias-vidas): m=200/23=25m = 200/2^3 = 25 mg.
    Macete: três meias-vidas sempre reduzem a quantidade a 1/81/8 do valor original. Conte as meias-vidas antes de calcular — economiza álgebra.
  30. Ex. 91.30ModelingAnswer key

    Aplicação financeira rende Selic de 14,75% a.a. com capitalização contínua. Em quantos anos o capital dobra?

    Show solution
    Selic 14,75% a.a., juros contínuos: S=0,1475SS' = 0{,}1475\,S. Tempo para dobrar: 2=e0,1475t    t=ln2/0,14754,702 = e^{0{,}1475\,t} \implies t = \ln 2/0{,}1475 \approx 4{,}70 anos. Pelo "regra de 70": 70/14,754,770/14{,}75 \approx 4{,}7 — confirmado.
  31. Ex. 91.31Modeling

    Epidemia simplificada: I=rI(1I/N)I' = rI(1 - I/N) (equação logística). Identifique os equilíbrios e descreva o comportamento da solução.

    Show solution
    EDO logística: N=rN(1N/K)N' = rN(1 - N/K). Equilíbrios: N=0N = 0 (instável) e N=KN = K (estável). Para NKN \ll K, o fator (1N/K)1(1 - N/K) \approx 1 e o crescimento é aproximadamente exponencial. Para NKN \to K, o crescimento desacelera até zero. O modelo captura a saturação pela capacidade de suporte.
  32. Ex. 91.32Modeling

    Queda com resistência do ar: mv˙=mgkvm\dot{v} = mg - kv (k>0k > 0). Calcule a velocidade terminal vv_\infty (quando a aceleração cessa).

    Show solution
    Queda livre com resistência: mv˙=mgkvm\dot v = mg - kv. Velocidade terminal: quando v˙=0\dot v = 0, mg=kv    v=mg/kmg = kv_\infty \implies v_\infty = mg/k. Fisicamente: o arrasto aumenta com vv até equilibrar o peso.
  33. Ex. 91.33Modeling

    Iodo-131 (τ1/2=8\tau_{1/2} = 8 dias). Escreva a EDO e calcule quanto resta de 100 g após 24 dias.

    Show solution
    Iodo-131: τ1/2=8\tau_{1/2} = 8 dias, λ=ln2/8\lambda = \ln 2/8. Com N(0)=100N(0) = 100 g: N(t)=100e(ln2/8)tN(t) = 100\,e^{-(\ln 2/8)t}. Após 24 dias (três meias-vidas): N(24)=100/8=12,5N(24) = 100/8 = 12{,}5 g. Clinicamente relevante: determina quando o paciente pode ser liberado do isolamento radioativo.
  34. Ex. 91.34Modeling

    Um ativo desvaloriza 3% ao ano de forma contínua. Escreva a EDO e expresse o valor após 5 anos em termos do valor inicial P0P_0.

    Show solution
    EDO: P=0,03PP' = -0{,}03\,P (deflação de 3% a.a.). Solução: P(t)=P0e0,03tP(t) = P_0\,e^{-0{,}03 t}. Em t=5t = 5 anos: P(5)=P0e0,150,861P0P(5) = P_0\,e^{-0{,}15} \approx 0{,}861\,P_0. O preço cai a 86,1% do valor inicial — perda de 13,9% em poder de compra.
  35. Ex. 91.35ModelingAnswer key

    (Forense) Corpo encontrado às 22 h com temperatura 32 °C. Sala a 21 °C, k=0,374k = 0{,}374 h1^{-1}. Temperatura normal do corpo: 37 °C. Escreva a EDO e encontre kk usando as condições dadas.

    Show solution
    EDO: T=k(T21)T' = -k(T - 21). Dados: T(0)=37T(0) = 37, T(1)=32T(1) = 32. Solução geral: T(t)=21+16ektT(t) = 21 + 16\,e^{-kt}. Usando T(1)=32T(1) = 32: ek=11/16    k=ln(16/11)0,374e^{-k} = 11/16 \implies k = \ln(16/11) \approx 0{,}374 h1^{-1}. Temperatura de vida: 37 °C (t=0t = 0), então tmortet_\text{morte} é incógnita — para fins forenses, resolve-se T(t)=TmedidaT(t) = T_\text{medida} para tt.
  36. Ex. 91.36UnderstandingAnswer key

    O que o Teorema de Picard-Lindelöf garante sobre a EDO y=f(x,y)y' = f(x, y), y(x0)=y0y(x_0) = y_0?

    Select the correct option
    Select an option first
    Show solution
    O Teorema de Picard-Lindelöf exige que ff e f/y\partial f/\partial y sejam contínuas perto de (x0,y0)(x_0, y_0). A continuidade de f/y\partial f/\partial y é equivalente a Lipschitz local em yy. Apenas continuidade de ff garante existência, mas não unicidade (contra-exemplo: y=yy' = \sqrt{\lvert y \rvert}). A opção (C) é incorreta: soluções podem existir atravessando descontinuidades de ff em xx.
  37. Ex. 91.37Understanding

    Por que a solução geral de uma EDO de ordem nn tem exatamente nn constantes arbitrárias? Como isso se relaciona com o número de condições iniciais necessárias?

    Show solution
    Uma EDO de ordem nn requer nn integrações para chegar à solução; cada integração introduz uma constante arbitrária. Portanto a solução geral tem exatamente nn constantes. Para determinar todas elas, precisamos de nn condições: y(x0)y(x_0), y(x0)y'(x_0), …, y(n1)(x0)y^{(n-1)}(x_0).
  38. Ex. 91.38Understanding

    Explique por que y=yy' = \sqrt{\lvert y \rvert}, y(0)=0y(0) = 0 tem infinitas soluções. Qual hipótese do Picard-Lindelöf é violada?

    Show solution
    Para y=yy' = \sqrt{\lvert y \rvert}, y(0)=0y(0) = 0: a função f(y)=yf(y) = \sqrt{\lvert y \rvert} é contínua em y=0y = 0, mas f/y=1/(2y)\partial f/\partial y = 1/(2\sqrt{\lvert y \rvert}) é descontínua (vai a ++\infty) em y=0y = 0. A condição de Lipschitz é violada. Soluções múltiplas: y0y \equiv 0 é solução, e também y=(xc)2/4y = (x - c)^2/4 para x>cx > c e zero antes de cc, para qualquer c0c \geq 0.
  39. Ex. 91.39Proof

    Demonstre que y=Cekxy = Ce^{kx} é a única família de soluções de y=kyy' = ky (a menos de escolha de CC). Dica: considere z(x)=y(x)ekxz(x) = y(x)\,e^{-kx}.

    Show solution
    Seja y(x)=Cekxy(x) = Ce^{kx}. Derive: y=kCekx=kCekx=kyy' = kCe^{kx} = k \cdot Ce^{kx} = ky. Portanto toda função da forma y=Cekxy = Ce^{kx} satisfaz y=kyy' = ky, com CC determinado pela CI. Para mostrar que não há outras soluções: se yy satisfaz y=kyy' = ky, considere z=yekxz = ye^{-kx}. Então z=yekxkyekx=(yky)ekx=0z' = y'e^{-kx} - kye^{-kx} = (y' - ky)e^{-kx} = 0. Logo zz é constante: z=Cz = C, portanto y=Cekxy = Ce^{kx}.
    Show step-by-step (with the why)
    1. Verificação: derive y=Cekxy = Ce^{kx} e substitua em y=kyy' = ky. Confirme que a igualdade vale para todo CC.
    2. Unicidade: seja yy qualquer solução de y=kyy' = ky. Defina z(x)=y(x)ekxz(x) = y(x) e^{-kx}.
    3. Calcule z=yekx+y(k)ekx=(yky)ekxz' = y'e^{-kx} + y \cdot (-k)e^{-kx} = (y' - ky)e^{-kx}.
    4. Como y=kyy' = ky, temos yky=0y' - ky = 0, logo z=0z' = 0 para todo xx.
    5. Uma função com derivada identicamente zero é constante: z(x)=Cz(x) = C.
    6. Portanto y(x)=Cekxy(x) = Ce^{kx} — forma única a menos de CC.
    Curiosidade: a mesma técnica de multiplicar por ekxe^{-kx} (fator integrante) generaliza para EDOs lineares não-homogêneas de 1ª ordem na Lição 93.
  40. Ex. 91.40Proof

    Demonstre que a solução de y=kyy' = ky, y(0)=y0y(0) = y_0 é y(t)=y0ekty(t) = y_0 e^{kt}, usando a técnica de separação de variáveis (preview da Lição 92).

    Show solution
    Separação de variáveis (preview §1.3): escreva dy/y=kdxdy/y = k\,dx. Integre ambos os lados: lny=kx+C0\ln \lvert y \rvert = kx + C_0. Exponencie: y=eC0ekx\lvert y \rvert = e^{C_0} e^{kx}. Defina C=±eC0C = \pm e^{C_0} (ou C=0C = 0 se y0y \equiv 0): y=Cekxy = Ce^{kx}. Aplique CI y(0)=y0y(0) = y_0: C=y0C = y_0. Resultado: y(t)=y0ekty(t) = y_0 e^{kt}.

Fontes

  • Notes on Diffy Qs — Jiří Lebl · 2024 · v6.6 · EN · CC-BY-SA. §0.2–1.3: definição de EDO, classificação, modelagem, exemplos de decaimento radioativo e resfriamento. Fonte primária desta lição.
  • Calculus Volume 2 — OpenStax · 2016 · EN · CC-BY-NC-SA. §4.1–4.3: verificação de soluções, condições iniciais, modelos de crescimento e decaimento, equações separáveis.
  • Active Calculus — Matt Boelkins et al. · 2024 · EN · CC-BY-NC-SA. §7.1–7.2: introdução visual a EDOs, campos de direção, modelagem qualitativa.

Updated on 2026-05-06 · Author(s): Clube da Matemática

Found an error? Open an issue on GitHub or submit a PR — open source forever.