Mostrar código-fonte de Equação diferencial

{{mais notas|data=Maio de 2022}}
{{wikificação|data=fevereiro de 2024}}
[[Imagem:Equazione differenziale.png|thumb|Soluções de uma equação diferencial (a negro) e as respectivas condições iniciais (a vermelho)]]
Em [[matemática]], uma '''equação diferencial''' é uma [[equação]] cuja incógnita é uma [[função (matemática)|função]] que aparece na equação sob a forma das respectivas [[derivada]]s. Dada uma variável {{mvar|x}}, função de uma variável {{mvar|y}}, a equação diferencial envolve, {{mvar|x}}, {{mvar|y}}, derivadas de {{mvar|y}} e eventualmente também derivadas de {{mvar|x}}.<ref name=mendel>{{citar livro|url=http://books.google.com.br/books?id=4bsoo853eEgC|publicado= Bookman|isbn = 9788560031092|último = Mendelson|primeiro = Elliot|último2= Ayres Jr|primeiro2= Frank|título= Teoria e problemas de cálculo}}</ref> Por exemplo:

<math display="block">x^2 \frac{d^2 y}{dx^2} + x \frac{dy}{dx} + (x^2)y = 0</math>

Equações diferenciais têm propriedades intrinsecamente interessantes como:
* solução pode existir ou não;
* caso exista, a solução é única ou não.

A ordem da equação diferencial é a ordem da derivada de maior grau que aparece na equação.<ref name=mendel/> A solução de uma equação diferencial de ordem {{mvar|n}}, conterá {{mvar|n}} constantes.<ref name=mendel/>

As equações diferenciais são usadas para construir modelos matemáticos de fenómenos físicos tais como na [[dinâmica de fluidos]] e em [[mecânica celeste]]. Deste modo, o estudo de equações diferenciais é um campo extenso na [[matemática pura]] e na [[matemática aplicada]].

As equações diferenciais têm inúmeras aplicações práticas em [[medicina]], [[engenharia]], [[química]], [[biologia]] e outras diversas áreas do conhecimento. As soluções destas equações são usadas, por exemplo, para projetar [[ponte]]s, [[automóvel|automóveis]], [[avião|aviões]] e [[circuitos elétricos]].

O tópico das equações diferenciais é muito vasto, podendo ser abordado de diferentes maneiras, dependendo do objetivo do estudo. Do ponto de vista da matemática aplicada, tais equações são utilizadas para descrever situações do mundo real, tendo um papel relevante na ligação e interação da matemática com outras ciências: na [[Física]], para descrever as forças sofridas por corpos; na Química, para compreender a velocidade das [[Reação química|reações]]; na Biologia, para modelar o crescimento de populações; na [[Economia]], para formular modelos acerca do mercado; além dos vários ramos da Engenharia.<ref name= Bassanezi-11-12>{{Harvnb|Bassanezi|p=11-12}}</ref>

== Tipos ==
As equações diferenciais dividem-se em dois tipos:
* Uma [[equação diferencial ordinária]] (EDO) contém apenas funções de uma variável e derivadas daquela mesma variável;<ref name=silva>{{citar livro|url= http://books.google.com.br/books?id=BxeQiS_M2d0C|publicado= Editora Universitária UFPE|isbn = 9788573153255|último = Santos|primeiro = José Dias dos|primeiro2= Zanomi Carvalho da|ultimo2= Silva|título= Métodos Numéricos|data= 2006}}</ref>

* Uma [[equação diferencial parcial]] (EDP) contém funções com mais do que uma variável e suas [[derivada parcial|derivadas parciais]].

== Exemplos ==
'''Exemplo 1'''

Equações diferenciais são extremamente importantes para as ciências, pois nos informam como a variação de uma grandeza afeta outras grandezas relacionadas. A lei mais importante de [[Física]] Clássica, a [[segunda lei de Newton]]:
<math display="block"> \vec{F} = m  \vec{a} </math>
é na verdade uma equação diferencial de segunda ordem:
<math display="block"> \vec{F}(\vec{r},t) = m \frac{d^{2} \vec{r} }{dt^2} </math>

Equações diferenciais fazem parte de nosso dia a dia, mesmo que não nos demos conta disto.

No entanto, as equações diferenciais são mais difíceis de resolver do que as equações algébricas comuns. À exceção das equações separáveis, a resolução de cada tipo diferente de equação sem que se conheça a técnica é uma obra homérica. Por isso, cada avanço no campo das equações diferenciais em geral é creditado a um matemático diferente (exceto por [[Leonhard Euler]]).

'''Exemplo 2'''

Mostrando que as funções

<math>
y_1 (x) = e^{5x} \qquad \mathrm{e} \qquad y_2(x) = e^{-3x}
</math>

são soluções da equação diferencial e dos métodos numéricos após as únicas bidimensionais da resolução.

<math>
y''  - 2y' - 15y = 0
</math>

'''Resolução'''

por simples substituição da função e as suas derivadas
vê-se facilmente que cada uma das funções dada é solução:<ref name=Villate>{{citar livro|título=Equações Diferenciais e Equações de Diferenças|local=Porto|nome=Jaime E.|sobrenome=Villate|ano=2011|total-páginas=120|acesso=13 de julho de 2013|url=http://villate.org/publications/Villate_2001_Equacoes_Diferenciais.pdf}}</ref>

<math>
\begin{align}
  &&25 e^{5x} - 10 e^{5x} - 15 e^{5x} = 0 \\
  &&9 e^{-3x} + 6 e^{-3x} - 15 e^{-3x} = 0
\end{align}
</math>

'''Exemplo 3'''

Demonstre que a relação
<math>
    x + y + e^{xy} = 0
</math>

é solução implícita de
<math>
    \bigg(1 + x e^{xy}\bigg) \dfrac{dy}{dx} + 1 + y e^{xy} = 0
</math>

'''Resolução'''

<math>
\begin{align}
  &\dfrac{d}{dx} (x + y + e^{xy}) = 0\\
  &1 + y' + e^{xy} {d(xy) \over dx} = 0\\
  &1 + y' + (y + xy') e^{xy} = 0\\
  &(1 + x e^{xy}  ) {dy \over dx} + 1 + y e^{xy} = 0 \qquad
\end{align}
</math>

== Classificação ==
=== Equações de primeira ordem ===
As equações diferenciais ordinárias de primeira ordem são da forma <math>F(x,y,y') = 0,</math> mas geralmente por meio de simples manipulação algébrica conseguem-se re-escrever na forma de uma ou mais equações:

<math>
  {dy \over dx} = f(x,y)
</math>

A chamada '''forma inversa''' da equação anterior é

<math>
  {dx \over dy} = \frac{1}{f(x,y)}
</math>

Qualquer solução implícita de uma das duas equações é solução da outra, e se a inversa de uma solução explícita <math>y(x)</math> da primeira equação existir, será solução (<math>x(y)</math>) da equação inversa. A equação pode ser também escrita na chamada '''forma diferencial'''.<ref name=Villate />

<math>
  f(x,y) d x - d y = 0
</math>

Existem em geral muitas soluções de uma equação diferencial de primeira ordem. Dado um valor inicial <math>y(x_0) = y_0,</math> é possível calcular a derivada <math>y'</math> no ponto <math>x_0</math> (igual a <math>f(x_0, y_0)</math> segundo a equação diferencial), e geralmente é possível encontrar uma curva (curva [[integral]]) que passe pelo ponto <math>(x_0, y_0)</math> e com derivada igual a <math>f(x,y)</math> em cada ponto. O problema de valores iniciais:

<math>
  {dy \over dx} = f(x,y) \qquad y(x_0) = y_0
</math>

consiste em encontrar a curva integral (ou curvas integrais) que passa pelo ponto <math>(x_0, y_0).</math><ref name=Villate />

=== Existência e unicidade da solução ===
As condições suficientes para a existência de uma solução única de uma equação diferencial de primeira ordem são definidas pelo [[teorema de Picard]]:

==== Teorema de Picard ====
Considere o [[problema de valor inicial]]

<math>
  {dy \over dx} = f(x,y) \qquad g(x_0) = y_0
</math>

se a função <math>f</math> e a derivada parcial de <math>f</math> em função de <math>y</math> são contínuas
numa vizinhança do ponto <math>(x_0, y_0),</math> existe uma solução única <math>y = g(x)</math> em
certa vizinhança do ponto <math>(x_0, y_0)</math> que verifica a condição inicial <math>g(x_0)= y_0.</math>

O intervalo onde existe a solução única pode ser maior ou menor que o intervalo onde a função <math>f</math> e a sua derivada parcial <math>\partial f/\partial y</math> são contínuas (o teorema não permite determinar o tamanho do intervalo).<ref name=Villate />

As condições do teorema de Picard são condições suficientes, mas não necessárias para a existência de solução única. Quando <math>f</math> ou a sua derivada parcial <math>\partial f/\partial y</math> não sejam contínuas, o teorema não nos permite concluir nada: provavelmente existe solução única apesar das duas condições não se verificarem.

'''Exemplo'''

Demonstre que a relação

<math>
    x^2 + y^2 - c^2 = 0
</math>

onde <math>c</math> é uma constante positiva, é solução implícita da equação

<math>
    {dy \over dx} = -\frac{x}{y}
</math>

o que pode se concluir a partir do teorema de Picard?

'''Resolução'''

<math>
\begin{align}
      &2x + 2y y' = 0\\
      &y' = -\frac{x}{y}
\end{align}
</math>

a função <math>f = -x/y</math> e a sua deriavada parcial <math>\partial f/\partial y = x/y^2</math>
são contínuas em quaisquer pontos fora do eixo dos <math>x.</math> A solução implícita dada conduz às soluções únicas:

<math>
      y_1  = \sqrt{c^2 - x^2} \qquad y_2 = -\sqrt{c^2 - x^2}
</math>

no intervalo <math>-c < x < c.</math> O teorema de Picard nada permite concluir nos pontos <math>y = 0,</math> mas segundo o resultado obtido acima vemos que em cada ponto <math>y = 0</math> existem duas soluções, <math>y_1</math> e <math>y_2.</math><ref name=Villate />

== Métodos para resolução de equações diferenciais ==
Existem diferentes métodos para resolução de equações diferenciais e aqui serão apresentados alguns deles.

=== Variáveis separáveis ===
Equações diferenciais com variáveis separáveis são as mais simples de se resolver, utilizamos o seguinte método:

==== Solução por integração ====
Considere a equação diferencial de primeira ordem <math>
    {dy \over dx} = f(x,y).
</math> Quando f não depende de y, isto é <math>
    f(x,y)=g(x),
</math> a equação diferencial <math>
    {dy \over dx} = g(x)
</math> pode ser resolvida por integração. Se g(x) é uma [[função contínua]], então integrando ambos os lados teremos <math>
    y = \int g(x)dx=G(x)+c,
</math> em que G(x) é a antiderivada (integral indefinida) de g(x).<ref name=":0">{{citar livro|titulo=Equações Diferenciais|editora=Cengage Learning}}</ref>

==== Definição ====
O caso apresentado acima é um caso especial de quando f em <math>
    {dy \over dx} = f(x,y)
</math> é um produto de uma função de x por uma função de y. A definição de uma equação diferencial de primeira ordem que é separável ou de variáveis separáveis é <math>
    {dy \over dx} = g(x)h(y),
</math> pois nesse caso podemos manipular a equação com o intuito de deixar as funções que dependem da mesma variável do mesmo lado da equação e a partir daí usar o método da integração de ambos os lados.<ref name=":0" />

==== Exemplo 1: Resolvendo uma equação diferencial separável ====
'''Resolva''' <math>
(1+x)dy-ydx=0
</math>.

'''Solução:''' Dividindo por (1+x)y, podemos escrever <math>
    {dy \over y} = {dx \over (1+x)}
</math>, da qual segue

<math>
    \int{dy \over y} = \int{dx \over (1+x)}
</math>

<math>
    ln|y| = ln|x+1|+c
</math>

<math>
y = e^{ln|x+1|+c} = e^{ln|x+1|}.e^{c}
</math>

<math>
y = |1+x|e^{c} = \pm e^{c}(1+x)
</math>

Podemos escrever <math>
\pm e^{c}
</math> como uma constante D e obtemos <math>
y = D(1+x)
</math>.

'''Solução alternativa:''' Uma vez que cada integral resulta em logaritmo, uma escolha judiciosa da constante de integração <math>\ln|c|</math> em vez de <math>c</math>. Reescrevendo a segunda linha da solução como <math>\ln|y|=\ln|1+x|+\ln|c|</math>, podemos combinar os termos no lado direito usando as propriedades de logaritmo. De <math>\ln|y|=\ln|c(1+x)|</math> obtemos imediatamente <math>y=c(1+x)</math>. Mesmo que as integrais indefinidas não sejam ''todas'' logaritmos, pode ainda ser vantajoso usar <math>\ln|c|</math>. Porém, nenhuma regra geral pode ser dada.<ref name=":0" />

'''Nota''': Não há necessidade de usar duas constantes na integração de uma equação separável, pois se escrevermos <math>
H(y)+c = G(x)+d
</math>, a diferença entre <math>
d-c
</math> poderá ser substituída por uma única constante <math>
a
</math>. Em várias ocasiões renomearemos constantes, múltiplos de constantes ou combinações de constantes por uma única constante.<ref name=":0" />

==== Exemplo 2: Resolvendo um problema de valor inicial ====
'''Resolva:''' <math>(e^{2y}-y)\cos(x){dy \over dx}=e^{y}\sin(2x)</math> , sendo <math>y(0)=0</math>.

'''Solução:''' Dividindo a equação por <math>e^{y}\cos(x)</math> temos:

<math>{e^{2y}-y \over e^{y}}dy={\sin(2x) \over \cos(x)}dx</math>

Antes de integrar, usamos a divisão termo a termo, no lado esquerdo, e a identidade trigonométrica <math>\sin(2x)=2\sin(x)\cos(x)</math> no lado direito. Então da integração por partes segue:

<math>\int (e^{y}-ye^{-y})dy=2\int \sin(x)dx</math>

<math>e^{y}+ye^{-y}+e^{-y}=-2\cos(x)+c</math>

A condição inicial <math>y=0</math> quando <math>x=0</math> implica <math>c=4</math>. Assim, a solução do problema de valor inicial é:

<math>e^{y}+ye^{-y}+e^{-y}=4-2\cos(x)</math>.<ref name=":0" />

=== Equações lineares ===
==== Definição ====
Uma equação diferencial de primeira ordem é chamada equação linear na variável dependente <math>y</math> se for da forma:

<math>a_1(x){dy \over dx}+a_0(x)y=g(x)</math>

Quando <math>g(x)=0</math>, a equação linear é chamada homogênea; do contrário, é não homogênea.

Dividimos ambos os lados da equação acima pelo coeficiente dominante <math>a_1(x)</math>, obtemos uma forma mais conveniente, a forma padrão de uma equação linear, que segue

<math>{dy \over dx}+P(x)y=f(x)</math>

Procuramos uma solução dessa equação em um intervalo <math>I</math> no qual as funções <math>P</math> e <math>f</math> são contínuas.

==== Propriedade ====
Uma propriedade da ED é que a soma de duas soluções: <math>y=y_c+y_p</math>, onde <math>y_c</math> é uma solução da equação homogênea correspondente

<math>{dy \over dx}+P(x)y=0</math>

e <math>y_p</math> é uma solução particular para a equação não homogênea.<ref name=":0" />

==== Procedimento ====
Para resolução de uma equação linear seguimos esses passos.
# Ponha e equação linear dada na forma padrão.
# Identifique <math>P(x)</math> na forma padrão  e então encontre o fator integrante <math>e^{\int P(x)dx}</math> para a ED visando resolver por integração.
# Multiplique a forma padrão da ED pelo fator integrante. O lado esquerdo da equação resultante é automaticamente a derivada do produto do fator integrante por <math>y</math>: <math>{d \over dx}[e^{\int P(x)dx}y]=e^{\int P(x)dx}f(x)</math>.
# Integre ambos os lados dessa última equação.<ref name=":0" />

==== Exemplo 1: Resolvendo uma ED linear homogênea ====
'''Resolva <math>{dy \over dx}-3y=0</math>'''

'''Solução:''' Essa equação pode ser resolvida pelo método das separáveis. Alternativamente, uma vez que a equação já está na forma padrão, vemos que <math>P(x)=-3</math>, e portanto, o fator integrante é <math>e^{ \int (-3) dx}=e^{-3x}</math>. Multiplicamos a equação por esse fator e observamos que:

<math>e^{-3x}{dy \over dx}-3e^{-3x}y=0</math> é o mesmo que <math>{d \over dx}[e^{-3x}y]=0</math>

Integrado ambos os lados da última equação, obtemos <math>e^{-3x}y=c</math>. Resolvendo para <math>y</math>, obtemos a solução explícita <math>y=ce^{3x}</math>, <math>-\infty < x < \infty</math>.<ref name=":0" />

== Equações Diferenciais Reais Ordinárias Homogêneas Lineares a Coeficientes Constantes ==
==== Definição ====
Uma equação diferencial real ordinária homogênea linear a coeficientes constantes de ordem N é qualquer equação da forma.

<math>\sum_{n=0}^N a_n y^{(n)} = 0</math>

Sendo os coeficientes constantes reais e y uma função real dependente da variável real x. Agora, o novo método de solução: se supormos y = b exp(p x), obtemos:<ref>{{citar livro|título=Equações Diferenciais com Aplicações em Modelagem|ultimo=ZILL|primeiro=Dennis|editora=Cengage Learning|ano=2016|local=São Paulo|páginas=141|}}</ref>

<math>\sum_{n=0}^N a_n b  p^{n}exp(px) = 0</math>

<math>\sum_{n=0}^N a_n p^{n} = 0</math>

O que nada mais é que um polinômio de grau N em função de p. As raízes desse polinômio são os valores de p que fazem a afirmação acima ser verdadeira, tal que a solução geral da equação diferencial será uma [[combinação linear]] da solução b exp(p x) utilizando todos os valores de b e p adequados, ou seja:

<math>y= \sum_{n=0}^N b_n e^{p_nx}</math>

Frequentemente, algumas dessas raízes são complexas, tal que p = u + vi, com u e v reais. Pois bem,

<math>e^{x(u_n+v_ni)}=e^{xu_n}e^{xiv_n} = e^{xu_n}(cos(v_nx)+isen(v_nx))=e^{xu_n}cos(v_nx) + ie^{xu_n}sen(v_nx))
</math>

É garantido por teorema que sempre que a solução de uma equação diferencial real resultar em uma função complexa, ambos o módulo da parte real e o módulo da parte complexa serão soluções reais para a equação. Assim sendo, a solução pode ser escrita de uma forma mais geral por:

<math>y= \sum_{n=0}^N b_n [e^{u_nx}cos(v_nx)+e^{u_nx}sen(v_nx)]</math>

=== Equações não lineares ===
Existem muitas equações diferenciais, especialmente não lineares, que não são suscetíveis à solução analítica de algum modo razoavelmente conveniente.<ref>{{citar livro|titulo=Nonlinear Systems|ultimo=P. G.|primeiro=Drazin|editora=|ano=1992|local=Cambridge University|paginas=|acessodata=}}</ref> Métodos numéricos são uma forma de tratar estas equações. Uma outra abordagem possível analisa o caráter geométrico e nos leva a uma compreensão qualitativa do comportamento das soluções, ao invés de uma informação quantitativa de forma detalhada. Considerando um sistema linear homogêneo de segunda ordem com coeficientes constantes, ou seja, um sistema da forma <math>x'=Ax</math>, onde <math>A</math> é uma matriz 2x2 e <math>x</math> é um vetor 2x1, onde supomos que o determinante da matriz de coeficientes não se anula. A análise deste sistema é importante para a determinação dos autovalores da matriz de coeficientes, determinando assim o comportamento de soluções com seus respectivos autovetores, ou seja, os autovalores e autovetores fornecem uma forma de analisar o chamado plano de fase, que mostra o comportamento das soluções, de forma qualitativa (geométrica), em um plano ortogonal, sem necessariamente fornecer uma solução analítica para o problema.

Existe também o caso onde o sistema se comporta como não linear, entretanto próximo aos pontos críticos comporta-se como aproximadamente linear.<ref>{{citar livro|titulo=Computer Aplications to DifferentialEquations|ultimo=J.M.A.|primeiro=Danby|editora=Englewood Cliffs, NJ: Prentice-Hall|ano=1985|local=|paginas=|acessodata=}}</ref> São os chamados ''Sistemas Localmente Lineares.'' Vamos considerar, como um exemplo desta situação, um sistema autônomo bidimensional não linear da forma <math>x'=f(x)</math>. Nosso objetivo principal é investigar o comportamento das trajetórias deste sistema perto de um ponto crítico. Vamos considerar, primeiro, o que significa um sistema não linear estar "próximo" de um sistema linear. Suponhamos então que <math>x'=Ax+g(x)</math>, em que x=0 é um caso isolado do sistema. Isso significa que existe algum círculo em torno da origem para o qual não existem outros pontos críticos. Suponhamos, também, que o determinante da matriz de coeficientes não se anula, de modo que x=0 também é um ponto crítico isolado no sistema linear <math>x'=Ax</math>. Para que o sistema não linear e o sistema linear anteriores estejam próximos, temos que supor que a função <math>g(x)</math> é pequena.<ref>{{citar livro|titulo=Equações Diferenciais Elementares|ultimo=Boyce|primeiro=Willian|editora=LTC|ano=|local=|paginas=|acessodata=}}</ref> Matematicamente, isso equivale a dizer que <math>||g(x)||/||x|| \longrightarrow 0</math> quando <math>x\longrightarrow0</math>, ou seja, o módulo da função é pequeno quando comparado com o módulo de X em torno da origem. Este tipo de sistema é chamado de sistema localmente linear na vizinhança do ponto crítico X=0.

==== Exemplo 3: Resolvendo equações localmente lineares ====
Determine se o sistema <math>\binom{x}{y}' = \begin{pmatrix} 1 & 0 \\ 0 & 0.5 \end{pmatrix}\binom{x}{y} + \binom{-x^2 - xy}{-0.75xy - 0.25y^2}</math> é localmente linear perto de uma vizinhança em torno da origem.

'''Solução:''' Neste problema, é conveniente fazermos uso de coordenadas polares, com <math>x=r\cos(t)</math> e <math>y=rsen(t)</math>. Feito isso, temos que <math>g1(x,y)/r\longrightarrow (-x^2-xy)/r\longrightarrow(-r^2cos^2(t)-r^2sen(t)cos(t))/r = -r(cos^2(t) + sen(t)cos(t))\longrightarrow0</math>, quando <math>r\longrightarrow0</math>. Analogamente, podemos mostrar que a segunda função também tende a 0, sob a mesma análise. Portanto, o sistema dado é localmente linear próximo à origem.

=== Sistemas de Equações Diferenciais Ordinárias Lineares ===
Um sistema de equações diferenciais ordinárias lineares é um sistema, por exemplo, da forma

<math>x'=Ax+By</math>

<math>y'=Cx+Dy</math>

Sendo A, B, C e D constantes reais dadas, e x e y duas funções incógnitas dependentes de uma variável t. Este é o caso mais simples que pode ser chamado de sistema: o 2x2, com duas equações diferenciais dadas e duas funções a serem determinadas. O sistema pode ser resolvido através da seguinte sequência de passos:

1) Isole a função x na segunda equação, e derive os dois lados com respeito a t.

<math>x={y'-Dy \over C}</math>

<math>x'={y''-Dy' \over C}</math>

2) Substitua na primeira equação.

<math>{y''-Dy' \over C}=A{y'-Dy \over C}+By</math>

3) Agora, temos uma equação diferencial ordinária linear de segunda ordem, que trata apenas da função y. Resolva-a.

<math>({1 \over C})y''+({-D-A \over C})y'+({AD \over C}+B)y=0</math>

Digamos que y = exp(rt) para algum r real:

<math>y=e^{(r)t}</math>

<math>y'=(r)e^{(r)t}</math>

<math>y''=(r)^2e^{(r)t}</math>

Substituindo na equação diferencial:

<math>({1 \over C})(r)^2e^{(r)t}+({-D-A \over C})(r)e^{(r)t}+({AD \over C}-B)e^{(r)t}=0</math>

<math>({1 \over C})(r)^2+({-D-A \over C})(r)+({AD \over C}-B)=0</math>

Um polinômio de segundo grau em função de r, que terá duas soluções dadas por

<math>r={-({-D-A \over C})\pm\surd(({-D-A \over C})^2-4({1 \over C})({AD \over C}-B) \over 2({1 \over C})}</math>

A solução geral para y será, portanto, a combinação linear das duas soluções:

<math>y=Uexp({\displaystyle {-({-D-A \over C})+ \surd (({-D-A \over C})^{2}-4({1 \over C})({AD \over C}-B) \over 2({1 \over C})}}t)+Vexp({\displaystyle {-({-D-A \over C})- \surd (({-D-A \over C})^{2}-4({1 \over C})({AD \over C}-B) \over 2({1 \over C})}}t  )</math>

Sendo U e V duas constantes reais arbitrárias. Isso também pode ser escrito como:

<math>P=Uexp(-{({-D-A \over C}) \over 2({1 \over C})})</math>

<math>Q=Vexp(-{({-D-A \over C}) \over 2({1 \over C})})</math>

<math>y=Pexp({\displaystyle { \surd (({-D-A \over C})^{2}-4({1 \over C})({AD \over C}-B) \over 2({1 \over C})}}t)+Qexp({\displaystyle {- \surd (({-D-A \over C})^{2}-4({1 \over C})({AD \over C}-B) \over 2({1 \over C})}}t  )</math>

4) Agora, precisamos determinar a função x.  Substituindo y na primeira equação diferencial do sistema:

<math>x'=Ax+B[{\displaystyle Pexp({\displaystyle {\surd (({-D-A \over C})^{2}-4({1 \over C})({AD \over C}-B) \over 2({1 \over C})}}t)+Qexp({\displaystyle {-\surd (({-D-A \over C})^{2}-4({1 \over C})({AD \over C}-B) \over 2({1 \over C})}}t)} ]</math>

<math>x'-Ax=B[{\displaystyle Pexp({\displaystyle {\surd (({-D-A \over C})^{2}-4({1 \over C})({AD \over C}-B) \over 2({1 \over C})}}t)+Qexp({\displaystyle {-\surd (({-D-A \over C})^{2}-4({1 \over C})({AD \over C}-B) \over 2({1 \over C})}}t)} ]</math>

Multiplicando os dois lados por uma função arbitrária de t:

<math>\mu x'-\mu Ax=B\mu[{\displaystyle Pexp({\displaystyle {\surd (({-D-A \over C})^{2}-4({1 \over C})({AD \over C}-B) \over 2({1 \over C})}}t)+Qexp({\displaystyle {-\surd (({-D-A \over C})^{2}-4({1 \over C})({AD \over C}-B) \over 2({1 \over C})}}t)} ]</math>

Isso pode ser imaginado como uma sentença da forma

<math>\mu x'+\mu'x=(\mu x)'</math>

De forma que

<math>-A \mu x=\mu'x</math>

<math>-A\mu=\mu'={d\mu \over dt}</math>

<math>-Adt={d\mu \over \mu}</math>

Integrando-se os dois lados:

<math>-A\int dt=\int {d\mu \over \mu}</math>

<math>-At+I=ln(\mu)+J</math>

Sendo I e J constantes reais arbitrárias. Se I - J = K, temos

<math>-At+K=ln(\mu)</math>

<math>\mu=exp(-At+K)=Hexp(-At+K)</math>

Se H = exp(K). Substituindo no lado direito da equação diferencial em x, temos

<math>BHexp(-At)[{\displaystyle Pexp({\displaystyle {\surd (({-D-A \over C})^{2}-4({1 \over C})({AD \over C}-B) \over 2({1 \over C})}}t)+Qexp({\displaystyle {-\surd (({-D-A \over C})^{2}-4({1 \over C})({AD \over C}-B) \over 2({1 \over C})}}t)}]=(Hxexp(-At))'</math>

<math>\int BHexp(-At)[{\displaystyle Pexp({\displaystyle {\surd (({-D-A \over C})^{2}-4({1 \over C})({AD \over C}-B) \over 2({1 \over C})}}t)+Qexp({\displaystyle {-\surd (({-D-A \over C})^{2}-4({1 \over C})({AD \over C}-B) \over 2({1 \over C})}}t)}]dt=Hxexp(-At)</math>

<math>BH\int exp(-At){\displaystyle Pexp({\displaystyle {\surd (({-D-A \over C})^{2}-4({1 \over C})({AD \over C}-B) \over 2({1 \over C})}}t)+Qexp({\displaystyle {-\surd (({-D-A \over C})^{2}-4({1 \over C})({AD \over C}-B) \over 2({1 \over C})}}t)}dt=Hxexp(-At)</math>

<math>x={{P\exp(({\displaystyle {\surd (({-D-A \over C})^{2}-4({1 \over C})({AD \over C}-B) \over 2({1 \over C})}}-A)t) \over {\displaystyle {\surd (({-D-A \over C})^{2}-4({1 \over C})({AD \over C}-B) \over 2({1 \over C})}}}+{Q\exp(({\displaystyle {-\surd (({-D-A \over C})^{2}-4({1 \over C})({AD \over C}-B) \over 2({1 \over C})}}-A)t) \over {\displaystyle {-\surd (({-D-A \over C})^{2}-4({1 \over C})({AD \over C}-B) \over 2({1 \over C})}}} \over \exp(-At)}</math>

Assim, as duas funções estão determinadas, e o sistema está resolvido.

== Equações Diferenciais Reais Parciais Lineares Homogêneas a Coeficientes Constantes ==
==== Definição ====
Uma equação diferencial parcial linear homogênea de ordem M que depende de N variáveis independentes é uma equação da forma:

<math>\sum_{n=0}^N \sum_{m=0}^M a_{m,n} {d^{m}y \over dx_n^m} = 0</math>

Sendo a(m,n) coeficientes constantes. O caso particular M = N = 1 e a(0,n) = 0 se reduz a uma expressão do tipo:

<math>a_0 {dy\over dx_0}+a_1 {dy\over dx_1}=0</math>

<math>{dy\over dx_0}=k {dy\over dx_1}</math>

Se k = - a1/a0. Neste ponto, é aplicado o principal e mais simples método para resolução deste tipo de equação diferencial: o Método da Separação de Variáveis, que supõe que:

<math>y=f(x_0)g(x_1)</math>

Logo,

<math>{dy\over dx_0}=f'g</math>

<math>{dy\over dx_1}=fg'</math>

Substituindo na equação, temos:

<math>f'g=kfg'</math>

Sendo x0 e x1 variáveis independentes, então os problemas sobre f e g são independentes, tal que:

<math>f'=kf</math>

<math>g=g'</math>

Ou, alternativamente,

<math>g=kg'</math>

<math>f=f'</math>

O teorema de existência e unicidade nos garante que essas duas formas chegarão na mesma solução, que é a única solução verdadeira a menos de uma constante por equação (o que prova que a hipótese inicial de que y = fg estava correta). Assim, o método transformou uma equação diferencial parcial em duas equações diferenciais ordinárias de fácil resolução, e a solução do problema será o produto das soluções dessas duas equações.

== Software ==
* [[ExpressionsinBar]]:<ref>{{Citar web|titulo=ExpressionsinBar|url=http://www.alelvisoftware.com/Expressions/ExpressionsinBar.html|obra=www.alelvisoftware.com|acessodata=2020-05-17}}</ref> desolve(y'=b*y-a,y)
* [[Maple]]:<ref>{{Citar web|titulo=dsolve - Maple Programming Help|url=https://www.maplesoft.com/support/help/Maple/view.aspx?path=dsolve}}</ref> dsolve
* [[SageMath]]:<ref>{{Citar web|titulo=Basic Algebra and Calculus — Sage Tutorial v9.0|url=http://doc.sagemath.org/html/en/tutorial/tour_algebra.html|obra=doc.sagemath.org|acessodata=2020-05-17}}</ref>  desolve()
* [[Xcas]]:<ref>{{citar web|url=http://www-fourier.ujf-grenoble.fr/~parisse/giac/cascmd_en.pdf|titulo=Symbolic algebra and Mathematics with Xcas|data=|acessodata=|publicado=|ultimo=|primeiro=}}</ref> desolve(y'=k*y,y)

== Lista de equações diferenciais ==
*[[Segunda lei de Newton]]
*[[Equações de Hamilton]]
*[[Decaimento radioativo]]
*[[Equação do pêndulo]]
*[[Equação de Lane-Emden]]
*[[Equação da onda]]
*[[Equações de Maxwell]]
*[[Equação do calor]]
*[[Equação de Laplace]]
*[[Equação de Poisson]]
*[[Equação de Schrödinger]]
*[[Equações de Navier-Stokes]]
*[[Equação de Lotka-Volterra]]
*[[Equações de Cauchy–Riemann|Equações de Cauchy-Riemann]]

== Ver também ==
*[[Equação diferencial linear]]
*[[Equação diferencial de Bernoulli|Equação de Bernoulli]]

{{Referências}}

== Bibliografia ==
* {{citar livro|nome=Rodney Carlos|sobrenome=Bassanezi|título=Equações Diferenciais Ordinárias: Um curso introdutório|editora=Coleção BC&T - UFABC|ano=|local=São Bernardo do Campo, São Paulo|url=http://gradmat.ufabc.edu.br/disciplinas/listas/iedo/notasdeaulas/equacoes-diferenciais-ordinrias-rodney.pdf|ref=harv}}

{{Equações diferenciais}}
{{Análise}}
{{Matemática industrial e aplicada}}
{{Áreas da matemática}}
{{Controle de autoridade}}

[[Categoria:Equações diferenciais| ]]