kotejebekomet

import numpy as np
from matplotlib import *
from pylab import *
import math

print "komet"

x0=0.1
x=0.000000001
y=0.000000001
vx0=0.30
vy0=3.0
y0=10.0
vy=0.0

GM=-100.0

xx=[]
yy=[]

for i in np.arange(0, 10, 0.1):

    r=math.sqrt(((x0**2)+(y0**2)))
    x=(1.0/2.0)*(GM/(r**2))*x/abs(x)*0.1**2+vx0*0.1+x0
    x0=x
    vx=(GM/(r**2))*0.1*x/abs(x)+vx0
    vx0=vx

    y=(1.0/2.0)*(GM/(r**2))*y/abs(y)*0.1**2+vy0*0.1+y0
    vy=(GM/(r**2))*y/abs(y)*0.1+vy0
    y0=y
    vy0=vy

    print x, vx, y, vy
    xx.append(x)
    yy.append(y)

plot(xx, yy, "g.")
show()