VPF::[Pascal] Решение системы линейных уравнений методо

Rodoss

Дата 28.1.2007, 22:41 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 2
Регистрация: 8.6.2006

Репутация: нет
Всего: нет

Методом Гаусса, извиняюсь за кривое название темы, это случайно так получилось. Нашел сабж на C++, а может ли кто-нибудь реализовать на Паскале? Вводится количество уравнений и все коэффициенты каждого из них. Система будет выглядеть примерно так:

2х+3y+4z=5
3x+2y+5z=10
х+4y+7z=4
3x+2y+3z=6

В программе должен быть массив:
2 3 4 5
3 2 5 10
1 4 7 4
3 2 3 6

Суть метода Гаусса в том, чтобы свести к треугольному массиву (первое уравнение с одним неизвестным, второе - с двумя и т.д.)

smile

Это сообщение отредактировал(а) Rodoss - 28.1.2007, 22:43

KpoHyc

Дата 28.1.2007, 23:59 (ссылка)

(нет голосов)

Загрузка ...

Шустрый

Профиль
Группа: Участник
Сообщений: 129
Регистрация: 23.12.2006
Где: Санкт-Петербург

Репутация: 1
Всего: 5

Если еще нужно -выложи Сишный код - переожу на паскалевкий...если еще сам не напишу...

--------------------

AScript + Pascal + C -> C++ ->C#Adobe Photoshop 7.0/CS 2.0 + GIMP+ Visual Studio .NET(sp1)/2005 pro(sp1)

Rodoss

Дата 29.1.2007, 19:58 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 2
Регистрация: 8.6.2006

Репутация: нет
Всего: нет

Вообще очень даже нужно)

http://forum.vingrad.ru/topic-37776/15.html#entry285720 - вот в этой теме нашел c-шный код, только не рублю в нем совершенно.

Код

//////////////////////////////////////////////////////////////////////////////
//
//  Solving system of linear equations (Gauss method)
//  (c) Johna Smith, 1996
//
//  Method description:
//    This methos based on excluding variables from equations.
//    Using first equation we exclude x1 from all other equations
//    Using second equation equation we exclude x2 from all equations
//   (excluding first) and so on. So we have triangular matrix from which
//    we can easily get vector X. (AX=B)
//
//////////////////////////////////////////////////////////////////////////////

#include <stdio.h>
#include <math.h>

#define N     4     // size of matrix
#define N1    N+1

float matrix[N][N1]=
                  {{2,5,4,1,28},
                   {1,3,2,1,17},
                   {2,10,9,7,77},
                   {3,8,9,2,54}
                  };

void ShowMatrix(void)
{
  for (int i=0;i<N;i++)
  {
    for (int j=0;j<N;j++)
      printf("%+f*x%d",matrix[i][j],i+1);
    printf("=%f\n",matrix[i][N]);
  }
}

void main(void)
{
  // Variables declaration
  float tmp,x[N1];
  register short int i,j,k;

  // Printing given matrix
  ShowMatrix();
  // Main loop
  for (i=0;i<N;i++)
  {
    // Excluding variable x[i] from equations
    tmp=matrix[i][i];
    for (j=N;j>=i;j--) matrix[i][j]/=tmp;
    for (j=i+1;j<N;j++)
    {
      tmp=matrix[j][i];
      for (k=N;k>=i;k--)
        matrix[j][k]-=tmp*matrix[i][k];
    }
  }

  // Calculating vector x
  x[N-1]=matrix[N-1][N];
  for (i=N-2;i>=0;i--)
  {
    x[i]=matrix[i][N];
    for (j=i+1;j<N;j++) x[i]-=matrix[i][j]*x[j];
  }

  // Printing solution
  printf("\nSolution:\n");
  for (i=0;i<N;i++)
    printf("x%d=%f\n",i+1,x[i]);
}

M
Alexeis

Модератор: правильно используйте подсветку кода http://forum.vingrad.ru/index.php?showtopic=126445

Это сообщение отредактировал(а) Alexeis - 30.1.2007, 15:33

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| Центр помощи \| Следующая тема »