metodos numericos

[eli_zabet]

ayuda urgente alguien me puede enseñar o dar algun programa que realizace los metodos numericos como el metodode eliminación de Gauss, el metoodo de Gauus- Jordan , el metodo de Gauss Seidel, el metodo de matriz inversa, ya que me lo pidieron pero no tengo ni idea
se los agradesco

[johnny]

Hola...

Antes que nada devo suponer que conoces la teoria del metodo de Gaus_Jordan, o el metodo del Simplex,

Hace tiempo hice un programa que resolvia Sistemas de ecuaciones lineales de hasta 8 variables, que usaba el metodo del Simplex.

Ahí te envio el esqueleto del programa para que tú puedas implementar sus funciones.

El programa esta echo en el modo consola (MS_DOS).

Código: Seleccionar todo

///////////////////////////////////////////////////////////////////
//Programa que resuelve sistemas de ecuaciones lineales
#include <iostream.h>
#include <stdlib.h>
#include <string.h>
#include <iomanip.h>
#include <conio.h>


#define IMF 8         // Nº maximo de filas
#define IMC IMF+1      // Nº maximo de columnas

/*
      x1   x2   x3   b
      2   -5   6   12
      7   6   -5   15
      6   -2   1   -4
*/

char csino(char *pszmen)
{
   char cop;
   cout<<"\n\t\t"<<pszmen<<" (S/N): "<<flush;
   do{
      cop=toupper(getch());
   }while((cop!='S')&&(cop!='N'));
   putch(cop);
   cout<<endl;
   return cop;
}

// Declaramos dos matrices para almacenar los datos tipo float
static float fmatA[IMF][IMC],fmatB[IMF][IMC];

// Escribir los datos de la matrizA en la pantalla
static void vescribe_mat(int imaxf,int imaxc);

// Ingresar el nº_de_variables = nº_filas = nº_Columnas + 1
static void vent_parametros(int *pimf,int *pimc);

// Ingresar los datos en la matrizA
static void vlee_matriz(int imaxf,int imaxc);

// Transformar la matrizA, con pivote en (Fil= icont,Col= icont)
// Almacenando los datos transformados en la matrizB
static void vtrans_mat(int imaxf,int imaxc,int icont);

// Copiar los elementos de la matrizB en la matrizA
static void vcopia_mat(int imaxf,int imaxc);

// Imprime los resultados (la ultima columna de la matrizA)
static void vimp_result(int imxf,int imxc);

void main(void)
{
   int iop;   // Opcion ¿ continuar o no ?
   int imf;   // nº de filas
   int imc;   // nº de columnas
   int icont;   // contador de transformaciones
   do{
      icont=0;
      system("cls");

      // Inicializar las matrices
      memset(fmatA,0,sizeof(fmatA));
      memset(fmatB,0,sizeof(fmatB));

      // Ingresar los parametros Nº_Fil y Nº_Col
      vent_parametros(&imf,&imc);

      // Ingresar los datos (Cohef. y terminos independ.)
      vlee_matriz(imf,imc);

      // Escribir la matriz con los datos ingresados
      vescribe_mat(imf,imc);

      // Hacer mientras que el cont  <  al Nº_Filas
      while(icont<imf){
         vtrans_mat(imf,imc,icont);
         vcopia_mat(imf,imc);
         vescribe_mat(imf,imc);
         icont++;
      }
      vimp_result(imf,imc);
      iop=csino("Desea continuar");
   }while(iop=='S');
   return;
}


Espero haberte ayudado.
johnny

[eli_zabet]

hola te queria agradecerte por darme el codigo lo voy ha poner en practica y despues te aviso pero sinceramente muchas gracias.
chau, eli