[Visual Python] View of Jupiter from Earth

As the title says I am attempting to build a view of Jupiter as it would appear if I looked at the night sky every night for a given time frame using the North Star (Polaris) as a reference point.

I am posting the current code I have utilized. This is a work in progress (I am not good at Math). Based on the code from here – http://www.phyast.pitt.edu/~akl2/Phys0174/files/EarthSun.py

from visual import *
import random


#Constants
G = 6.7e-11
day = 24*60*60
year = 365*day
pscale = 1e11/1e29/3
fscale = 1e11/1e22/4



#set the window length and center

#scene = display(width = 800, height = 800, center = (0,5,0))
Sun = sphere(pos=(0,0,0), radius=10e9, color=color.yellow)

#sphere objects
Star = sphere(
    pos=(0,1.9e11,0), 
    radius=3e9, 
    materials = materials.rough,
    color = color.white
    )

Earth = sphere(
    pos = vector(1.5e11,0,0),
    radius = 5e9,
    color=color.cyan
    )

Jupiter = sphere(
    pos=(1.6e11,0,0), 
    radius = 5e9, 
    color = color.orange,
    make_trail = true, 
    )


Starlabel = label(
    pos=Star.pos,
    text='North Star', 
    xoffset=20,
    yoffset=12, 
    height=10, 
    border=6,
    font='sans'
    )



parr = arrow(color=color.red)
farr = arrow(color=color.green)

Sun.m = 2e30
Earth.m = 6e24
Jupiter.m = 1.89813e+27

#Initial values
Sun.p = vector(0,0,0)
Earth.p = Earth.m*vector(0,2*pi*Earth.pos.x/year,0)
Earth.trail = curve(color=Earth.color)
Earth.trail.append(pos=Earth.pos)

Jupiter.p = Jupiter.m*vector(0,2*pi*Jupiter.pos.x/year,0)
Jupiter.trail = curve(color=Jupiter.color)
Jupiter.trail.append(pos=Jupiter.pos)


deltat = day
t=0
count = 0

for i in range(365):
    rate(50)
    r = Earth.pos - Sun.pos
    r1 = Jupiter.pos - Sun.pos
    rmag = mag(r)
    rmag1 = mag(r1)
    rhat = r/rmag
    rhat1 = r1/rmag1
    Fnet = (-G*Earth.m*Sun.m/rmag**2)*rhat
    Fnet1 = (-G*Jupiter.m*Sun.m/rmag1**2)*rhat1
    
    Earth.p = Earth.p + Fnet*deltat
    Jupiter.p = Jupiter.p + Fnet1*deltat
    Earth.pos = Earth.pos+Earth.p/Earth.m*deltat
    Jupiter.pos = Jupiter.pos+Jupiter.p/Jupiter.m*deltat
    Earth.trail.append(pos=Earth.pos)
    Jupiter.trail.append(pos=Jupiter.pos)
    
##    parr.pos = Earth.pos
##    parr.axis = pscale*Earth.p
##    farr.pos = Earth.pos
##    farr.axis = fscale*Fnet
    
    t = t + deltat
    count += 1



    
    if count == 12 or count == 24 or count == 36 or count == 48 or count == 60:
        arrow(pos=Jupiter.pos, axis=(Star.pos - Jupiter.pos), shaftwidth=1e9)
        arrow(pos=Jupiter.pos, axis=(Jupiter.pos - Earth.pos), shaftwidth=1e9)
        arrow(pos=Earth.pos, axis=(Star.pos - Earth.pos), shaftwidth=1e9)


Earthlabel = label(
    pos=Earth.pos,
    text='Earth', 
    xoffset=20,
    yoffset=12, 
    height=10, 
    border=6,
    font='sans'
    )

Juplabel = label(
    pos=Jupiter.pos,
    text='Jupiter', 
    xoffset=20,
    yoffset=12, 
    height=10, 
    border=6,
    font='sans'
    )

North