Решение СЛУ Аx=b с нижней треугольной матрицей
 

Помогите пожалуйста составить программу решения СЛУ Аx=b с нижней треугольной матрицей

такая задачка сложная и за неделю ни одного ответа(

Телепатов здесь нет.

dasha131, я честно пытался понять ее смысл, но не смог :(
Нам бы схемку, аль чертеж, мы б затеяли вертёж (с)
Пример входных/выходных данных, Система линейных уравнений - это еще более менее понятно.
А что делать с нижней треугольной матрицей и как она соотносится с СЛУ, не совсем понятно.

dasha131, решение этой задачи описано практически в любом учебнике линейной алгебры. Уверен, что если забить текст твоего сообщения в гугле, найдётся не меньше 10-ка статей, в которых объясняется, как эту задачу решить. Чем ты хочешь, чтобы тебе помогли?

Смотрика, и впрям надо было спрашивать Google полностью. Разложение_Холецкого

Здесь нужно использовать метод исключения гаусса и lu-разложение

То есть имея нижнюю треугольную матрицу, нужно прямой подстановкой решить систему,
формула прямой подстановки по-моему имеет вот такой вид

n
х(i)= (b(i) - E a(ij)x(j))/a(ii) и i получается равно 1,2,3...n-1,n
j=i+1

E это сумма от j=i+1 до n

Помогите пожалуйста все это в программу связать...

Берёте эту книгу и разбираете --- там и теория есть и отлаженный исходный код на FORTRAN’е 66. Элегантно сделано --- и изложено и закодировано красиво...
Численное решение больших разреженных систем уравнений /
А. Джордж, Дж. Лю М. : Мир, 1984. 333 с.
http://depositfiles.com/files/r6zia7r4r

Нужен другой язык --- переведёте. Хотя для численных методов лучше Fortran’а ничего нет.

Че, незнает никто как делать??

dasha131, а приведенных ссылок недостаточно? Или вы думаете, что все горят желанием расписывать Холецкого на языке программирования? Мы можем поправить существующий код, указать на ошибки и недочеты, но не решать задачу за других. Для этого есть платные услуги в газетах.

Помогите пожалуйста найти ошибку

program Gauss_Classic;
uses crt;

type vec=array[1..10] of real;
mas=array[1..10] of vec;

var a:mas;b,x:vec;
h:real;
i,j,k,n:integer;

Begin {Начало основной программы}
Clrscr;
read(n);
for i:=1 to n do begin
for j:=1 to n do begin
write('a[',i,',',j,']=');
read(a[i,j]);
end;
write('b[',i,']=');
read(b[i]); end;
Clrscr;
writeln('Исходная матрица');
for i:=1 to n do begin
write(i);
for j:=1 to n do
write(a[i,j]:10:5);
writeln(b[i]:10:5) end;
writeln; writeln; writeln;
Обратный ход - нахождение корней}
for i:=n-1 downto 1 do begin
h:=b[i];
for j:=i+1 to n do h:=h-x[j]*a[i,j];
x[i]:=h/a[i,i]; end;
writeln('Корни уравнения');
for i:=1 to n do writeln('x(',i,')=',x[i]:10:5);
readkey
end.