Assasin, Вы направили начальную скорость почти на "большое тело", поэтому при сильном сближении увеличивается ускорение, выбранная временная дискретность (t = 0.1;) становится слишком велика и вычисления становятся очень далеки от действительности (и не говорите, что Вас не предупреждали ) (Кстати - что это за тело массой 4Е12; маловато будет, Земля=5,976E+24 кг, не говоря о Солнце - 1,991E+30 кг )
Поставте начальную скорость:
t1.v.x := -2;;
t1.v.y := 2; и получите вполне приличный эллипс (не забудте раскоментировать и приращение положения a.x * t*t / 2; )