from numpy *
N = 1000 # Number of iterations
r = zeros((N,3))# 3darray
v = zeros((N,3)) # 3darray

for i in range(N-1):
    F = ... # Forces between planet and sun
    for j in range(3):
        ... #Euler solver