Форум программистов, компьютерный форум, киберфорум
Наши страницы
Free Pascal
Войти
Регистрация
Восстановить пароль
 
Рейтинг 5.00/3: Рейтинг темы: голосов - 3, средняя оценка - 5.00
vladimir1092
0 / 0 / 0
Регистрация: 21.05.2013
Сообщений: 19
#1

Решение нелинейных уравнений методом Ньютона. Не выдает результат, когда решение = 0

26.12.2013, 22:07. Просмотров 588. Ответов 1
Метки нет (Все метки)

Есть программа,Решение нелинейных уравнений методом ньютона,но есть проблема,Не выдает результат, когда уравнение имеет корень 0!Помогите пожалуйста исправить проблему
Pascal
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
uses crt;
 
var a,d,t,n,z,p:integer;
i: byte;
mas: array [0..25] of integer;
x1,eps,x0,xn,h,q,w:real;
 
function f(x:real):real;
begin
f:= mas[1]+mas[2]*x+mas[3]*x*x+mas[4]*x*x*x+mas[5]*x*x*x*x+
mas[6]*x*x*x*x*x+mas[7]*x*x*x*x*x*x+mas[8]*x*x*x*x*x*x*x+mas[9]*x*x
*x*x*x*x*x*x+mas[10]*x*x*x*x*x*x*x*x*x+mas[11]*x*x*x*x*x*x*x*x*x*x+
mas[12]*x*x*x*x*x*x*x*x*x*x*x+mas[13]*x*x*x*x*x*x*x*x*x*x*x*x+
mas[14]*x*x*x*x*x*x*x*x*x*x*x*x*x+mas[15]*x*x*x*x*x*x*x*x*x*x*x*x*x*x;
end;
 
function f1(x:real):real;
begin
f1:=(mas[2])+(mas[3])*2*x+(mas[4])*3*x*x+(mas[5])*4*x*x*x+
(mas[6])*5*x*x*x*x+(mas[7])*6*x*x*x*x*x+(mas[8])*7*x*x*x*x*x*x+(mas[9])*8*x
*x*x*x*x*x*x+(mas[10])*9*x*x*x*x*x*x*x*x+(mas[11])*10*x*x*x*x*x*x*x*x*x+
(mas[12])*11*x*x*x*x*x*x*x*x*x*x+(mas[13])*12*x*x*x*x*x*x*x*x*x*x*x+
(mas[14])*13*x*x*x*x*x*x*x*x*x*x*x*x+(mas[15])*14*x*x*x*x*x*x*x*x*x*x*x*x*x;
end;
 
begin
clrscr;
 
writeln('Vvedite stepen polinoma:');
read(d);
i:=1; a:=d+1;
for t:=1 to a do
begin
z:=i-1;
writeln('Vvedite koeficient c[', z,']: ');
read(mas[i]);
i:=i+1;
end;
 
Writeln('Vvedite levuy granicu: ');
read(x0);
writeln('Vvedite pravuy granicu: ');
read(xn);
writeln('Vvedite shag: ');
read(h);
writeln('Tochnost:');
read(eps);
n:=1;
p:=0;
q:=x0; 
repeat
w:=q+h;
if f(q)*f(w)<=0 then
begin
writeln ('a= ',q:4:4, ' b= ',w:4:4);
repeat
x1:=q-f(q)/f1(q);
n:=n+1;
if abs(f(x1))>=eps then q:=x1;
if abs(f(x1))<=eps then p:=1;
if n>100 then p:=1;
until p=1;
writeln ('x= ',x1:3:3);
end;
q:=w;
until q>xn; writeln('konec poiska');
read(t);
end.
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
26.12.2013, 22:07
Ответы с готовыми решениями:

Решение нелинейных уравнений методом Ньютона
Решить уравнение методом Ньютона. y=log3(x+2)-1 На интервале x0=0 Соблюдая пункты: 1....

Решение нелинейных алгебраических уравнений методом простой итерации
Для нахождения корня используем бесконечный цикл с послеусловием и выходом из цикла с поиощью...

Приближенное решение нелинейных уравнений
На заданном отрезке найти все корни уравнения f(x)=0 с точность ε {e}^{-x}-3 sin x-0,4=0

Решение нелинейных уравнений методом Ньютона
Помогите решить 0,1sin(x)+x^2-1=0 С точностью Е=10^-4 Интервал изоляции Зарание спасибо!

Решение нелинейных уравнений методом Ньютона
решить методом касательных(метод ньютона). на промежутке (1;3) 2*(ln(x))^2+6*ln(x)-5

1
volvo
Супер-модератор
Эксперт Pascal/DelphiАвтор FAQ
25968 / 17412 / 6906
Регистрация: 22.10.2011
Сообщений: 30,656
Записей в блоге: 6
27.12.2013, 18:07 #2
Можно посмотреть какие значения ты вводишь в программу?
0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
27.12.2013, 18:07

Решение нелинейных уравнений методом Ньютона.
проверьте пожалуйста. непойму в чем ошибка program zadacha; uses Crt; var x0,x1,a,b,e:real; ...

Решение нелинейных уравнений методом Ньютона
Помогите...плиз...очень...очень...нужно..... уравнение: 2xsinx-cosx=0; Отрезок, содержащий...

Решение системы нелинейных уравнений методом Ньютона
Решить систему методом Ньютона с точностью ε= 0,00001. x=tg(xy)=x2 y=0.5*x2+2*y2=1 Написать...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
2
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Рейтинг@Mail.ru