Przy czym dobrano taki układ odniesienia, w którym suma dwóch oddziałujących ze sobą mas wynosi 1 oraz stała grawitacji jest również równa 1. close all clear all clc % rozwiązanie równania d^r/dt^2 = - r/r^3 metodą leapfrog dt = 0.001; % krok czasowy % warunki początkowe r = [0; 0.5; 0.75]; v = [0.25; 0; 1]; r2 = r'*r; a = -r/(r2*sqrt(r2)); % czas symulacji T = 10; res = NaN(round(T/dt)+1,6); idx = 0; for t=0:dt:T idx = idx + 1; v = v + 0.5*a*dt; r = r + v*dt; r2 = r'*r; a = -r/(r2*sqrt(r2)); v = v + 0.5*a*dt; res(idx,:) = [r' v']; end comet3(gca,res(:,1),res(:,2),res(:,3))
Zobacz też metoda Eulera algorytm Rungego-Kutty algorytm Verleta
Przypisy
Linki zewnętrzne , Drexel University
Kategoria:Algorytmy rozwiązywania równań różniczkowych Kategoria:Równania różniczkowe zwyczajne Kategoria:Fizyka matematyczna