Metody numeryczne w fizyce Laboratorium 1
1. Rozwiąż numerycznie równanie zaniku promieniotwórczego:
U U
dN N
dt ,
implementując algorytm Eulera. Obliczenia przeprowadź dla stałej zaniku τ = 1 s oraz różnych wartości początkowej liczby atomów NU(0) i wielkości kroku czasowego Δt.
Porównaj rozwiązania numeryczne z rozwiązaniem analitycznym sporządź wykresy ilustrujące błąd bezwzględny i względny w funkcji czasu.
2. Rozwiąż to zagadnienie przyjmując stałą zaniku równą τ = 7 × 108 lat -- wprowadź odpowiednie skalowanie.
3. Zapoznaj się z funkcją Matlaba ode23, wykorzystaj ją do rozwiązania tego zagadnienia.
Karol Tarnowski Wrocław, 2019