2211.sagews

Ue=6730000*pi #earth girt

def durchmesser(u):
    return u/pi

def planetenoffset(ubase,uadd):
    uadd+=ubase
    return (durchmesser(uadd)-durchmesser(ubase))/2.0
    
plot(planetenoffset(Ue,x),(0,5))

import numpy
m=numpy.array([0],dtype=float)
m[0]=planetenoffset(Ue,1)
print(m[0])

0.159154943191