Solutore di EDO del Primo Ordine

Solutore di EDO del Primo Ordine

Il Solutore di EDO del Primo Ordine accetta un'equazione differenziale ordinaria nella forma dy/dx = f(x, y), ne classifica automaticamente la struttura (separabile, lineare, autonoma, esatta o generale) e produce sia una soluzione simbolica in forma chiusa ove possibile, sia una soluzione numerica ad alta precisione ovunque. Una visualizzazione live del campo di direzioni con la curva di soluzione sovrapposta rende immediatamente evidente il significato geometrico dell'equazione: le soluzioni sono esattamente le curve tangenti ad ogni freccia.

Cos'è un'EDO del Primo Ordine?

Un'equazione differenziale ordinaria del primo ordine coinvolge solo una funzione incognita y(x) e la sua derivata prima y'(x). La forma esplicita standard è:

Combinata con una condizione iniziale y(x₀) = y₀, definisce un problema ai valori iniziali (IVP). Il teorema di Picard-Lindelöf garantisce una soluzione unica in un intorno di x₀ purché f sia lipschitziana in y vicino a (x₀, y₀). Geometricamente, l'IVP richiede l'unica curva passante per (x₀, y₀) la cui pendenza in ogni punto corrisponde a f in quel punto — esattamente la curva tangente al campo di direzioni.

Sei Classi Riconosciute dal Risolutore

Metodo in Forma Chiusa: Lineare a Coefficienti Costanti

Quando il membro destro si semplifica in dy/dx = a·y + b con costanti a e b, il fattore integrante μ(x) = e^(-a·x) fornisce una soluzione esatta. La soluzione generale è:

L'applicazione della condizione iniziale y(x₀) = y₀ determina la costante C e produce l'unica soluzione particolare. Questa singola classe copre un numero enorme di problemi scolastici:

• Crescita esponenziale — dy/dx = k·y, soluzione particolare y(t) = y₀·e^(k·t).
• Decadimento esponenziale — dy/dx = -k·y, emivita ln 2 / k.
• Legge del raffreddamento di Newton — dy/dx = -k·(y - T_amb), la temperatura del corpo si stabilizza esponenzialmente verso quella ambiente.
• Carica di un circuito RC — dV/dt = (1/RC)·(V_in - V), la tensione del condensatore si avvicina alla sorgente.
• Clearance dei farmaci — farmacocinetica del primo ordine con tasso di eliminazione k.

Leggere un Campo di Direzioni

In ogni punto della griglia (x, y) lo strumento traccia un breve segmento la cui pendenza è pari a f(x, y). Tre osservazioni utili:

• Gli equilibri sono punti in cui f(x, y) = 0 — il campo di direzioni è orizzontale. Per le equazioni autonome, questi sono punti fissi y* che soddisfano f(y*) = 0; le traiettorie vicine si avvicinano (stabili) o si allontanano (instabili) da y*.
• Le isocline sono curve dove f(x, y) è uguale a una costante c, quindi tutte le frecce lungo la curva hanno la stessa pendenza c.
• Le curve di soluzione non si incrociano mai (quando f è Lipschitziana) — visivamente ovvio perché due curve che si incrociano dovrebbero avere pendenze diverse nel punto di intersezione.

Metodo Numerico: Runge-Kutta Classico (RK4)

Dato (x_n, y_n), il valore successivo viene calcolato mediando quattro stime della pendenza:

RK4 ha un errore di troncamento locale O(h⁵) e un errore globale O(h⁴), fornendo circa sei cifre di precisione al numero di passi predefinito per equazioni non-stiff. Il risolutore integra verso l'esterno dal punto iniziale in entrambe le direzioni x e si ferma in modo pulito se l'ampiezza di y supera 10¹⁵ — tipico delle soluzioni che esplodono in tempo finito, come dy/dx = y².

Come Usare Questo Calcolatore

1. Inserisci il membro destro nel campo dy/dx = ... . Usa x e y come variabili, * per la moltiplicazione, ^ o ** per le potenze e le funzioni standard come sin, cos, exp, log, sqrt. Le costanti pi ed e sono riconosciute.
2. Specifica la condizione iniziale (x₀, y₀) — l'unica curva di soluzione passerà per questo punto.
3. Scegli l'intervallo x su cui tracciare il campo di direzioni e la curva di soluzione. L'intervallo y viene adattato automaticamente dalla soluzione integrata.
4. Clicca su Risolvi e Visualizza. Il classificatore viene eseguito per primo; se la tua equazione corrisponde a un modello in forma chiusa (lineare a coefficienti costanti), otterrai la risposta simbolica. Il campo di direzioni e la curva di soluzione vengono sempre generati.
5. Attiva o disattiva il campo di direzioni per concentrarti sulla curva di soluzione, oppure riproduci l'animazione del tracciamento della curva per vedere come progredisce l'integrazione dal punto iniziale.

