Итерационные методы решения СЛАУ (Якоби и min невязок)

Нужна программа на Pascal. Ооооочень!

Задание. Решить задания с помощь итерационных методов . Вывести значения приближенных решени, значение номера итерарации, при котором достигнута заданая точность eps, и графики зависимости нормы невязки от номера итерации.

Рассмотрим СЛАУ Ау=f, где А- матрица 10x10 (задается самостоятельно), f -вектор-столбец правых частей, у -вектор-столбец неизвестных. Пусть A=L+D+U? где L-нижнетреугольная, U-верхнетреугольная, D-диагональные части матрицы А. Усл. ||r(k)||<=eps, где r(k)=Ay(k)-f. Можно использовать ||x||2=(sum(xi^2))^(1/2)

1. Метод Якоби(простой итерации)
yi(k+1)=(fi-sum(aij*yj(k)))/aii, i=1,..,n
метод сходится для любого у(о), если А-положительно определенная симметричная матрица с диагональным преобладанием.

2. Метод минимальных невязок
E*(y(k+1)-y(k))/((r(k),Ar(k)/(Ar(k),Ar(k)))+A*y(k)=f