Показать сообщение отдельно
Старый 10.01.2019, 12:02   #2 (permalink)
xamelione25
Новичок
 
Регистрация: 09.01.2019
Сообщений: 4
Сказал(а) спасибо: 0
Поблагодарили 0 раз(а) в 0 сообщениях
Репутация: 10
По умолчанию

Я сделал все как в посте №2... но у меня вылетели почему-то не правильные ответы

Код:
Const
 e=0.0001;
Var
 x_old, x_new, y, z:real;

Function Fun_1(x:real):real;
begin
 Fun_1:=(2*sqrt(5/(1-x)))+x-1;
end;

Function Fun_2(x:real):real;
begin
 Fun_2:=(-2*sqrt(5/(1-x)))+x-1;
end;

Function Der_1(x:real):real;
begin
 Der_1:=sqrt(5)*power((1/(1-x)), (3/2))+1;
end;

Function Der_2(x:real):real;
begin
 Der_2:=1-sqrt(5)*power((1/(1-x)), (3/2));
end;

Begin    
 x_old:=-1.5;
 Repeat
  z:=x_old;
  x_new:=x_old-Fun_1(x_old)/Der_1(x_old);
  x_old:=x_new;
 Until Abs(x_old-z)<e;
 y:=-Sqrt(5*(1-x_old));
 Writeln('x1= ',x_old:0:5);
 Writeln('y1= ',y:0:5);
 Writeln;
 x_old:=1.5;
 Repeat
  z:=x_old;
  x_new:=x_old-Fun_2(x_old)/Der_2(x_old);
  x_old:=x_new;
 Until Abs(x_old-z)<e;
 y:=Sqrt(5*(1-x_old));
 Writeln('x2= ',x_old:0:5);
 Writeln('y2= ',y:0:5);
 Readln;
End.
Изображения
  
xamelione25 вне форума   Ответить с цитированием
Ads

Яндекс

Member
 
Регистрация: 31.10.2006
Сообщений: 40200
Записей в дневнике: 0
Сказал(а) спасибо: 0
Поблагодарили 0 раз(а) в 0 сообщениях
Репутация: 55070