Метод Ньютона для решения системы m нелинейных уравнений

[Viktoria_S]

Метод Ньютона для решения системы m нелинейных уравнений

Здравствуйте. Помогите, пожалуйста, с созданием программы на языке Pascal в среде lazarus. На Вашем форуме нашла пример решения подобной задачи:

Var
Eps,Res:Real;
i:Integer;

Function F(z:real):Real;
begin
F:=Sqr(Sin(z)+Cos(z))/Exp(Ln(33.5)*2/3)+Sqrt(3/7)-z;
end;

Function G(z:real):Real;
begin
G:=(Sin(z)+Cos(z))*2*(Cos(z)-Sin(z))/Exp(Ln(33.5)*2/3)-1;
end;

Procedure EQRoot(e:real; var R:real);
var
x1,x2,d:real;
begin
x1:=0.0;
Repeat
x2:=x1-F(x1)/G(x1);
d:=abs(x1-x2);
x1:=x2;
Until d<e;
R:=x2;
end;

Begin
Eps:=0.1;
for i:=1 to 10 do
begin
EQRoot(Eps,Res);
Writeln(Res:12:10);
Eps:=Eps/10;
end;
Readln
End.


Но здесь, как я понимаю, решения одного нелинейного уравнения, а мне нужна система уравнений. Подскажите, как это реализовать? Заранее спасибо.

 

[Vladimir_S]

Но здесь, как я понимаю, решения одного нелинейного уравнения, а мне нужна система уравнений. Подскажите, как это реализовать? Заранее спасибо.

Так, как расписано здесь: Метод Ньютона — Википедия в разделе "Многомерный случай".