Planetary Motion,Feynman Lectures

[@run]

Homework Statement

This question is a part of independent study.
The problem is based on Feynman lectures volume 1,chapter 9,Newton's laws of dynamics,section 9-7 planetary motion. To compute numerically path of motion of a planet around sun using Newtons laws of kinematics and gravitation Feynman took initial conditions as x=0.5,y=0,v(x)=0,v(y)=1.63 & e=0.1.I understand that v(x)=0 and v(y) should have some value initially,but cannot understand why he take the value 1.63 for v(y) and how he computed 1.63.

Homework Equations

F = -G*M*m/r^2
r=sqrt(x^2+y^2)
x=x(0) + e* v(x) , v(x) = v(0) + e*a(x) , a(x) = -x/(r^3)
y=y(0) + e* v(y) , v(y) = v(0) + e*a(y) , a(y) = -x/(r^3)

The Attempt at a Solution

Two attempts are made
1
.x = a cos $, y = b sin $ -------(1)
v(x) = -a(d$/dt) sin $ -------(2),
v(y) = -b(d$/dt) cos $ = -(bx/a)(d$/dt)--------(3)
y = 0 means $ = 0, x=a.
so, v(y) = -b(d$/dt) = -bw
now there is a relation connecting the parameters a,b and e of an ellipse
b*b = a*a(1-e*e).
taking e = 0.8(a guss) value of b is calculated.
If I can calculate w in some way, i can get v(y). Some attempts made using the w*r^2 is conserved.
But couldn't calculate w.
2.
This attempt is made first for understanding that v(y) = 1.63 have any importance.
A simple program in c++ is made and plotted x-y graph.This analysis revealed that 1.63 have importance,then only perfect ellipse will be formed. If v(y) is taken as 0.8 planet will not revolve
around sun,the same thing happens when a value greater than 1.63 is taken ,for example 2.2
Analyis also says that the value of v(y) depend on x.When i take x=1 ,v(y) should have a value
near 1.4.This can be expect naturally because when x value changes from 0.5 to 1 the
gravitational force decreased,so less v(y) required.

 

[member 553083]

So it can be concluded that v(y) depend on initial x value. But how feynman computed the value 1.63 remain a mystery.