Αρχείο:Instability of Euler's method.svg
Μετάβαση στην πλοήγηση
Πήδηση στην αναζήτηση
Το μέγεθος αυτής της PNG προεπισκόπησης αυτού του SVG το αρχείο: 720 × 540 εικονοστοιχεία. Άλλες αναλύσεις: 320 × 240 εικονοστοιχεία | 640 × 480 εικονοστοιχεία | 1.024 × 768 εικονοστοιχεία | 1.280 × 960 εικονοστοιχεία | 2.560 × 1.920 εικονοστοιχεία.
Πρωτότυπο αρχείο (Αρχείο SVG, ονομαστικό μέγεθος 720 × 540 εικονοστοιχεία, μέγεθος αρχείου: 23 KB)
Σύνοψη
| Περιγραφή |
English: This graph shows that the instability of Euler's method when solving the equation . The black curve shows the exact solution, the blue squares the numerical approximation with step size , and the red circles the numerical solution with . |
| Ημερομηνία | |
| Πηγή | Έργο αυτού που το ανεβάζει |
| Δημιουργός | Jitse Niesen |
Αδειοδότηση
Εγώ, ο κάτοχος των πνευματικών δικαιωμάτων αυτού του έργου, το δημοσιεύω δια του παρόντος υπό την εξής άδεια χρήσης:
 
|
Το αρχείο αυτό έχει διατεθεί με Creative Commons CC0 1.0 Παγκόσμια Εκχώρηση Κοινού Κτήματος. |
| Το πρόσωπο που συσχέτισε ένα έργο με αυτή την πράξη έχει απελευθερώσει αυτό το έργο στην δημόσια σφαίρα παραιτούμενος από όλα τα δικαιώματά του σε αυτό το έργο παγκοσμίως υπό τη νομοθεσία των πνευματικών δικαιωμάτων και όλα τα σχετικά ή παρεμφερή νόμιμα δικαιώματα που είχε στο έργο, στο εύρος που νόμος ορίζει. Έργα υπό την CC0 δεν χρειάζονται απόδοση. Όταν παραθέτετε το έργο, δε χρειάζεται να υπαινιχθείτε έγκριση από το συγγραφέα.
|
Source code
import matplotlib.pyplot as plt
from numpy import linspace, exp
from math import floor
la = -2.3
end_time = 5
exact_times = linspace(0, end_time, 100)
exact_solution = exp(la * exact_times)
euler1_h = 1;
euler1_N = int(floor(end_time / euler1_h))
euler1_times = [ k * euler1_h for k in range(euler1_N + 1) ]
euler1_solution = [ (1 + euler1_h * la) ** k for k in range(euler1_N + 1) ]
euler2_h = 0.7;
euler2_N = int(floor(end_time / euler2_h))
euler2_times = [ k*euler2_h for k in range(euler2_N + 1) ]
euler2_solution = [ (1 + euler2_h * la) ** k for k in range(euler2_N + 1) ]
plt.plot(exact_times, exact_solution, '-k',
euler1_times, euler1_solution, ':sb',
euler2_times, euler2_solution, ':or', markersize = 10)
plt.gca().tick_params(labelsize = 20)
plt.savefig('euler-instability.svg')
Ιστορικό αρχείου
Πατήστε σε μια ημερομηνία/ώρα για να δείτε το αρχείο όπως εμφανιζόταν εκείνη την χρονική στιγμή.
| Ημερομηνία/Ώρα | Μικρογραφία | Διαστάσεις | Χρήστης | Σχόλιο | |
|---|---|---|---|---|---|
| τρέχον | 18:02, 10 Μαρτίου 2012 | | 720 × 540 (23 KB) | wikimediacommons>Jitse Niesen |
Χρήση αρχείου
Η ακόλουθη σελίδα χρησιμοποιεί προς αυτό το αρχείο:
