Для нахождения

обычно решают систему диф. уравнений (4) и (3а).
Для этого решаем (4) относительно

:
(5)Уравнение (3а) в обозначениях

:
(6)Подставляем

из
(5) в (6) и решаем диф. уравнение относительно


Получаем :
(8)где A - интеграл движения.
Решаем квадратное уравнение (8) относительно

и имеем 2 решения:

Тогда соответственно 2 выражения для

:
(9)
(10)Координатная скорость :Найдем координатную скорость для двух радиальных геодезических:

Получаем после упрощений:
(11а)
(12а)Собственно эти два выражения нам в дальнейшем и понадобятся. Если они верные, то далее можно попытаться понять, что творится на горизонте.
Проведем оценку того, что творится с интервалом

для частиц двигающихся радиально , для этого подставим

из (11а) и (12а) в первоначальные выражение (1) метрики.
Получим:
(18)Это убегающая геодезическая . И вторая "падающая":
(19)На горизонте: при

стремящемся к

:
для первой и :
для второй.
Поскольку времениподобная не может обратиться в изотропную , для второй геодезической существует некая особенность на горизонте. Чтобы понять это, рассмотрю также как и со Шварцшильдом два частных решения.
-- 13.08.2013, 11:46 --Частные решения.3. A=1.Физический смысл постоянной
A я не буду пока рассматривать. Дело в том, что Новиков-Фролов ссылаются на Богородского, но у него, и в подобных решениях у Брумберга , у Гильберта, фигурирует просто интеграл движения А. Почему у них появилось вместо этого

, пусть остается на совести авторов. Скорее надо было сослаться на Ландау-Лифшица. При переходе к метрике Шварцшильда , у меня (может случайно) получилось , что

. Но я не буду этим вопросом озабачиваться.
Итак, пусть
A=1. Тогда Первая (падающая) геодезическая:

Вторая:

После интегрирования получим 2 геодезические на одном рисунке в координатах

:

рис. 3. (A=1)

, красная - "падающая", синяя - "Убегающая" геодезические в координатах Эддингтона-Финкельштейна. (вместо

иметь в виду

)
Видно , что красная-"падающая" спокойно пересекает горизонт, в то время как вторая " убегающая" (синяя)
расщепилась на 2 части. Собственно этот рисунок нам потребуется для дальнейшего анализа.
Вне горизонта частица попадает на правую синюю ветку и может улететь совсем, если ее скорость близка к скорости света. Внутри горизонта она падает по любому в центр к истинной сингулярности.
4. Построим также графики для другого случая

Приведу окончательные графики:

рис. 4

,

Начальные условия выбраны так, чтобы точка остановки для второй геодезической

была при

.
Вывод-вопросы.1.Таким образом, особенность на горизонте осталась. Это особенность физическая, а не координатная.В точке

скорость второй геодезической равна скорости света.
2. Можно ли корректно поставить задачу Коши на горизонте? Наверное нет. Это неустойчивое решение.
3. Приходит в голову такая задача: изотоп распадается на горизонте ( в случае вечной ЧД). Один осколок летит радиально в центр , а второй в зависимости от своего положения: чуть выше горизонта - летит по одной траектории ( сначала синяя , потом - красная), чуть ниже горизонта - по нижней синей ветке на рис. 4.
4. Качественно картина будет похожей и для другой метрики , например для Леметра, поскольку они связаны допустимыми координатными преобразованиями.
PS. Черновики залил сюда в формате docx:
http://yadi.sk/d/7MIoEz967rYUmНе знаю как приклеить файл.