Autor Tema: Programar el método de eliminación gaussiana en Mathematica

0 Usuarios y 1 Visitante están viendo este tema.

28 Noviembre, 2012, 04:39 pm
Leído 1472 veces

xamo

  • Aprendiz
  • Mensajes: 373
  • Karma: +0/-0
  • Sexo: Femenino
Buenas, me requieren varios apartados. El primero es:

1- Programar una función llamada solGaussSimple que resuelva un sistema de ecuaciones mediante el método de Gauss sin algoritmo de pivoteo. Esta función tomará dos parámetros: la matriz del sistema y el vector segundo miembro y devolverá un vector que represente la solución, x, del sistema de ecuaciones Ax = b.


El código que he escrito es:

solGaussSimple[A_, b_] :=
 Module[{m = 0, i, k, n = Length[A] },
  For[k = 1, k = n - 1, k++,
   For[i = k + 1, i = n, i++,
     m = -A[[i, k]]/A[[k, k]]
        A[] = A[] + m*A[[k]]
          b[] = b[] + m*A[[k]]
     ];
   ];
  Print[MatrixForm[A], MatrixForm[,b]];
  LinearSolve[A, b]]


sin embargo, no me hace las transformaciones necesarias. A ver si me podéis ayudar pronto, por favor.

PD: La coma en rojo es para que no salga un pedazo de texto en negrita.

28 Noviembre, 2012, 08:10 pm
Respuesta #1

xamo

  • Aprendiz
  • Mensajes: 373
  • Karma: +0/-0
  • Sexo: Femenino
Ya lo resolví. Gracias a todo el que se hubiera dignado en ayudar.  :D