Интерполяция многочленом ньютона Опции

[Dark_san]

Вот такое вот у меня было лицо
Нужен алгоритм для решения програмным путём
Входные данные - 12 X|Y и 4 точки x1 x2 x3 x4.

[Dark_san]

Спасибо но..
Во первых выдает Eror 207 Invalid Floating point operation.
Во вторых или я не выспалась, или я не вижу где используется массив с Y. И чуть чуть не поняла твоих констант...


Добавлено через 10 мин.

 Program Nyuton;
uses crt;
var
		str: string;
		s, s1, p, p1, xx,  a: real;
		k, i, n, j, ii, e1, f1, n1: integer;
		x, y, z : array [1..50] of real;



begin
  repeat
    clrscr;
    write('Vvedite kolichestvo tochek: ');
    readln(n);
    if n<=0 then
      begin
        writeln;
        write('Vvodit nado polozhitelnoe chislo! Nazhmite Enter');
        readln;
      end;
  until n>0;
		writeln;
  for i:=1 to n do
    begin
      Write('vvedite x',i,': '); Readln(x[i]);
      Write('vvedite y',i,': '); Readln(y[i]);
      writeln;
    end;
	write('  Vvedite x*:'); Readln(xx);
	s:=y[1];
  for i:=2 to n do
    begin
      p:=1; s1:=0;
      for j:=1 to i do
        begin
          if j<>i then p:=p*(xx-x[j]);
          
          p1:=1;
          for k:=1 to i do
            begin
              if j<>k then
                p1:=p1*(x[j]-x[k]);
            end;
          p1:=y[j]/p1; s1:=s1+p1;
        end;
      p:=p*s1; s:=s+p;
    end;
  writeln;
  for i:=1 to n do
	begin
			writeln('X[',i,']',x[i]4,'  Y[',i,']',y[i]4);
	end;
		write('  Otvet f(x*)=');
		textcolor(18);
		writeln(s:5:3);
		writeln;
		textcolor(15);
		write('  Nazhmite Enter');
  
  readln;
  end

Воооооооо ^_^ все оки теперь.
Спасиб ))) скоро прийду с схемой ейткина)))

[Lapp]

Во первых выдает Eror 207 Invalid Floating point operation.
Во вторых или я не выспалась, или я не вижу где используется массив с Y. И чуть чуть не поняла твоих констант...

Используется? В конечных разностях, ессно, по формуле. То бишь в функции FinDif. А задается он константой (я не стал вводить точно, только два знака). Это для простоты, поскольку не имеет отношения к главному. Если тебе нужно вводить вручную или из файла - сделай, как надо.

Теперь скажи: в какой строке ошибка? У меня нет Турбо, я использую только FreePascal. В нем все проходит, вроде.. Сейчас еще раз проверю.. Да, все идет и на этом компе. Вот ответ (с приведенными данными):
0.65459

Откомпилируй с информацией о строках и скажи, в какой строке сбой.

Добавлено через 6 мин.

Воооооооо ^_^ все оки теперь.

Фиии! "Nyuton"... кашмар..

Спасиб ))) скоро прийду с схемой ейткина)))

ну, это уже не ко мне.. мою прекрасную программу забраковала - успехов тебе..

[Dark_san]

0.65459

Откомпилируй с информацией о строках и скажи, в какой строке сбой.

Натупил..кхм... это как ? точнее где)

Добавлено через 2 мин.
Вобщем вот тут, на n-ой итерации p:=p+Prod(t,k-1)*FinDif(k,0)/f выбивает, Eror 207 Invalid Floating point operation. т.е переполнение... вот)