Resolva a seguinte equação diferencial ou problema de valor inicial conforme seja o caso:

(b) 2 x y + x 2 y + y 3 3 d x + x 2 + y 2 d y = 0

Passo 1

Temos aqui uma EDO que está com d x e d y separados e multiplicados por funções de x e y . Ou seja, temos um problema com cara de EDO exata na forma: M x , y d x + N x , y d y = 0

Sendo: M x , y = 2 x y + x 2 y + y 3 3 , N x , y = x 2 + y 2

OK, mas para saber se é mesmo uma EDO exata, precisamos aplicar a condição de Euler: ∂ M ∂ y = ∂ N ∂ x

Vamos aplicá-la então. Derivando M x , y em relação a y , obtemos

∂ M ∂ y = 2 x + x 2 + y 2 .

Já derivando N x , y em relação a x ,

∂ N ∂ x = 2 x .

Opa! As derivadas deram diferentes, então não temos aqui uma EDO exata...

Passo 2

Quem poderá nos salvar agora? Bem, não temos uma EDO exata logo de cara, mas podemos ter uma EDO quase exata que admite um fator integrante. Então, vamos agarrar essa esperança e tentar achá-lo. Existem dois candidatos:

Candidato 1 :

∂ M ∂ y - ∂ N ∂ x N = h x

Candidato 2:

- ∂ M ∂ y - ∂ N ∂ x M = z y

Em que só adotaremos o candidato se h x for função apenas de x ou z y for função apenas de y . OK, já calculamos as derivadas, então é só substituir. No candidato 1 vai ficar

∂ M ∂ y - ∂ N ∂ x N = 2 x + x 2 + y 2 - 2 x x 2 + y 2 = x 2 + y 2 x 2 + y 2 = 1

O resultado deu uma função constante. Pode não parecer intuitivo, mas é perfeitamente aceitável h x = 1 como uma função apenas de x , então a condição para admitir esse candidato é satisfeita e ficaremos com o candidato 1 mesmo.

Passo 3

Encontrado o candidato, o fator integrante I é determinado a partir dele:

I x = e ∫ h x d x

Substituindo o que achamos: I x = e ∫ d x = e x

Assim mesmo, sem a constante por enquanto.

Passo 4

Encontrado o fator integrante, vamos multiplicar a EDO por ele: 2 x y + x 2 y + y 3 3 d x + x 2 + y 2 d y = 0

e x 2 x y + x 2 y + y 3 3 d x + e x x 2 + y 2 d y = 0

E temos novos M e N que vamos chamar de M - e N - :

M - x , y = e x 2 x y + x 2 y + y 3 3 , N - = e x x 2 + y 2

Então, testamos a condição de Euler novamente: ∂ M - ∂ y = ∂ N - ∂ x

∂ M - ∂ y = e x 2 x y + x 2 y + y 3 3 ' = e x 2 x + x 2 + y 2

∂ N - ∂ x = e x x 2 + y 2 ' = e x x 2 + y 2 + 2 x e x = e x 2 x + x 2 + y 2

OPA NOVAMENTE! Só que agora, temos exatamente o que queríamos:

∂ M - ∂ y = ∂ N - ∂ x

Transformamos a nossa equação quase exata numa equação exata. Agora podemos resolvê-la como tal.

Passo 5

A equação exata tem solução implícita na forma:

u x , y = ∫ M - x , y d x = C

Então vamos calcular essa integral:

∫ M - x , y d x = ∫ e x 2 x y + x 2 y + y 3 3 d x

Como estamos integrando em relação a x , tudo que for y vai se comportar como constante. Além disso, vai aparecer uma função de y ao final no lugar da tradicional constante da integral indefinida. Show, então vamos separar essa integral para facilitar: ∫ e x 2 x y + x 2 y + y 3 3 d x = 2 y ∫ x e x d x + y ∫ x 2 e x d x + y 3 3 ∫ e x d x

Vamos resolver integral por integral com bastante calma. Vem comigo!

Para a primeira,

∫ x e x d x ,

vamos usar a integração por partes.

∫ u d v = u v - ∫ v d u

Seguindo a regra do LIATE para escolher u , temos que a função algébrica x ganha da exponencial. Logo, teremos: u = x

d u = d x

O d v é aquilo que sobra no integrando:

d v = e x d x

E o v sai integrando d v sem a constante:

v = e x

Show. Assim vai ficar:

∫ x e x d x = x e x - ∫ e x d x

∫ x e x d x = x e x - e x = e x x - 1

Uma já foi. Vamos agora para outra:

∫ x 2 e x d x

Seguimos a mesma ideia da integral anterior: u = x 2

d u = 2 x d x

d v = e x d x

v = e x

Logo:

∫ x 2 e x d x = x 2 e x - 2 ∫ x e x d x

A integral que apareceu acabou de ser calculada. Então temos: ∫ x 2 e x d x = x 2 e x - 2 e x x - 1

∫ x 2 e x d x = e x x 2 - 2 x + 2

Por último:

∫ e x d x = e x

Assim, vamos ter: ∫ e x 2 x y + x 2 y + y 3 3 d x = 2 y ∫ x e x d x + y ∫ x 2 e x d x + y 3 3 ∫ e x d x

= 2 y e x x - 1 + y e x x 2 - 2 x + 2 + y 3 3 e x + g y

= 2 y x e x - 2 y e x + y x 2 e x - 2 y x e x + 2 y e x + y 3 3 e x + g y

∫ M - x , y d x = y 3 3 e x + g y

Passo 6

Maneiro, mas não podemos escrever a solução final sem encontrar a expressão de g y . Para achá-lo, usamos a relação: ∂ ∂ y ∫ M - x , y d x = N - x , y

Então derivamos o que encontramos anteriormente em relação a y :

∂ ∂ y ∫ M - x , y d x = y 2 e x + g ' y

E comparamos com N - x , y :

y 2 e x + g ' y = e x x 2 + y 2

y 2 e x + g ' y = y 2 e x + x 2 e x

g ' y = x 2 e x

Passo 7

Encontramos g ' y . Para achar g y , fazemos o que? Integramos!

g y = ∫ g ' y d y = ∫ x 2 e x d y = x 2 e x y

Novamente sem a constante. Agora podemos substituir na fórmula geral:

u x , y = ∫ M - x , y d x = C

u x , y = y 3 3 e x + g y = C

u x , y = y 3 3 e x + x 2 e x y = C

Como o fator integrante e x multiplicou ambos lados da equação, a solução da equação exata que encontramos agora é a mesma da equação inicial do problema. 😊

Resposta

u x , y = y 3 3 e x + x 2 e x y = C

Exercícios de Livros Relacionados

Em cada um dos problemas de 1 e 12:(a) Desenhe um campo de direções para a equação diferencial dada.(b) Com base em uma análise do campo de direções, descreva o comportamento das soluções para valores
Em cada um dos problemas de 1 e 12:(a) Desenhe um campo de direções para a equação diferencial dada.(b) Com base em uma análise do campo de direções, descreva o comportamento das soluções para valores
Em cada um dos problemas de 1 e 12:(a) Desenhe um campo de direções para a equação diferencial dada.(b) Com base em uma análise do campo de direções, descreva o comportamento das soluções para valores
Em cada um dos problemas de 1 e 12:(a) Desenhe um campo de direções para a equação diferencial dada.(b) Com base em uma análise do campo de direções, descreva o comportamento das soluções para valores
Em cada um dos problemas de 1 e 12:(a) Desenhe um campo de direções para a equação diferencial dada.(b) Com base em uma análise do campo de direções, descreva o comportamento das soluções para valores