Solucionador de EDO de Primeira Ordem

Solucionador de EDO de Primeira Ordem

O Solucionador de EDO de Primeira Ordem recebe uma equação diferencial ordinária na forma dy/dx = f(x, y), classifica automaticamente sua estrutura (separável, linear, autônoma, exata ou geral) e produz tanto uma solução simbólica de forma fechada, onde possível, quanto uma solução numérica de alta precisão em todos os casos. Uma visualização do campo de inclinação em tempo real com a curva de solução sobreposta torna o significado geométrico da equação imediatamente óbvio — as soluções são exatamente as curvas tangentes a cada seta.

O Que é uma EDO de Primeira Ordem?

Uma equação diferencial ordinária de primeira ordem envolve uma função desconhecida y(x) e apenas sua primeira derivada y'(x). A forma explícita padrão é:

Combinada com uma condição inicial y(x₀) = y₀, isso define um problema de valor inicial (PVI). O teorema de Picard-Lindelöf garante uma solução única em alguma vizinhança de x₀, desde que f seja Lipschitz contínua em y perto de (x₀, y₀). Geometricamente, o PVI busca a curva única que passa por (x₀, y₀) cuja inclinação em cada ponto corresponde a f naquele ponto — exatamente a curva tangente ao campo de inclinação.

Seis Classes que o Solucionador Reconhece

Método de Forma Fechada: Linear com Coeficientes Constantes

Quando o lado direito se simplifica para dy/dx = a·y + b com constantes a e b, o fator integrante μ(x) = e^(-a·x) fornece uma solução exata. A solução geral é:

Aplicar a condição inicial y(x₀) = y₀ define a constante C e produz a solução particular única. Esta classe única cobre um número enorme de problemas de livros didáticos:

• Crescimento exponencial — dy/dx = k·y, solução particular y(t) = y₀·e^(k·t).
• Decaimento exponencial — dy/dx = -k·y, meia-vida ln 2 / k.
• Lei de resfriamento de Newton — dy/dx = -k·(y - T_amb), a temperatura do corpo relaxa exponencialmente em direção à ambiente.
• Carregamento de um circuito RC — dV/dt = (1/RC)·(V_in - V), a voltagem do capacitor aproxima-se da fonte.
• Depuração de fármacos — farmacocinética de primeira ordem com taxa de eliminação k.

Como Ler um Campo de Inclinação

Em cada ponto da grade (x, y), a ferramenta desenha um pequeno segmento de reta cuja inclinação é igual a f(x, y). Três observações úteis:

• Equilíbrios são pontos onde f(x, y) = 0 — o campo de inclinação é horizontal. Para equações autônomas, estes são pontos fixos y* que satisfazem f(y*) = 0; trajetórias próximas se aproximam (estáveis) ou se afastam (instáveis) de y*.
• Isoclinas são curvas onde f(x, y) é igual a uma constante c, portanto todas as setas ao longo da curva têm a mesma inclinação c.
• As curvas de solução nunca se cruzam (quando f é Lipschitz) — visualmente óbvio, pois duas curvas que se cruzassem precisariam de inclinações diferentes na interseção.

Método Numérico: Runge-Kutta Clássico (RK4)

Dado (x_n, y_n), o próximo valor é computado pela média de quatro estimativas de inclinação:

O RK4 possui erro de truncamento local O(h⁵) e erro global O(h⁴), proporcionando aproximadamente seis dígitos de precisão na contagem de etapas padrão para equações não rígidas. O solucionador integra para fora do ponto inicial em ambas as direções de x e para de forma limpa se a magnitude de y exceder 10¹⁵ — típico para soluções que explodem em tempo finito, como dy/dx = y².

Como Usar Esta Calculadora

1. Insira o lado direito no campo dy/dx = .... Use x e y como variáveis, * para multiplicação, ^ ou ** para potências e funções padrão como sin, cos, exp, log, sqrt. Constantes pi e e são reconhecidas.
2. Especifique a condição inicial (x₀, y₀) — a curva de solução única passará por este ponto.
3. Escolha o intervalo de x sobre o qual plotar o campo de inclinação e a curva de solução. O intervalo de y é ajustado automaticamente a partir da solução integrada.
4. Clique em Resolver & Visualizar. O classificador é executado primeiro; se sua equação corresponder a um padrão de forma fechada (linear com coeficientes constantes), você obterá a resposta simbólica. O campo de inclinação e a curva de solução são sempre renderizados.
5. Alterne o campo de inclinação para ligar ou desligar e focar na curva de solução, ou repita a animação de desenho da curva para ver como a integração progride a partir do ponto inicial.

Exemplo Resolvido: Lei de Resfriamento de Newton

Uma xícara de café a 80 °C esfria em uma sala a 20 °C. A taxa de transferência de calor é proporcional à diferença de temperatura:

Esta é linear com coeficientes constantes (a = -0.1, b = 2). A forma fechada é:

Após 30 minutos: T(30) = 20 + 60·e⁻³ ≈ 22.99 °C. A visualização do campo de inclinação torna o comportamento limite óbvio — cada curva de solução, independentemente da temperatura inicial, aproxima-se assintoticamente da linha horizontal T = 20.

Aplicações Comuns

• Dinâmica populacional — modelos exponenciais, logísticos e de efeito Allee.
• Farmacocinética — absorção e eliminação de fármacos, cálculos de meia-vida.
• Transferência de calor — lei de resfriamento de Newton, modelos de capacitância global.
• Circuitos RC e RL — transientes elétricos lineares de primeira ordem.
• Decaimento radioativo — cadeias de decaimento de isótopo único.
• Tanques de mistura — concentração de um soluto sob fluxo de entrada/saída.
• Objeto em queda com arrasto — análise de velocidade terminal dv/dt = g - kv.

Perguntas Frequentes

O que é uma equação diferencial ordinária de primeira ordem?

Uma equação diferencial ordinária (EDO) de primeira ordem é uma equação da forma dy/dx = f(x, y) que envolve a função desconhecida y(x) e sua primeira derivada. Resolver a EDO significa encontrar a função y(x) cuja derivada corresponde ao lado direito. Com uma condição inicial y(x₀) = y₀, a solução é única sob suposições de regularidade moderadas (teorema de Picard-Lindelöf).

O que é um campo de inclinação?

Um campo de inclinação (ou campo de direção) plota um pequeno segmento de linha em cada ponto da grade (x, y) cuja inclinação é igual a f(x, y). As curvas de solução da EDO são exatamente as curvas que são tangentes a esses segmentos em todos os pontos. O campo de inclinação oferece uma intuição visual instantânea para o comportamento global das soluções sem resolver a equação simbolicamente.

Quais classes de EDOs de primeira ordem esta ferramenta resolve?

A ferramenta classifica automaticamente a equação em uma de: integrável (depende apenas de x, resolvida por integração direta), linear com coeficientes constantes y' = a·y + b (forma fechada completa fornecida), autônoma (depende apenas de y), separável (fatores como g(x)·h(y)), linear com coeficientes variáveis (P(x)·y + Q(x)) ou geral. Para cada classe, uma solução numérica de Runge-Kutta de alta precisão e uma visualização de campo de inclinação são produzidas.

Qual método numérico é utilizado?

O método clássico de Runge-Kutta de quarta ordem (RK4) é aplicado com 300 subetapas em cada direção a partir do ponto inicial. O RK4 possui erro de truncamento local O(h⁵) e é o padrão para EDOs não rígidas nesta escala. O solucionador detecta divergência (overflow ou NaN) e interrompe a integração de forma limpa para que o gráfico permaneça válido.

O que é o método do fator integrante para EDOs lineares?

Para uma EDO linear de primeira ordem y' + P(x)·y = Q(x), multiplique ambos os lados pelo fator integrante μ(x) = e^∫P(x) dx. O lado esquerdo torna-se a derivada exata d/dx[μ·y], então y(x) = (1/μ(x)) · (∫μ(x)·Q(x) dx + C). Quando P e Q são constantes, isso colapsa na forma fechada y = -b/a + C·e^(a·x), que a ferramenta retorna automaticamente.

Esta ferramenta pode lidar com equações rígidas ou sistemas de EDOs?

Este solucionador destina-se a EDOs escalares de primeira ordem não rígidas. Problemas muito rígidos (onde a solução possui múltiplas escalas de tempo diferindo por muitas ordens de magnitude) podem precisar de um método implícito como Euler regressivo ou Rosenbrock; sistemas acoplados requerem um solucionador de valor vetorial. Para esses casos, use um pacote dedicado como o solve_ivp do SciPy ou solucionadores de EDO rígidas especializados.

Leitura Adicional

• Equação diferencial ordinária — Wikipédia
• Slope field — Wikipédia (Inglês)
• Métodos de Runge-Kutta — Wikipédia
• Fator integrante — Wikipédia
• Teorema de Picard-Lindelöf — Wikipédia