Решение СЛАУ методом Монте-Карло Опции

[MagicPRO]

Народ у меня расчетно-графическая робота, решение системы линейных уравнени методом Монте-Карло!
Может кто знает или у кого-то есть код решения этой задачи, принимается код на С++!
Заранее благодарен!

[klem4]

В поиске смотрел ?

Монте-Карло

Сообщение отредактировано: klem4 - 6.06.2007 17:58

[MagicPRO]

Помогите пределать код с Basica на Pascal!

10 REM ~~~МОНТЕ_КАРЛО~~~
20 REM ВЫЧИСЛЕНИЯ КОРНЯ СИСТЕМЫ ЛИНЕЙНЫХ
УРАВНЕНИЙ
30 DIM A(20,2l),B(20,20),W(20)
40 INPUT "ВВЕДИТЕ ПОРЯДОК РЕШАЕМОЙ СИСТЕМЫ";N
50 PRINT "ВВЕДИТЕ ПО СТРОКАМ КОЭФФИЦИЕНТЫ
СИСТЕМЫ"
60 FOR I=1 ТО N
70 INPUT A(I,I):B(I,I)=ABS(A(I,I)) 80 FOR J=2 TO N
90 INPUT A(I,J):B(I,J)=B(I,J-I)+ABS(A(I,J))
100 NEXT J
110 INPUT A(I,N+1)
120 NEXT I
130 FOR I=1 TO N
140 W(I)=A(I,N+1)/(1-B(I,N))
150 NEXT I
160 INPUT "ВВЕДИТЕ НОМЕР РЕШАЕМОГО УРАВНЕНИЯ -В И КОЛИЧЕСТВО РЕАЛИЗАЦИИ СЛУЧАЙНОГО ПРОЦЕССА - М ";В,М
170 T=1:Y=0
180 S=B:Y=1
190 IF T>M THEN 260
200 C=RND(1)
210 FOR J=N TO 1 STEP -1
220 IF C<=B(S,J) THEN 240
230 IF J=N THEN T=T+1:Y=Y+V*W(S):GOTO 180 ELSE V=V*SNG(A(S,J+1)):S=J+1:GOTO 200
240 NEXT J
250 V=V*SGN(A(S,l)):S=1:GOTO 200
260 X=Y/M
270 PRINT "КОРЕНЬ ";B;" - ГО УРАВНЕНИЯ X=";X
280 END

[MagicPRO]

Ура, я код переделал, но теперь есть ошибки с вводом, помогите их решить!
я думаю не одному мне этот метод понадобится, в инете я все обискал, только одну прогу платную нашел!

uses crt;

  function SGN(a: real): integer;
  begin
    if a > 0 then
      SGN := 1;
    if a = 0 then
      SGN := 0;
    if a < 0 then
      SGN := -1;
  end;

var
  a:    array[1..20, 1..21] of real;
  b:    array[1..20, 1..20] of real;
  w:    array[1..20] of real;
  n, d, m, t, s: byte;
  i, j: word;
  x, c, v, y: real;
label
  label1, label2, label3;
begin
  v := 0;
  writeln('vvedite poradok system');
  readln(n);
  writeln('vvedite coeffs');
  for I := 1 to N do
  begin
    readln(A[I, I]);
    B[I, I] := ABS(A[I, I]);
    for J := 2 to N do
    begin
      readln(A[I, J]);
      B[I, J] := B[I, J - I] + ABS(A[I, J]);
    end;
    readln(A[I, N + 1]);
  end;
  readln;
  for I := 1 to N do
    W[i] := A[I, N + 1] / (1 - B[I, N]);
  writeln('vvedite nomer uravnenia');
  readln(d, m);
  T   := 1;
  Y   := 0;
  label2:
    S := d;
  Y   := 1;
  if T > M then
    goto label1;
  label3:
    C := random(1);
  for J := N downto 1 do
  begin
    if C <= B[S, J] then
      continue;
    if J = N then
    begin
      T := T + 1;
      Y := Y + V * W[S];
      goto label2;
    end
    else
    begin
      V := V * SGN(A[S, J + 1]);
      S := J + 1;
      goto label3;
    end;
  end;
  V := V * SGN(A[S, 1]);
  S := 1;
  goto label2;
  label1:
    X := Y / M;
  writeln('coren    ', d, ' - go uravnenia X=', X);
  readln;
end.