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

[Dark_san]

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

[Dark_san]

Нужна помощь, скрипя зубами, пытаюсь написать процедуру для Ньютона, но нужен хелп вот в ожном кусочке кода.

 var
  s, s1, p, p1, xx, H, J1, a: real;
  k, i, n, j, ii, e1, f1, n1: integer;
  x, y, z, z1: 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
            goto l;
          p:=p*(xx-x[j]);
          l:
          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;

____________________________________________

if j=i then
            goto l;
          p:=p*(xx-x[j]);
          l:

Вот кусочек когда, можно как то это оформить без goto? Меня препод съест

Сообщение отредактировано: Lapp - 21.04.2010 3:02

[Lapp]

if j=i then
            goto l;
          p:=p*(xx-x[j]);
          l:

Вот кусочек когда, можно как то это оформить без goto?

Гм?..

if j<>i then p:=p*(xx-x[j]);

Меня препод съест

И правильно сделает, приятного ему аппетита..

[Dark_san]

Гм?..

if j<>i then p:=p*(xx-x[j]);

И правильно сделает, приятного ему аппетита..

Т.е

if j=i then
            goto l;
          p:=p*(xx-x[j]);
          l:
          p1:=1;

можно заменить на

if j<>i then p:=p*(xx-x[j]) else
 p1:=1;

[Lapp]

Т.е

if j=i then
            goto l;
          p:=p*(xx-x[j]);
          l:
          p1:=1;

можно заменить на

if j<>i then p:=p*(xx-x[j]) else
 p1:=1;

Ни в коем разе!! Ты не передергивай - я "else" не писал!
В твоем посте был фрагмент с обходом оператора p:=p*(xx-x[j]) в случае, когда i=j. Все. Никаких других условий не было. То, что стоит за меткой l, будет выполняться всегда. Значит, никаких else. Ты уж определись, что именно тебе нужно, дорогая негритяночка..

Я посмотрю попозже твою прогу сегодня. И, может, напишу ченть. Мне нравится, что ты хоть что-то делаешь. Скажи препу - пусть оставит мне кусочек от обеда )). А вообще - ужасные у вас порядки..