Цитата:
Сообщение от felarl
Сам ход решения хоть правильный?
|
Не могу сказать, ибо хода не понимаю. Вы что, пытаетесь в одну программу оба метода засадить? Допустим. Но, во-первых, там F1 (первая производная) сосчитана неверно (откуда Вы взяли коэффициент 64? Должен быть 16), а во-вторых - ну на кой там вторая производная?
Цитата:
Сообщение от felarl
Да, да вы правы там +1 !!! Корень помоему где-то 0.14 должен быть. На счет того что не правильно задал функции, это да, просто так объяснили нам...
|
Совсем другое дело. Теперь график выглядит так:
Но корень там вовсе не "где-то 0.14", а около -0.1. Между прочим, в случае одного вещественного корня кубическое уравнение решается точно
методом Кардано, и это можно использовать для контроля правильности результата работы итерационных программ. Так вот, формула Кардано даёт х=-0.112085...
Теперь решим задачу методом Ньютона:
Код:
Const
e=0.001;
Var
x0_old,x0_new,D:Real;
Function F(z:real):real;
begin
F:=z*z*z*2-z*z*8+z*8+1;
end;
Function F1(z:real):real;
begin
F1:=z*z*6-z*16+8;
end;
Begin
x0_old:=0;
Repeat
x0_new:=x0_old-F(x0_old)/F1(x0_old);
D:=Abs(x0_old-x0_new);
x0_old:=x0_new;
Until D<e;
Writeln('Result: x = ',x0_new:0:3);
Readln
End.
Результат работы этой программы: х=-0.112.