-
Path: news-archive.icm.edu.pl!news.icm.edu.pl!plix.pl!newsfeed1.plix.pl!goblin2!gobli
n.stu.neva.ru!feeder1.cambriumusenet.nl!feed.tweaknews.nl!209.197.12.242.MISMAT
CH!nx01.iad01.newshosting.com!newshosting.com!newsfeed.neostrada.pl!unt-exc-01.
news.neostrada.pl!unt-spo-b-01.news.neostrada.pl!news.neostrada.pl.POSTED!not-f
or-mail
From: "slawek" <s...@h...pl>
Newsgroups: pl.comp.programming
Subject: Adams-Bashforth - wiki?
Date: Thu, 29 Jul 2010 17:35:02 +0200
MIME-Version: 1.0
Content-Type: text/plain; format=flowed; charset="iso-8859-2"; reply-type=original
Content-Transfer-Encoding: 8bit
X-Priority: 3
X-MSMail-Priority: Normal
Importance: Normal
X-Newsreader: Microsoft Windows Live Mail 14.0.8089.726
X-MimeOLE: Produced By Microsoft MimeOLE V14.0.8089.726
Lines: 62
Message-ID: <4c519fa4$0$19166$65785112@news.neostrada.pl>
Organization: Telekomunikacja Polska
NNTP-Posting-Host: 62.69.202.124
X-Trace: 1280417700 unt-rea-a-01.news.neostrada.pl 19166 62.69.202.124:63921
X-Complaints-To: a...@n...neostrada.pl
Xref: news-archive.icm.edu.pl pl.comp.programming:186191
[ ukryj nagłówki ]Czy przypadkiem w Wikipedii nie ma błędu odnośnie metody Adamsa-Bashforth'a?
http://en.wikipedia.org/wiki/Linear_multistep_method
http://wapedia.mobi/en/Adams_Bashforth
Wzorki po frazie "The Adams-Bashforth methods with s = 1, 2, 3, 4, 5 are"
(ten dla y_{n+5} szczególnie).
Tak samo jest w skrypcie Pawła Góry,
http://student.eldoras.com/UJ/metody_numeryczne/Gora
/wyklad07.pdf, niemniej
jednak nie pytałem go, czy to nie przypadkiem z tego samego źródła (Wiki?) -
może w nim jest błąd? Trzeba będzie pokopać dalej, i/lub wyliczyć owe b[n]
samemu.
Działało u kogoś z s = 5 dobrze?
Mój program ślicznie mi liczy co trzeba... pod warunkiem, że nie daję tego
tasiemca dla n+5 - czyżby ktoś pomylił się przepisując współczynniki?!
Jeżeli jest s = 1, 2, 3, 4 --- to jest dobrze.
Dam s = 5 --- idzie w krzaki.
Ogólnie metoda bardzo przydatna, bo mam przezatentego duże koszty liczenia
f, to takie całkowanie jest brrr...
W A-B mogę użyć f-ów archiwalnych, czysta oszczędność. Tylko dziwi, czemu
danie "lepszego" ciut wzoru prowadzi aż tak w krzaki - powinno ciut
prowadzić wentualnie.
Czy całkuję spline'ami czy trapezami - wisi to rozwiązaniu. Czyli to nie
problem z ciągłością czegoś itd. itp. - nic finezyjnego, wygląda na
literówkę we współczynnikach. Ale żeby w Wikipedii? Wątpię.
Ma ktoś kreatywne pomysły?
slawek
P.S. fragmencik copy/paste - jak wyrzucić pierwsze trzy linie, to chodzi jak
trzeba. A jak jest jak jest, to robi sruuu i wywala rozwiązanie w
nieskończoność.
if (j .gt. 4) then
y1(j+1) = y1(j) + h *
(1901.*f1(j)-2774.*f1(j-1)+2616.*f1(j-2)-1274.*f1(j-
3)*251.*f1(j-4))/720.
y2(j+1) = y2(j) + h *
(1901.*f2(j)-2774.*f2(j-1)+2616.*f2(j-2)-1274.*f2(j-
3)*251.*f2(j-4))/720.
elseif (j .gt. 3) then
y1(j+1) = y1(j) + h *
(55.*f1(j)-59.*f1(j-1)+37.*f1(j-2)-9.*f1(j-3))/24.
y2(j+1) = y2(j) + h *
(55.*f2(j)-59.*f2(j-1)+37.*f2(j-2)-9.*f2(j-3))/24.
elseif (j .gt. 2) then
y1(j+1) = y1(j) + h * (23.*f1(j)-16.*f1(j-1)+5.*f1(j-2))/12.
y2(j+1) = y2(j) + h * (23.*f2(j)-16.*f2(j-1)+5.*f2(j-2))/12.
elseif (j .gt. 1) then
y1(j+1) = y1(j) + h * (3.*f1(j)-f1(j-1))/2.
y2(j+1) = y2(j) + h * (3.*f2(j)-f2(j-1))/2.
else
y1(j+1) = y1(j) + h * f1(j)
y2(j+1) = y2(j) + h * f2(j)
endif
Następne wpisy z tego wątku
- 29.07.10 15:41 slawek
Najnowsze wątki z tej grupy
- We Wrocławiu ruszyła Odra 5, pierwszy w Polsce komputer kwantowy z nadprzewodzącymi kubitami
- Ada-Europe - AEiC 2025 early registration deadline imminent
- John Carmack twierdzi, że gdyby gry były optymalizowane, to wystarczyły by stare kompy
- Ada-Europe Int.Conf. Reliable Software Technologies, AEiC 2025
- Linuks od wer. 6.15 przestanie wspierać procesory 486 i będzie wymagać min. Pentium
- ,,Polski przemysł jest w stanie agonalnym" - podkreślił dobitnie, wskazując na brak zamówień.
- Rewolucja w debugowaniu!!! SI analizuje zrzuty pamięci systemu M$ Windows!!!
- Brednie w wiki - hasło Dehomag
- Perfidne ataki krakerów z KRLD na skrypciarzy JS i Pajton
- Instytut IDEAS może zacząć działać: "Ma to być unikalny w europejskiej skali ośrodek badań nad sztuczną inteligencją."
- Instytut IDEAS może zacząć działać: "Ma to być unikalny w europejskiej skali ośrodek badań nad sztuczną inteligencją."
- Instytut IDEAS może zacząć działać: "Ma to być unikalny w europejskiej skali ośrodek badań nad sztuczną inteligencją."
- U nas propagują modę na SI, a w Chinach naukowcy SI po kolei umierają w wieku 40-50lat
- C++. Podróż Po Języku - komentarz
- "Wuj dobra rada" z KDAB rozważa: Choosing the Right Programming Language for Your Embedded Linux Device
Najnowsze wątki
- 2025-06-30 Kraków => Koordynator Produkcji / Przedstawiciel ds. rozwoju produktu
- 2025-06-30 Środa Wielkopolska => Konsultant wewnętrzny SAP FI/CO <=
- 2025-06-30 Białystok => Programista Mainframe (z/OS, Assembler) <=
- 2025-06-30 Warszawa => International Freight Forwarder <=
- 2025-06-30 Bieruń => Spedytor Międzynarodowy (handel ładunkami/prowadzenie flo
- 2025-06-30 Warszawa => Spedytor Międzynarodowy <=
- 2025-06-30 Lublin => Delphi Programmer <=
- 2025-06-30 Lublin => Programista Delphi <=
- 2025-06-30 Wrocław => Controlling systems Consultant <=
- 2025-06-30 Nowa tarcza do telefonu
- 2025-06-29 Spotkania z Ariane De Rotschild, szefową Iluminatów, Księżniczką Hiszpanii Leonor
- 2025-06-29 Re: Dr. Kontek (już od paru lat nie SGH) odkrył odchylenia statystyczne [PO EKSPERCIE?]
- 2025-06-28 Upadłość i zwolnienia [w Diorze, która była pol prod. głośników - przyp. JMJ]
- 2025-06-28 Taśma izolacyjna do prac elektrycznych
- 2025-06-27 Recenzja 3.1A ;) w 6 gniazdach...