Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




 [Mathematica] Построение Re[y]
Аватара пользователя
Имеется уравнение F[x,y]==0, и хотелось бы построить график в координатах x, Re[y] где х принимает значения из некоторого известного промежутка. Я пыталась использовать ContourPlot, поскольку функция задана неявно, но ничего не получается.
Спасибо за помощь!

 Re: Построение Re[y]
можно попытаться решить это уравнение относительно у, а потом взять вещественную часть

 Re: Построение Re[y]
уравнение сильно нелинейное, я уже думал о том, чтобы решать его численно. Но увы из всех пересмотренных методов только метод мюллера сходится и то не при всех начальных условиях. из самой постановки задачи начальное приближение неизвестно.

 Re: Построение Re[y]
Аватара пользователя
volchenok (или MultiVac; я что-то запутался ;-), можете привести ваше уравнение, чтобы знать, от чего плясать?

 Re: Построение Re[y]
Аватара пользователя
Со слов автора:
Код:
a = Sqrt[x^2 - (1 - 1/y^2) y^2]
b = Sqrt[x^2 - 2 y^2]
c = Sqrt[x^2 - y^2]

F[x, y] := -E^(-2 a) (-1 + (c (1 - 1/y^2))/a) (-1 + (b (1 - 1/y^2))/(2 a)) + (1 + (c (1 - 1/y^2))/a) (1 + (b (1 - 1/y^2))/(2 a))

(Оффтоп)

MultiVac писал(а):
Я напишу не в синтаксисе LaTeX, а так, как оно в Mathematica выглядит, чтобы, если что Вам не пришлось перпечатывать всё в неё:
F[x, y] := -E^(-2 a) (-1 + (c (1 - 1/y^2))/a) (-1 + (b (1 - 1/y^2))/(2 a)) + (1 + (c (1 - 1/y^2))/a) (1 + (b (1 - 1/y^2))/(2 a))
где
a = Sqrt[x^2 - (1 - 1/y^2) y^2]
b = Sqrt[x^2 - 2 y^2]
c = Sqrt[x^2 - y^2]
Я её строю при помощи ContourPlot, и получается вот такое:
Изображение
А хотелось бы для наглядности продлить верхний график до начала координат, вот так:
Изображение

Может быть, нужно действительно писать уравнение на обе части аргумента, я даже не знаю, как это определить из функции.
Ещё я попробовала построить уравнение типа Re[F[x,y]]==0, но там получается ещё одна ветвь, которая, наверное, соответствует каким-то комбинациям, где оба аргумента мнимые, или что-то такое:
Изображение


Проблемы начинаются около нуля (пунктир от руки), когда мнимая часть «шалит», а действительная колеблется около нуля.

 [ Сообщений: 5 ] 


Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group