ChrisAndre

import PyKEP
import random
import math

r_ang = random.uniform(0, math.pi * 2)
rx = math.cos(r_ang)
ry = math.sin(r_ang)
rz = 0

vx = -math.sin(r_ang)

ChrisAndre

import PyKEP as pk
import matplotlib as mpl
from mpl_toolkits.mplot3d import Axes3D
import matplotlib.pyplot as plt

r1 = [149.0e9, 0.0, 0.0]
r2 = [300.0e9, -100.0e9, 0.0]
mu = 1.32e20
tof = 86400 * 550.0
lw = True