eskutheone/lab-4

Code and report from Lab 4

lab-4's Issues

Komentarz

Na pierwszy rzut oka kod wyglądał dobrze. Wyjątek stanowi użycie niepotrzebnych zmiennych wskaźnikowych, przez co kod traci na czytelności, a rozdzielanie zmiennej i wskaźnika do niej może prowadzić do błędów.

Lab-4/main.cpp

Line 23 in 5406d23

double* t0_ = &t0, * tm_ = &tm, * y0_ = &y0, * l_ = &l; //wskazniki by latwiej wczytywac zmienne bez ciaglego copy-paste

Wszystko jednak rozbija się o formułę:

Lab-4/main.cpp

Lines 110 to 111 in 5406d23

	rzade = round(log(fabs((ye - anali(y0, t, t0)))) / log(fabs((yee[k-1] - anali(y0, t, y0)))));
	rzadrk = round(log(fabs((yr - anali(y0, t, t0)))) / log(fabs((yrr[k-1] - anali(y0, t, y0))))); //szacowanie rzędu zbieżności

co tu nie działa:

1. rząd metody jest funkcją h i nie możemy oszacować go bez uwzględnienia tego jak zmieniało się h między iteracjami.
2. yee[k-1] jest równe ye, gdyż w międzyczasie k zostało zwiększone o 1.
3. funkcja anali jest wywoływana ze złymi argumentami
4. zamiast porównywać dwa rozwiązania z różnym h porównuje się (gdyby pominąć pkt. 2) rozwiązania dla różnego t z tym samym h

jest to o tyle ciekawy przypadek, gdyż błąd 3 maskuje błąd 2, a błąd 2 ukrywa błąd 4.

Co w rezultacie prowadzi to do błędnych wniosków.