Не выходит получить правильный ответ методом Ньютона (касательных)

@Katachresis · Регистрация: 18.06.2014

Студворк — интернет-сервис помощи студентам

Для уравления cos*(ln(x)-cos^2(x))-1⁄3=0 найти решения на отрезке [-5, -4].
Методом касательных в 15 итераций x0 = -4.9293.
Методом Ньютона уже на 2 итерации выдает x0 = -4.7548 + 0.0852i
Код программы под спойлером

Кликните здесь для просмотра всего текста

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
clear 
clc
xa = -5;
xb = -4;
y = @(x) cos(x).*(log(x)-(cos(x)).^2)-1/3;
dy = @(x) cos(x).*(2*cos(x).*sin(x) + 1./x) - sin(x).*(log(x) - cos(x).^2);
ddy = @(x) - cos(x)*(log(x) - cos(x)^2) - 2*sin(x)*(2*cos(x)*sin(x) + 1/x) - cos(x)*(2*sin(x)^2 - 2*cos(x)^2 + 1/x^2);
res = 1;
%метод Ньютона
d = @(x) y(x).*ddy(x);
if d(xa) > 0
    x0 = xa;
elseif d(xb) > 0
    x0 = xb;
else
    res = 0;
end
 
for i=1:15
    i
    x0 = x0 - y(x0)/dy(x0)
end
res
x0

Помогите найти ошибки в коде, пожалуйста

@SSC · 06.09.2016, 15:14

Сообщение от Katachresis

cos*(ln(x)-cos^2(x))-1⁄3=0

Как минимум после cos пропущено (x).
К тому же в 5 строке квадрат берется от разности log(x)-cos(x), а в функции в квадрате только cоs.
Уточняйте задачу

Добавлено через 6 минут
Кроме этого log из отрицательного числа только комплексный.

Добавлено через 1 минуту

Сообщение от Katachresis

Методом касательных в 15 итераций x0 = -4.9293.

И как Вам удалось получить действительный корень?

Новые блоги и статьи Все статьи Все блоги /
Новый CodeBlocs. Версия 25.03 palva 04.01.2026 Оказывается, недавно вышла новая версия CodeBlocks за номером 25. 03. Когда-то давно я возился с только что вышедшей тогда версией 20. 03. С тех пор я давно снёс всё с компьютера и забыл. Теперь. . .	Модель микоризы: классовый агентный подход anaschu 02.01.2026 Раньше это было два гриба и бактерия. Теперь три гриба, растение. И на уровне агентов добавится между грибами или бактериями взаимодействий. До того я пробовал подход через многомерные массивы,. . .	Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Programma_Boinc 28.12.2025 Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Налог на собак: https:/ / ********/ gallery/ V06K53e Финансовый отчет в Excel: https:/ / ********/ gallery/ bKBkQFf Пост отсюда. . .	Кто-нибудь знает, где можно бесплатно получить настольный компьютер или ноутбук? США. Programma_Boinc 26.12.2025 Нашел на реддите интересную статью под названием Anyone know where to get a free Desktop or Laptop? Ниже её машинный перевод. После долгих разбирательств я наконец-то вернула себе. . .	Thinkpad X220 Tablet — это лучший бюджетный ноутбук для учёбы, точка. Programma_Boinc 23.12.2025 Рецензия / Мнение/ Перевод Нашел на реддите интересную статью под названием The Thinkpad X220 Tablet is the best budget school laptop period . Ниже её машинный перевод. Thinkpad X220 Tablet —. . .
PhpStorm 2025.3: WSL Terminal всегда стартует в ~ and_y87 14.12.2025 PhpStorm 2025. 3: WSL Terminal всегда стартует в ~ (home), игнорируя директорию проекта Симптом: После обновления до PhpStorm 2025. 3 встроенный терминал WSL открывается в домашней директории. . .	Как объединить две одинаковые БД Access с разными данными VikBal 11.12.2025 Помогите пожалуйста !! Как объединить 2 одинаковые БД Access с разными данными.	Новый ноутбук volvo 07.12.2025 Всем привет. По скидке в "черную пятницу" взял себе новый ноутбук Lenovo ThinkBook 16 G7 на Амазоне: Ryzen 5 7533HS 64 Gb DDR5 1Tb NVMe 16" Full HD Display Win11 Pro	Музыка, написанная Искусственным Интеллектом volvo 04.12.2025 Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .	От async/await к виртуальным потокам в Python IndentationError 23.11.2025 Армин Ронахер поставил под сомнение async/ await. Создатель Flask заявляет: цветные функции - провал, виртуальные потоки - решение. Не threading-динозавры, а новое поколение лёгких потоков. Откат?. . .