Нужна помошь в создании програмы для решения систем нелинейных уровнений (матрицей Якоби) методом Ньютона на Delphi

Не практиковала язык Pascal поэтому не могу решить задачу. Дано матрицу Якоби которую нужно решить с помощью метода Ньютона. Есть только блок схема и процедура решения матрицы.

`Procedure JAC (Var Ja: Matrix; X: Vec; Q: Real; N: Integer);
Var i, j :Integer;
Fp : Vec; 
Begin
Q := Eps / 1000.0;
FM (F, X, N); 
For j:=1 to N do
Begin
X[j] := X[j] + Q; 
FM (Fp, X, N);
For i:=1 to N do
Ja[i, J] := (Fp[i] - F[i]) / Q;
X[j] := X[j] – Q 
End
End;
`
`Procedure FM( Var F: Vec; X:Vec; N: Integer);
Begin
F[1]:=x[1] + Exp (x[1]-1)+(x[2]+x[3])*(x[2]+x[3])-27;
F[2]:=x[1] * Exp (x[2] - 2) + x[3] * x[3] - 10;
F[3]:=x[3] + Sin (x[2] - 2) + x[2] * x[2] - 7
End;
`