Вычислить значение выражения

var F,L:real;
    i,j,n:integer;
    x:array[1..10] of real;
    y:array[1..10] of real;
 
begin
     write('n=');readln(n);
     FOR i:=1 TO n DO
     begin
          write('x[',i,']=');readln(x[i]);
          write('y[',i,']=');readln(y[i]);
     end;
     begin
          write('x[',n+1,']=');readln(x[n+1]);
     end;
     y[n+1]:=0;
     F:=0;
     FOR j:=1 TO n DO begin
     L:=1;
     FOR i:=1 TO n DO
     begin
         IF i<>j THEN
         begin
         L:=L*(x[n+1]-x[i])/(x[j]-x[i]);
         end;
     end;
     y[n+1]:=y[n+1]+y[J]*L;end;
     writeln('y[',n+1,']=',y[n+1]:1:0);
     FOR i:=1 TO n DO
     begin
          writeln('x[',i,']=',x[i]:10:10,' y[',i,']=',y[i]:10:10);
     end;
     begin
          writeln('x[',n+1,']=',x[n+1]:10:10,' y[',n+1,']=',y[n+1]:10:10);
     end;
     readln;
end.

{
f(x)=a*x1^2+b*x2^2->min
g(x)=c*x1+d*x2+e=0
}
program L;
 
const
  a = 1;
  b = 1;
  c = 1;
  d = 1;
  e = -2;
 
  function f(x1, x2: real): real;
  begin
    f := a * x1 * x1 + b * x2 * x2;
  end;
 
var
  x1, x2, La: real;
begin
  {после составления системы линейных уравнений решим её}
  La := 2 * a * b * e / (a * d * d + b * c * c);
  x1 := -c * La / (2 * a);
  x2 := -d * La / (2 * b);
  writeln('Решив систему уравнений получим:');
  writeln('Lambda=', La);
  writeln('x1=', x1);
  writeln('x2=', x2);
  {Значение целевой функции в точке условного экстремума}
  writeln('Значение целевой функции в точке условного экстремума:');
  writeln('f(', x1, ', ', x2, ')=', f(x1, x2));
end.