Отправляет email-рассылки с помощью сервиса Sendsay

RFpro.ru: Программирование на языке Pascal

Подписаться Бесплатная «Золотая» новостная рассылка . Подписчиков 2.943 RSS
←   Сентябрь 2011   →
 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30


За последние 60 дней ни разу не выходила


Сайт рассылки: http://prog.rusfaq.ru/pascal
Открыта: 14-03-2002

Автор
Калашников О.А.

Статистика

2.943 подписчиков
0 за неделю
←   Все выпуски   →

RFpro.ru: Программирование на языке Pascal


Хостинг портала RFpro.ru:
Московский хостер
Профессиональный ХОСТИНГ на базе Linux x64 и Windows x64

РАССЫЛКИ ПОРТАЛА RFPRO.RU

Лучшие эксперты по данной тематике

Асмик Гаряка
Статус: Академик
Рейтинг: 8523
∙ повысить рейтинг » 		Орловский Дмитрий
Статус: Советник
Рейтинг: 5925
∙ повысить рейтинг » 		lamed
Статус: Академик
Рейтинг: 5541
∙ повысить рейтинг »

/ КОМПЬЮТЕРЫ И СОФТ / Программирование / Pascal (Паскаль)

Номер выпуска:1207
Дата выхода:29.09.2011, 22:00
Администратор рассылки:Boriss (Академик)
Подписчиков / экспертов:162 / 171
Вопросов / ответов:1 / 1

Консультация # 184040: Здравствуйте, уважаемые эксперты! Прошу вас ответить на следующий вопрос: нужно реализовать решение заданной СЛАУ, учитывая тот факт, что на главной диагонали матрицы могут присутствовать нули. Решить методом Гаусса, матрицу выводить из файла заранее спасибо!...

Консультация # 184040:

Здравствуйте, уважаемые эксперты! Прошу вас ответить на следующий вопрос:
нужно реализовать решение заданной СЛАУ, учитывая тот факт, что на главной диагонали матрицы могут присутствовать нули.
Решить методом Гаусса, матрицу выводить из файла
заранее спасибо!

Дата отправки: 16.09.2011, 21:40
Вопрос задал: Посетитель - 381027 (Посетитель)
Всего ответов: 1
Страница онлайн-консультации »

Консультирует Сергей Бендер (Бакалавр):

Здравствуйте, Посетитель - 381027!

Текст программы:

Код :
const n=5;

var a,aa:array[1..n,1..n] of real;
    b,x:array[1..n] of real;
    i,j,k:integer;
    h:real;
    fa,fb:text;

begin
{Считывание матрицы A}
assign(fa,'a.txt');
reset(fa);
for i:=1 to n do
begin
     for j:=1 to n do
         read(fa,a[i,j]);
     readln(fa);
end;
close(fa);
{Копия матрицы для проверки}
aa:=a;

{Считывание вектора B}
assign(fb,'b.txt');
reset(fb);
for i:=1 to n do
    readln(fb,b[i]);
close(fb);
bb:=b;

{Прямой ход - исключение переменных}
for i:=1 to n-1 do
begin
     j:=i;
     {Перестановка строк при нулевом диагональном элементе}
     while (a[j,i]=0) and (j<=n) do inc(j);
     if j<>i then
     begin
          for k:=1 to n do
          begin
               h:=a[i,k];
               a[i,k]:=a[j,k];
               a[j,k]:=h;
          end;
          h:=b[i];
          b[i]:=b[j];
          b[j]:=h;
     end;
     {Обнуление поддиагональных коэффициентов A}
     for j:=i+1 to n do
     begin
          a[j,i]:=-a[j,i]/a[i,i];
          for k:=i+1 to n do
              a[j,k]:=a[j,k]+a[j,i]*a[i,k];
          b[j]:=b[j]+a[j,i]*b[i]
     end;
end;

{Обратный ход - нахождение корней}
for i:=n 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('Ответ: X=');
for i:=1 to n do
    writeln(x[i]:0:3);

writeln('Проверка: A*x=');
for i:=1 to n do
begin
     b[i]:=0;
     for j:=1 to n do
         b[i]:=b[i]+aa[i,j]*x[j];
     writeln(b[i]:0:3);
end;
readln;

end.


Исходный файл с матрицей A -- a.txt
Код :
1 2 3 4 5
5 4 3 2 1
3 2 4 1 5
5 5 1 1 2
2 2 3 3 4


Исходный файл с вектором b -- b.txt
Код :
10
11
12
13
14


Результат выполнения:
Код :
Ответ: X=
15.200
-14.000
-8.900
6.500
4.700

Проверка: A*x=
10.000
11.000
12.000
13.000
14.000

Консультировал: Сергей Бендер (Бакалавр)
Дата отправки: 27.09.2011, 22:20

5
Спасибо большое!!
-----
Дата оценки: 28.09.2011, 23:20
Рейтинг ответа:

НЕ одобряю 0 одобряю!

Оценить выпуск | Задать вопрос экспертам

главная страница  |  стать участником  |  получить консультацию
техническая поддержка  |  восстановить логин/пароль

Дорогой читатель!
Команда портала RFPRO.RU благодарит Вас за то, что Вы пользуетесь нашими услугами. Вы только что прочли очередной выпуск рассылки. Мы старались. Пожалуйста, оцените его. Если совет помог Вам, если Вам понравился ответ, Вы можете поблагодарить автора - для этого в каждом ответе есть специальные ссылки. Вы можете оставить отзыв о работе портале. Нам очень важно знать Ваше мнение. Вы можете поближе познакомиться с жизнью портала, посетив наш форум, почитав журнал, который издают наши эксперты. Если у Вас есть желание помочь людям, поделиться своими знаниями, Вы можете зарегистрироваться экспертом. Заходите - у нас интересно!
МЫ РАБОТАЕМ ДЛЯ ВАС!

© 2001-2011, Портал RFPRO.RU, Россия
Авторское право: ООО "Мастер-Эксперт Про"
Калашников О.А. | Гладенюк А.Г.
Хостинг: Компания "Московский хостер"
Версия системы: 2011.6.33 от 10.09.2011
В избранное