Движение окружности по трактрисе

[No_Comments]

Привет, форумчане!
Нужна помощь!
Необходимо визуализировать траекторию движения окружности по трактрисе!
ПРоблем с реализацией на ЯВУ нет.
Есть проблема с формулами
Может кто подскажет, по каким формулам координаты рассчитываются?!

[Ezhik Kihze]

В вики же формула есть... Допустим по трактрисе двигается центр, стерли старую окружность, вычислили новый центр и нарисовали.

[No_Comments]

А эта формула правильная?!

[Ezhik Kihze]

Сомневаетесь? В гугле куча результатов с формулой - сравните, откройте учебник...

[No_Comments]

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

 public partial class Form1 : Form
    {
        private int x1, y1, x2, y2;
        private double a, t, fi;
        private Pen pen = new Pen(Color.DarkRed, 2);

        public Form1()
        {
            InitializeComponent();
        }

        private void Form1_Paint(object sender, PaintEventArgs e)
        {
            Graphics g = e.Graphics;
            g.DrawEllipse(pen, x2, y2,20,20);
        }

        private void timer1_Tick(object sender, EventArgs e)
        {
            fi += 0.01;
            t = Math.Tan(fi);
            x2 =(int) (a * (Math.Log(t, Math.E) + Math.Cos(fi)));
            y2 = (int)(a * Math.Sin(fi));
            Invalidate(); 

        }

        private void Form1_Load(object sender, EventArgs e)
        {
            x1 = ClientSize.Width / 2;
            y1 = ClientSize.Height / 2;
            a = 150;
            fi = -0.5;
            t = Math.Tan(fi);
            x2 = (int)(a*(Math.Log(t, Math.E)+Math.Cos(fi))) ;
            y2 = (int)(a*Math.Sin(fi));

        }

[No_Comments]

Вот параметрическое описание трактрисы

[phomm]

Начните с того, что tan(-0.5) отрицателен и log от него - NaN, в итоге падают эксепшены.

[No_Comments]

Спасибо за совет
Теперь на форме появляется окружность и движется от левого края формы к правой.
Но, когда окружность пропадает за правым краем формы, через какое-то время появляется на форме красный крест, как на этой картинке.


Перед вызовом Invalidate() поставила условие if(fi < Math.PI) чтобы избежать вызова метода, т.к. tg(PI)=0... Думала, поможет.