Решение систем литейных уравнений матричным методом

[mashulya]

Товарищи программисты помогите пожалуйста! Нужно написать прогу которая будет решать систему линейный уравнений с 4мя неизвестными. Программа должна реализоваваться матричным способом. Может кто нибудь уже делал что нибудь подобное? пожалуйста помогите! или может быть у кого нибудь есть хотя бы алгоритм?

[StarWorm]

Приведи пример... А то так сходу сложно во все это вникнуть!!!

[Min]

вобщем тут методом Гаусса проще всего...... гдето у меня уже готовый исходничек валялся.... как найду-скину

[mashulya]

Сообщение от Min

вобщем тут методом Гаусса проще всего...... гдето у меня уже готовый исходничек валялся.... как найду-скину

Это было бы очень здорово! Я буду очень благодарна!!!
Хотя это как я представляю не совсем то что мне надо...

[mashulya]

Сообщение от StarWorm

Приведи пример... А то так сходу сложно во все это вникнуть!!!

Дана система неравенств линейного вида с четырьма неизвестными.
{a1x1+b1x2+c1x3+d1x4<f
a2x1+b2x2+c2x3+d2x4<f
a3x1+b3x2+c3x3+c3x4<f
L=x2+bx4( целевая функция)
}
Вот и надо написать программкукоторая будет это решать матричным способом.

[Min]

ну вот метод Гаусса...... только теперь уже неравенства???
метод Гаусса считает именно систему линейных уравнений, а для неравенств нужен другой способ.....
коррекцию на линейную зависимость и прочую лабуду я не делал:

#include <stdio.h>
#include <conio.h>

int a[100][100];
int N;
float x[100];

void swapFunc(int i1,int i2)
{
     for(int i=0;i<=N;i++)
     {
      int temp=0;
      temp=a[i][i1];
      a[i][i1]=a[i][i2];
      a[i][i2]=temp;
     }
}

int FindNotNull(int w)
{
    for(int i=w;i<N;i++)
     if(a[w][i]!=0)
      return i;
    return -1;
}

void setOther(int MainI,int w)
{
    for(int i=MainI+1;i<=N;i++)
     a[i][w]=a[i][w]*a[MainI][MainI]-a[MainI][w]*a[i][MainI];
    a[MainI][w]=0;
}

void FindAllAndSet(int w)
{
    for(int i=w+1;i<N;i++)
     if(a[w][i]!=0)
      setOther(w,i);
}

void Accum()
{
    x[N-1]=(float)a[N][N-1]/a[N-1][N-1];
    for(int i=N-2;i>=0;i--)
    {
     float summ=0;
     for(int j=i+1;j<N;j++)
      summ+=a[j][i]*x[j];
     x[i]=(a[N][i]-summ)/a[i][i];
    }
}

int main()
{
    scanf("%d",&N);
    for(int i=0;i<N;i++)
     for(int j=0;j<=N;j++)
      scanf("%d",&a[j][i]);
    bool error=false;
    for(int i=0;(i<N)&!error;i++)
    {
     int q=FindNotNull(i);
     if (q==-1)
     {
      error=true;
      continue;
     }
     swapFunc(i,q);
     FindAllAndSet(i);
    }
    if(!error)
    {
     Accum();
     for(int i=0;i<N;i++)
      printf("%f ",x[i]);
    }
    else
    {
    printf("error!!!");
    }
     printf("\n");
     getch();
    return 0;
}

[Min]

ваше последнее условие очень похоже на симплекс-метод, но только если найти нужно конкретно наибольшее или наименьшее значение целевой функции...... в вашем примере значением функции будет являться промежуток значений, т.е. множество.... вы уверены, што найти нужно именно значение целевой функции, а не её наибольшее/наименьшее значение?

[mashulya]

Я уже не знаю... у меня препод не адекватный.... в общем мне нужно что то что решает задачу описанную выше... он сказал что это матричный метод... Я в ужасе...=(((

[mashulya]

Спасибо за код. Наверно попробую загнать ему это...