Witam
Mam taki układ równań:
Aij=1/3;
Aii=2.2+2*i
Bi=Aij*j ( suma po j )
Jaką metodę tutaj najlepiej użyć??
Witam
Mam taki układ równań:
Aij=1/3;
Aii=2.2+2*i
Bi=Aij*j ( suma po j )
Jaką metodę tutaj najlepiej użyć??
Wrzuć do Matlaba czy innego Octave'a, odwróć i pomnóż przez wektor wyników.
Chodzi wlasnie o to, ze mam to zaimpelementowac w c++ i szukam najlepszej metody bądź tricku.
Pozdrawiam
To wg mnie najlepszą metodą jest wykorzystanie jakiegoś bindinga do LAPACKa. Lepiej niż istniejące biblioteki Fortranowskie tego nie napiszesz.
Jak to jest na zaliczenie to bym się za bardzo nie starał i spróbował użyć zwykłego Gaussa (choć nie wiem, czy jakiś problemów natury numerycznej nie będzie). Przy równaniach liniowych to nawet metoda największego spadku działa (choć wolno).
Gaussa Seidla ( iteracyjna ) czy eliminacja Gaussa ( dokładna )? Bo myślałem nad LU