Метод Ньютона

[марианна]

Метод Ньютона

ребят, ну помогите еще с одним заданием пожалуйста

 

[Vladimir_S]

ребят, ну помогите еще с одним заданием пожалуйста

Пожалуйста.

Var
 x,y,b:Array[0..10] of Real;
 n,i,j,k:Byte;
 P,S,Q,z:Real;
Begin
 x[0]:=0.3240;
 for i:=1 to 10 do x[i]:=x[i-1]+0.1950;

 y[ 0]:=-0.3450;
 y[ 1]:=-0.0850;
 y[ 2]:= 0.0200;
 y[ 3]:= 0.0170;
 y[ 4]:=-0.0490;
 y[ 5]:=-0.1360;
 y[ 6]:=-0.1990;
 y[ 7]:=-0.1930;
 y[ 8]:=-0.0740;
 y[ 9]:= 0.2030;
 y[10]:= 0.6830;

 b[0]:=y[0];
 for i:=1 to 10 do
  begin
   S:=b[0];
   for j:=1 to i-1 do
    begin
     P:=1;
     for k:=0 to j-1 do P:=P*(x[i]-x[k]);
     S:=S+b[j]*P;
    end;
   P:=1;
   for k:=0 to i-1 do P:=P*(x[i]-x[k]);
   b[i]:=(y[i]-S)/P;
  end;

  z:=1.6960;
  Q:=b[0];
  for j:=1 to 10 do
   begin
    P:=1;
    for k:=0 to j-1 do P:=P*(z-x[k]);
    Q:=Q+b[j]*P;
   end;
  writeln('x= ',z:7:4,'    y= ',Q:7:4);

  z:=1.8860;
  Q:=b[0];
  for j:=1 to 10 do
   begin
    P:=1;
    for k:=0 to j-1 do P:=P*(z-x[k]);
    Q:=Q+b[j]*P;
   end;
  writeln('x= ',z:7:4,'    y= ',Q:7:4);

  z:=1.3400;
  Q:=b[0];
  for j:=1 to 10 do
   begin
    P:=1;
    for k:=0 to j-1 do P:=P*(z-x[k]);
    Q:=Q+b[j]*P;
   end;
  writeln('x= ',z:7:4,'    y= ',Q:7:4);
 Readln
End.

Несколько слов в пояснение.
Здесь, в отличие от темы про метод Лагранжа, я поленился писать заново и воспользовался своей старой программой отсюда:
http://www.tehnari.ru/f41/t85778/#post873309
Там же есть ссылки на описания самих алгоритмов, на базе которых писалась программа.

 

[марианна]

Пожалуйста.

Var
 x,y,b:Array[0..10] of Real;
 n,i,j,k:Byte;
 P,S,Q,z:Real;
Begin
 x[0]:=0.3240;
 for i:=1 to 10 do x[i]:=x[i-1]+0.1950;

 y[ 0]:=-0.3450;
 y[ 1]:=-0.0850;
 y[ 2]:= 0.0200;
 y[ 3]:= 0.0170;
 y[ 4]:=-0.0490;
 y[ 5]:=-0.1360;
 y[ 6]:=-0.1990;
 y[ 7]:=-0.1930;
 y[ 8]:=-0.0740;
 y[ 9]:= 0.2030;
 y[10]:= 0.6830;

 b[0]:=y[0];
 for i:=1 to 10 do
  begin
   S:=b[0];
   for j:=1 to i-1 do
    begin
     P:=1;
     for k:=0 to j-1 do P:=P*(x[i]-x[k]);
     S:=S+b[j]*P;
    end;
   P:=1;
   for k:=0 to i-1 do P:=P*(x[i]-x[k]);
   b[i]:=(y[i]-S)/P;
  end;

  z:=1.6960;
  Q:=b[0];
  for j:=1 to 10 do
   begin
    P:=1;
    for k:=0 to j-1 do P:=P*(z-x[k]);
    Q:=Q+b[j]*P;
   end;
  writeln('x= ',z:7:4,'    y= ',Q:7:4);

  z:=1.8860;
  Q:=b[0];
  for j:=1 to 10 do
   begin
    P:=1;
    for k:=0 to j-1 do P:=P*(z-x[k]);
    Q:=Q+b[j]*P;
   end;
  writeln('x= ',z:7:4,'    y= ',Q:7:4);

  z:=1.3400;
  Q:=b[0];
  for j:=1 to 10 do
   begin
    P:=1;
    for k:=0 to j-1 do P:=P*(z-x[k]);
    Q:=Q+b[j]*P;
   end;
  writeln('x= ',z:7:4,'    y= ',Q:7:4);
 Readln
End.

Несколько слов в пояснение.
Здесь, в отличие от темы про метод Лагранжа, я поленился писать заново и воспользовался своей старой программой отсюда:
http://www.tehnari.ru/f41/t85778/#post873309
Там же есть ссылки на описания самих алгоритмов, на базе которых писалась программа.

спасибо большое))))