У кого есть алгоритмы?

[Kirill Grishin]

А вот ещё один вопрос на засыпку: как избавиться от неопределённости при использовании арктангенса?
как правило рассчёты могут быть ошибочными из-за этого

[Kirill Grishin]

А вот ещё один вопрос на засыпку: как избавиться от неопределённости при использовании арктангенса?

  Ну, лорд Кирилл фон Гришин, Вы прямо разочаровываете такими вопросами...

Вы когда нибудь реализовывали алгоритмы, перед тем как резко переходить на личности?

[Kirill Grishin]

Вы когда нибудь реализовывали алгоритмы, перед тем как резко переходить на личности?

  Yes!

И как вы решили эту задачу?

[Kirill Grishin]

Ну не знаю, пишу параллельно языках на 5 разных для трёх разных аппаратных платформ и не имею особых сложностей. atan2 есть везде, но может отличаться немножко семантика.

Да, его можно встретить и в js и в Delphi
А вот ещё один интересный вопрос: как решается уравнение Кеплера для гиперболических орбит?

[Kirill Grishin]

А вот ещё один вопрос на засыпку: как избавиться от неопределённости при использовании арктангенса?
как правило рассчёты могут быть ошибочными из-за этого

Посмотри на страничке 49 в книжке Астрономия с калькулятором.
Внизу есть картинко, там всё расписано где какой знак арктангенса взависимости от
знаков  значений (x или y) стоящих под арктангенсом arctan(х/y), вычесленных зараннее.

atan2 не нужен. Ну твой вопрос хороший.

Да, но другое дело "научить программу определять нужный квадрат"

[Kirill Grishin]

Что касается e<1  то я уже давно написал код программы на дельфи:

Код: [Выделить]

  var m,e,ee,e1,m2:Extended;
  const rad=180/Pi;
begin
  m:=StrToFloat(edt1.Text);{средняя аномалия}
  e:=StrToFloat(edt2.Text);{эксцентриситет}
  e1:=m-e;{промежуток в котором лежит E}
  repeat
    e1:=e1+0.000001;{метод подстановки}
    m2:=e1-e*sin(e1/rad);{условие точности}
  until Abs(m2-m)<0.00001;{проверка условия точности}
  Label1.Caption:=FloatToStr(e1);
end;



[Kirill Grishin]

Пытался реализовать выше описанный алгоритм из книги Практическая астрономия с калькулятором, получается не тот ответ что нужен.
Acrtan2(x, y) оказался более полезным

[Kirill Grishin]

Негде не могу найти среднюю аномалию для планет(М) на эпоху J 2000.0
Есть устаревшие данные для J 1970
Может кто-нибудь подскажет ресурс