Esempio Svolto: Legge del Raffreddamento di Newton

Una tazza di caffè a 80 °C si raffredda in una stanza a 20 °C. Il tasso di trasferimento del calore è proporzionale alla differenza di temperatura:

Questa è lineare a coefficienti costanti (a = -0.1, b = 2). La forma chiusa è:

Dopo 30 minuti: T(30) = 20 + 60·e⁻³ ≈ 22.99 °C. La vista del campo di direzioni rende ovvio il comportamento limite — ogni curva di soluzione, indipendentemente dalla temperatura iniziale, tende asintoticamente alla retta orizzontale T = 20.

Applicazioni Comuni

• Dinamica delle popolazioni — modelli esponenziali, logistici e ad effetto Allee.
• Farmacocinetica — assorbimento ed eliminazione dei farmaci, calcoli dell'emivita.
• Trasferimento di calore — legge del raffreddamento di Newton, modelli a capacità termica concentrata.
• Circuiti RC e RL — transitori elettrici lineari del primo ordine.
• Decadimento radioattivo — catene di decadimento di singoli isotopi.
• Vasche di miscelazione — concentrazione di un soluto sotto afflusso / deflusso.
• Oggetto in caduta con attrito — analisi della velocità terminale dv/dt = g - kv.

Domande Frequenti

Cos'è un'equazione differenziale ordinaria del primo ordine?

Un'equazione differenziale ordinaria (EDO) del primo ordine è un'equazione della forma dy/dx = f(x, y) che coinvolge la funzione incognita y(x) e la sua derivata prima. Risolvere l'EDO significa trovare la funzione y(x) la cui derivata corrisponde al membro destro. Con una condizione iniziale y(x₀) = y₀ la soluzione è unica sotto lievi ipotesi di regolarità (teorema di Picard-Lindelöf).

Cos'è un campo di direzioni?

Un campo di direzioni (o slope field) traccia un piccolo segmento di retta in ogni punto della griglia (x, y) la cui pendenza è pari a f(x, y). Le curve di soluzione dell'EDO sono esattamente le curve tangenti a questi segmenti in ogni punto. Il campo di direzioni fornisce un'intuizione visiva istantanea del comportamento globale delle soluzioni senza risolvere l'equazione simbolicamente.

Quali classi di EDO del primo ordine risolve questo strumento?

Lo strumento classifica automaticamente l'equazione in una delle seguenti: integrabile (dipende solo da x, risolta per integrazione diretta), lineare a coefficienti costanti y' = a·y + b (forma chiusa completa fornita), autonoma (dipende solo da y), separabile (si scompone come g(x)·h(y)), lineare a coefficienti variabili (P(x)·y + Q(x)) o generale. Per ogni classe, vengono prodotte una soluzione numerica Runge-Kutta ad alta precisione e una visualizzazione del campo di direzioni.

Quale metodo numerico viene utilizzato?

Viene applicato il metodo classico di Runge-Kutta del quarto ordine (RK4) con 300 sotto-fasi in ogni direzione dal punto iniziale. RK4 ha un errore di troncamento locale O(h⁵) ed è lo standard per EDO non-stiff a questa scala. Il risolutore rileva la divergenza (overflow o NaN) e interrompe l'integrazione in modo pulito affinché il grafico rimanga valido.

Cos'è il metodo del fattore integrante per le EDO lineari?

Per una EDO lineare del primo ordine y' + P(x)·y = Q(x), moltiplica entrambi i membri per il fattore integrante μ(x) = e^∫P(x) dx. Il lato sinistro diventa la derivata esatta d/dx[μ·y], quindi y(x) = (1/μ(x)) · (∫μ(x)·Q(x) dx + C). Quando P e Q sono costanti, questo collassa nella forma chiusa y = -b/a + C·e^(a·x), che lo strumento restituisce automaticamente.

Questo strumento può gestire equazioni stiff o sistemi di EDO?

Questo risolutore è destinato a EDO scalari del primo ordine non-stiff. Problemi molto stiff (dove la soluzione ha scale temporali multiple che differiscono di molti ordini di grandezza) potrebbero richiedere un metodo implicito come Eulero all'indietro o Rosenbrock; i sistemi accoppiati richiedono un risolutore a valori vettoriali. Per questi casi, utilizza un pacchetto dedicato come solve_ivp di SciPy o risolutori specifici per EDO stiff.

Ulteriori Letture

• Equazione differenziale ordinaria — Wikipedia
• Campo di direzioni — Wikipedia
• Metodi di Runge-Kutta — Wikipedia
• Fattore integrante — Wikipedia
• Teorema di esistenza e unicità (Picard-Lindelöf) — Wikipedia