Решение системы нелинейных уравнений, Метод итераций и метод Зейделя... Опции

[Jill]

Помогите, ребята! Про линейные системы информацию нашла.
А вот про нелинейные - не нашла ни слова...

Задана система:
cos(x+0.5)-y=2
siny-2x=1

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

[klem4]

Хм вроде перевел с горем пополам, но что-то я не догоняю, программ практически при любых данных выдает один и тот-же результат, хотя возможно где-то я и ошибся :

Добавлено : гыгы надо просто ф-ю поменять ;)

uses crt;
const
   n = 2;
var
   f,x,w : array[1..n] of double;
   eps,d,e,g : double;
   i,k,m,fin,mis : integer;

procedure snur;
begin
    f[1] := sqrt(x[1]) - sqrt(x[2] -1) - x[1] + 8;
    f[2] := -sqrt(x[1]) - sqrt(9 - x[2]) - x[2] + 10;
end;

begin
   clrscr;
   write('Введите предельную погрешность : '); readln(eps);
   writeln('Введите начальное приближение');
   clrscr;
   for i := 1 to n do begin
       write('x[',i,']='); readln(x[i]);
   end;
   k := 0; d := 0.5; fin := 1; mis := 1;

   repeat
      inc(k);
      d := -d;
      gotoxy(4,60);
      writeln(k,' - e приближение');
      e := 0;
      for i := 1 to n do begin
         snur;
         x[i] := x[i] + f[i];
      end;
      for i := 1 to n do e := e + f[i] * f[i];
      e := sqrt(e);
      g := e; w[round(1.5 + d)] := e;
      gotoxy(5,60);
      write('Текущая невязка ');
      gotoxy(6,60);
      writeln(g:3:3);
      delay(100);
      if (w[round(1.5 - d)] - w[round(1.5 + d)] < 0) and (k > 9) then begin
         mis := 0;
         fin := 0;
      end;
      if e > eps then
       for i := 1 to n do x[i] := x[i] + f[i]
      else begin
         fin := 0;
         for i := 1 to n do x[i] := x[i] + eps;
      end;
      gotoxy(1, 1);
      for i := 1 to n do writeln('x[',i,']=',x[i]:3:2);
   until (fin = 0) or (k>1000);

   if mis = 1 then begin
      writeln('Окончательный результат : ');
      writeln('Точность удовлетворительная ');
   end
    else writeln('Процесс расходится ');
   readln;
end.

Сообщение отредактировано: volvo - 5.11.2006 21:27