интерполяция лагранжа

Ответ от Lika Snegina[новичек]
{Вычисление значения многочлена Лагранжа для функции, определенной таблично} uses Crt; const m=9; type Vector=array [0..m] of real; var X,Y: vector; {узлы интерполяции} i,j,N: integer; p,L,Xe,h: real; function step: real; {определение оптимального шага интерполяции} var i: integer; min: real; begin min:=Abs(X[1]-X[0]); for i:=1 to N-1 do begin if Abs(X[i+1]-X)<min>) end; step:=min; end; {основное тело программы} begin TextBackground(7); TextColor(0); ClrScr; {очистка экрана}; Writeln(' ************************************************** *****'); Writeln(' * Программа интерполяции функции многочленом Лагранжа *'); Writeln(' ************************************************** *****'); Write('Введите порядок многочлена Лагранжа N (>0): '); Read(N); WriteLn('Введите пары значений Xi,Yi (i=0,...N): '); for i:=0 to N do begin Write('X[',i,']= '); Read(X); Write('Y[',i,']= '); Read(Y); end; h:=step/2; {Вычисление шага интерполяции} WriteLn('Шаг интерполяции h=', h:0:7); WriteLn('Нажмите любую клавишу для вывода интерполяционной таблицы... '); readkey; ClrScr; WriteLn(' Данные интерполяции: '); WriteLn(' X Y '); WriteLn(' ---------------------------'); Xe:=X[0]; while Xe<=(X[N]+h) do begin L:=0; for i:=0 to N do begin p:=1; for j:=0 to N do if j<>i then p:=p*(Xe-X[j])/(X-X[j]); L:=L+y*p; end; WriteLn(' ':10,Xe:0:3,' ':4,L:12:10); Xe:=Xe+h; end; WriteLn(' Шаг интерполяции h=', h:0:7); WriteLn; WriteLn('Нажмите любую клавишу для выхода из программы... '); readkey; end.