Αρχείο:Instability of Euler's method.svg

Τα περιεχόμενα της σελίδας δεν υποστηρίζονται σε άλλες γλώσσες.
Αυτό το αρχείο προέρχεται από το Wikimedia Commons
Από τη Βικιπαίδεια, την ελεύθερη εγκυκλοπαίδεια

Εικόνα σε υψηλότερη ανάλυση(Αρχείο 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 CC-Zero Το αρχείο αυτό έχει διατεθεί με 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')

Λεζάντες

Δεν ορίστηκε λεζάντα

Items portrayed in this file

απεικονίζει

Ιστορικό αρχείου

Κλικάρετε σε μια ημερομηνία/ώρα για να δείτε το αρχείο όπως εμφανιζόταν εκείνη τη στιγμή.

Ώρα/Ημερομ.ΜικρογραφίαΔιαστάσειςΧρήστηςΣχόλια
τελευταία17:02, 10 Μαρτίου 2012Μικρογραφία για την έκδοση της 17:02, 10 Μαρτίου 2012720 × 540 (23 KB)Jitse Niesen

Τα παρακάτω λήμματα συνδέουν σε αυτό το αρχείο:

Καθολική χρήση αρχείου

Τα ακόλουθα άλλα wiki χρησιμοποιούν αυτό το αρχείο: