Risolutore di equazioni differenziali

Un'equazione differenziale (ED) è un'equazione che mette in relazione una funzione con le sue derivate. Un'equazione differenziale ordinaria (EDO) coinvolge una funzione di una sola variabile:

$F\left(x, y, y', y'', \ldots, y^{(n)}\right) = 0$

L'ordine di una ED è la derivata più alta che vi compare. Il grado è la potenza della derivata di ordine più alto (quando l'equazione è polinomiale nelle derivate).

EDO del primo ordine: $y' = f(x, y)$

EDO del secondo ordine: $y'' + p(x)y' + q(x)y = g(x)$

Una soluzione è una funzione $y(x)$ che soddisfa l'equazione su un certo intervallo. La soluzione generale contiene costanti arbitrarie (una per ordine). Un problema di Cauchy (PVI) specifica condizioni come $y(x_0) = y_0$ per determinare un'unica soluzione particolare.

Le equazioni differenziali modellano fenomeni reali: crescita di popolazioni, decadimento radioattivo, sistemi massa-molla, circuiti elettrici, conduzione del calore e flusso dei fluidi.

Metodo 1: Separazione delle variabili

Per equazioni della forma $\frac{dy}{dx} = f(x)g(y)$:

1. Separa: $\frac{dy}{g(y)} = f(x)\,dx$
2. Integra entrambi i membri: $\int \frac{dy}{g(y)} = \int f(x)\,dx$

Esempio: $\frac{dy}{dx} = 2xy$ → $\frac{dy}{y} = 2x\,dx$ → $\ln|y| = x^2 + C$ → $y = Ae^{x^2}$

Metodo 2: Fattore integrante (lineare del primo ordine)

Per $y' + P(x)y = Q(x)$, moltiplica per il fattore integrante $\mu(x) = e^{\int P(x)\,dx}$:

$\frac{d}{dx}[\mu(x) \cdot y] = \mu(x) \cdot Q(x)$

Poi integra entrambi i membri per trovare $y$.

Esempio: $y' + 2y = e^{-x}$. Qui $P(x) = 2$, quindi $\mu = e^{2x}$. Moltiplica: $(e^{2x}y)' = e^{x}$. Integra: $e^{2x}y = e^x + C$, quindi $y = e^{-x} + Ce^{-2x}$.

Metodo 3: Equazione caratteristica (coefficienti costanti)

Per $ay'' + by' + cy = 0$, risolvi l'equazione caratteristica $ar^2 + br + c = 0$:

Metodo 4: Coefficienti indeterminati

Per $ay'' + by' + cy = g(x)$ dove $g(x)$ è un polinomio, un esponenziale, un seno, un coseno o una combinazione:

1. Trova la soluzione generale dell'equazione omogenea
2. Ipotizza una forma per la soluzione particolare in base a $g(x)$
3. Sostituisci e risolvi per i coefficienti
4. Soluzione generale = omogenea + particolare

Metodo 5: Variazione delle costanti

Un metodo generale per $y'' + p(x)y' + q(x)y = g(x)$ quando le soluzioni omogenee $y_1, y_2$ sono note:

$y_p = -y_1 \int \frac{y_2 g}{W}\,dx + y_2 \int \frac{y_1 g}{W}\,dx$

dove $W = y_1 y_2' - y_2 y_1'$ è il Wronskiano.

Confronto dei metodi

• Dimenticare la costante di integrazione: nella separazione delle variabili, la costante deve essere inclusa prima di risolvere rispetto a $y$, poiché influenza la forma finale della soluzione.
• Fattore integrante errato: il fattore integrante per $y' + P(x)y = Q(x)$ è $e^{\int P(x)\,dx}$. Assicurati che l'equazione sia in forma standard (il coefficiente di $y'$ deve essere 1) prima di identificare $P(x)$.
• Trascurare il caso della radice doppia: quando l'equazione caratteristica ha una radice doppia $r$, la seconda soluzione è $xe^{rx}$, non di nuovo solo $e^{rx}$.
• Ipotesi errata per la soluzione particolare: se la tua ipotesi per $y_p$ è già una soluzione dell'equazione omogenea, moltiplica per $x$ (o $x^2$ se necessario) per ottenere una forma valida.
• Ignorare le condizioni iniziali: la soluzione generale ha costanti arbitrarie. Applica le condizioni iniziali solo dopo aver trovato la soluzione generale completa.