Esempio: soluzione di un sistema di equazioni ODE del primo ordine
Utilizzare un blocco di soluzione e la funzione odesolve per risolvere un sistema di equazioni differenziali ordinarie del primo ordine.
1. Definire il punto finale dell'intervallo della soluzione:
2. Definire il problema utilizzando le derivate e un insieme di condizioni iniziali:
3. Tracciare il grafico delle soluzioni su un singolo intervallo.
Utilizzo della funzione di Rkadapt
Utilizzare la funzione Rkadapt per risolvere lo stesso sistema di equazioni differenziali:
1. Definire una funzione che determina un vettore di valori della derivata per qualsiasi punto di soluzione (t,Y):
2. Definire argomenti supplementari per il solutore ODE:
a. Valore iniziale della variabile indipendente:
b. Vettore dei valori iniziali delle funzioni:
c. Numero di valori di soluzione su [t0, t1]:
3. Utilizzare la funzione Rkadapt per trovare la matrice di soluzione:
| È inoltre possibile utilizzare le funzioni rkfixed, Bulstoer o Radau. |
4. Estrarre i valori della variabile indipendente:
5. Estrarre i valori della prima funzione di soluzione:
6. Estrarre i valori della seconda funzione di soluzione:
7. Estrarre i valori della terza funzione di soluzione:
8. Tracciare i grafici delle tre soluzioni.
| Il grafico delle soluzioni di odesolve è quasi identico a quello delle soluzioni di Rkadapt. |