LINUX.ORG.RU
ФорумTalks

Решение уравнения cos x = x

 ,


0

3

Долго думал сегодня над этим уравнением, нашел вот такое решение:

http://storage4.static.itmages.ru/i/15/0304/h_1425473189_2549951_8554cca376.png

x=\frac1\pi \int_0^{\pi } \arctan\left(\tan \left(\frac{t-\sin t+\frac{\pi }{2}}2\right)\right) \, dt+\frac{1}{\pi }

Проверьте, пожалуйста

★★★★★

Последнее исправление: Nxx (всего исправлений: 2)

По приказу командования «военный косинус» на множестве комплексных чисел может быть равен х при любом х. © :)

quickquest ★★★★★
()

Если честно лень вручную проверять, а axiom-а под рукой нет.

Но любопытно как получил такое выражение: про полиномы Чебышева вспомнил или как-то иначе пришел?

praseodim ★★★★★
()
Ответ на: комментарий от post-factum

Тот случай, когда если все делать строго численное решение для проверки не годится. Тут надо подставить x и провести преобразования.

praseodim ★★★★★
()
Ответ на: комментарий от Medar

И кто так тип float сравнивает?

Точно: if (1679 * 9993 != 1679.0f * 9993.0f) Console.WriteLine("Внезапно!"); )

Sadler ★★★
()
Последнее исправление: Sadler (всего исправлений: 1)
Ответ на: комментарий от praseodim

или потому что погрешность такая?

Да погрешность, конечно. Маткад никогда дальше 7 знака интегралы верно и не выводил на моей памяти.

Sadler ★★★
()
Ответ на: комментарий от praseodim

Интеграл посчитан неточно.

Nxx ★★★★★
() автор топика

Черт, прочитал топик как «os x = x».

DeadEye ★★★★★
()
Ответ на: комментарий от Nxx

Возможно, это потому, что я юзаю урезанную версию.

Sadler ★★★
()

Разложи cos x в ряд Тейлора.

invy ★★★★★
()
Ответ на: комментарий от morse

Арктангенс тангенса ведет себя как y=x для x<-(-Pi/2;Pi/2), и обладает периодом, равным Pi. Эдакая пила. Я пытался разбирать решение ТС, там интеграл берется в элементарных, если его в нужном месте разбить на два интеграла, но вот нахождение этого самого места - задача по сложности равная решению cos(x)=x.

keyran ★★
()

интеграл будет равен

pi-x, где x решение уравнения cos(x)=x.

Так что не совсем верно. Надо где-то коэффициенты поправить в интеграле...

dikiy ★★☆☆☆
()
Последнее исправление: dikiy (всего исправлений: 1)
Ответ на: комментарий от Nxx

да. только что хотел написать, что ошибся :) Минус потерял по ходу дела. Теперь все сходится: да, решение верное.

dikiy ★★☆☆☆
()
Ответ на: комментарий от Nxx

не знаю. Но меня смущает (n/2)^(2m+....) может (1/2)^(2m+....)?

dikiy ★★☆☆☆
()
Вы не можете добавлять комментарии в эту тему. Тема перемещена в архив.