Pascal/Turbo Pascal - SISTEMA DE ECUACIONES NXN POR GAUSS JORDAN

   
Vista:

SISTEMA DE ECUACIONES NXN POR GAUSS JORDAN

Publicado por luis luis_nando_8@hotmail.com (1 intervención) el 13/04/2016 18:48:59
Buen dia, necesito ayuda urgentemente para resolver un sistema de ecuaciones nxn por gauss jordan en lenguaje pascal(delphi).
quien pueda ayudarme con el codigo le agradeceria. es de vida o MUERTE, gracias.
Valora esta pregunta
Me gusta: Está pregunta es útil y esta claraNo me gusta: Está pregunta no esta clara o no es útil
0
Responder

SISTEMA DE ECUACIONES NXN POR GAUSS JORDAN

Publicado por ramon (2072 intervenciones) el 17/04/2016 18:05:02
Mira si esto te ayuda algo.

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
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
program gauss_jordan;
  {$N+}
  uses
     crt;
 
  function rnd: double;
  const
     bias = 1023;
  var
    data : record
            b : byte;
            d : double;
          end;
    x : array[0..8] of byte absolute data;
    e, i, j : word;
 
  begin
    for i := 0 to 7 do
    x[i] := lo(random(256));
    e := bias;
    repeat
      j := 0;
      for i := 0 to 7 do
      begin
        j := (x[i] shl 1) + hi(j);
        x[i] := lo(j);
      end;
      e := e - 1;
      if (bias - e) mod 8 = 0 then
      x[0] := lo(random(256));
    until (x[7] and $10) = $10;
    x[7] := (x[7] and $0F) or lo(e shl 4);
    x[8] := lo(e shr 4);
    rnd := data.d;
  end;
 
  function gauss(a,d: double): double;
  const
    t : double = 0;
  var
    v1, v2, r : double;
  begin
    if t = 0 then
    begin
      repeat
        v1 := 2 * rnd - 1;
        v2 := 2 * rnd - 1;
        r := v1 * v1  + v2 * v2
      until r < 1;
      r := sqrt((-2 * ln(r)) / r);
      t := v2 * r;
      gauss := a + v1 * r * d;
    end
    else
     begin
      gauss := a + t * d;
      t := 0;
    end;
  end;
 
  begin
    clrscr;
    writeln('   El Resultado Es');
    writeln;
    write('   ',gauss(4,4):0:10);
    readln;
  end.
Valora esta respuesta
Me gusta: Está respuesta es útil y esta claraNo me gusta: Está respuesta no esta clara o no es útil
0
Comentar