Можно посмотреть
здесь.
Немного легче объяснение. За второй пункт: здесь всё дело в том, что при очень больших коэффициентах или сильно меняющихся коэффициентах трёхдиагональная матрица становится очень плохо обусловленной, поэтому даже типичное для вычислительных программ округление, скажем, в 16 разряде, влияет на ответ очень и очень плохо. За третье: тоже искал этим летом, было бесполезно. Проще уже самому написать и никому не рассказывать.
