Content-Type: text/html; charset=UTF-8 Double Star Calculator - Orbital Elements

Double Star Calculator - Orbit

- Revision: $Rev: 1216 $
- Report created: 2025-07-11 19:03:31

Orbital Elements for D19AB
WDS: D19AB WDS ID/Discoverer
P: 2858.000 Orbital period (years)
a: 2.95 Semi-major axis (arcseconds)
i: 100.0 Inclination (degrees)
Ω: 190.0 Longitude of ascending node (degrees)
T: 1903.000 Epoch of periastron passage (year)
e: 0.73 Eccentricity
ω: 162.0 Longitude of periastron (degrees)
Start-Year: 596.0 Year
End-Year: 3454.0 Year
Epoch: 2025.0 Besselian Epoch

Further data from the WDS Sixth Catalog of Orbits for D19AB
Position: 190820.20+551950.0 Epoch-2000 right ascension and declination
WDS-ID: 19083+5520 WDS designation
Discoverer Code: D19AB Discover designation and components
ADS number: 12104 Aitken Double Star catalog
HD number: 179142 HD catalog number
Hipparcos number: 94000 Hipparcos catalog number
V1: 7.67 Magnitude of the primary
V2: 9.51 Magnitude of the secondary
Equinox: 2000 Equinox
Date: 2009 Date of last observation
Orbit grade: 5 Orbit grade, ranging from 1 'definitive' to 5 'indeterminate'
Flags: n Flags

Double Star Orbit for D19AB
--> -->
 
 
ValueError
Python 3.12.3: /usr/bin/python3
Fri Jul 11 19:03:31 2025

A problem occurred in a Python script. Here is the sequence of function calls leading up to the error, in the order they occurred.

 /var/www/cgi-bin/ds_orbit.py in <module>
    891 
    892 # Smooth curve if necessary
=>  893 x_smooth, y_smooth = Q_get_smoothed_orbit(x_list, y_list)
    894 plt.plot(x_smooth, y_smooth, label=f'Orbit of {wds}', color='#528B8B')
    895 
x_smooth undefined, y_smooth undefined, Q_get_smoothed_orbit = <function Q_get_smoothed_orbit>, x_list = [1.0708088322562526, 1.0704040091513818, 1.0699981712300222, 1.0695913186187198, 1.0691834514429899, 1.0687745698273146, 1.0683646738951456, 1.0679537637689043, 1.0675418395699812, 1.067128901418739, 1.0667149494345112, 1.0662999837356024, 1.065884004439291, 1.0654670116618281, 1.0650490055184383, 1.0646299861233213, 1.0642099535896494, 1.0637889080295746, 1.0633668495542214, 1.0629437782736915, ...], y_list = [-4.8033494193843085, -4.803653930301478, -4.803953887024844, -4.804249288110718, -4.804540132109754, -4.8048264175669235, -4.805108143021531, -4.805385307007185, -4.805657908051797, -4.8059259446775755, -4.806189415401016, -4.806448318732886, -4.8067026531782275, -4.806952417236342, -4.8071976094007836, -4.807438228159342, -4.807674271994054, -4.8079057393811695, -4.808132628791162, -4.808354938688713, ...]
 /var/www/cgi-bin/ds_orbit.py in Q_get_smoothed_orbit(x_list=[1.0708088322562526, 1.0704040091513818, 1.0699981712300222, 1.0695913186187198, 1.0691834514429899, 1.0687745698273146, 1.0683646738951456, 1.0679537637689043, 1.0675418395699812, 1.067128901418739, 1.0667149494345112, 1.0662999837356024, 1.065884004439291, 1.0654670116618281, 1.0650490055184383, 1.0646299861233213, 1.0642099535896494, 1.0637889080295746, 1.0633668495542214, 1.0629437782736915, ...], y_list=[-4.8033494193843085, -4.803653930301478, -4.803953887024844, -4.804249288110718, -4.804540132109754, -4.8048264175669235, -4.805108143021531, -4.805385307007185, -4.805657908051797, -4.8059259446775755, -4.806189415401016, -4.806448318732886, -4.8067026531782275, -4.806952417236342, -4.8071976094007836, -4.807438228159342, -4.807674271994054, -4.8079057393811695, -4.808132628791162, -4.808354938688713, ...])
    272     # Nur wenn ausreichend Punkte vorhanden sind
    273     if len(x) >= 4:
=>  274         tck, u = splprep([x, y], s=0.0, per=True) # per=True für geschlossene Kurve (Orbit)
    275         u_fine = np.linspace(0, 1, 300)           # Feinere Parameterwerte
    276         return splev(u_fine, tck)
tck undefined, u undefined, splprep = <function splprep>, x = array([-0.19417028, -0.19330604, -0.19242457, ..., -0.19665573, -0.19584552, -0.1950169 ]), y = array([-0.00857006, -0.0202255 , -0.03187914, ..., 0.02639661, 0.01474207, 0.00308615]), s undefined, per undefined
 /usr/lib/python3/dist-packages/scipy/interpolate/_fitpack_py.py in splprep(x=[array([-0.19417028, -0.19330604, -0.19242457, ..., -0.19665573, -0.19584552, -0.1950169 ]), array([-0.00857006, -0.0202255 , -0.03187914, ..., 0.02639661, 0.01474207, 0.00308615])], w=None, u=None, ub=None, ue=None, k=3, task=0, s=0.0, t=None, full_output=0, nest=None, per=True, quiet=1)
    153     """
    154 
=>  155     res = _impl.splprep(x, w, u, ub, ue, k, task, s, t, full_output, nest, per,
=>  156                         quiet)
    157     return res
res undefined, global _impl = <module 'scipy.interpolate._fitpack_impl' from '...ist-packages/scipy/interpolate/_fitpack_impl.py'>, _impl.splprep = <function splprep>, x = [array([-0.19417028, -0.19330604, -0.19242457, ..., -0.19665573, -0.19584552, -0.1950169 ]), array([-0.00857006, -0.0202255 , -0.03187914, ..., 0.02639661, 0.01474207, 0.00308615])], w = None, u = None, ub = None, ue = None, k = 3, task = 0, s = 0.0, t = None, full_output = 0, nest = None, per = True, quiet = 1
 /usr/lib/python3/dist-packages/scipy/interpolate/_fitpack_impl.py in splprep(x=array([[-0.19417028, -0.19330604, -0.19242457, .... 0.02639661, 0.01474207, -0.00857006]]), w=array([1., 1., 1., ..., 1., 1., 1.]), u=array([0.00000000e+00, 9.93142417e-04, 1.9862416... 9.97021201e-01, 9.98013938e-01, 1.00000000e+00]), ub=0, ue=1, k=3, task=0, s=0.0, t=array([], dtype=float64), full_output=0, nest=2865, per=True, quiet=1)
    172     wrk = _parcur_cache['wrk']
    173     iwrk = _parcur_cache['iwrk']
=>  174     t, c, o = _fitpack._parcur(ravel(transpose(x)), w, u, ub, ue, k,
=>  175                                task, ipar, s, t, nest, wrk, iwrk, per)
    176     _parcur_cache['u'] = o['u']
t = array([], dtype=float64), c undefined, o undefined, global _fitpack = <module 'scipy.interpolate._fitpack' from '/usr/...polate/_fitpack.cpython-312-x86_64-linux-gnu.so'>, _fitpack._parcur = <built-in function _parcur>, global ravel = <function ravel>, global transpose = <function transpose>, x = array([[-0.19417028, -0.19330604, -0.19242457, .... 0.02639661, 0.01474207, -0.00857006]]), w = array([1., 1., 1., ..., 1., 1., 1.]), u = array([0.00000000e+00, 9.93142417e-04, 1.9862416... 9.97021201e-01, 9.98013938e-01, 1.00000000e+00]), ub = 0, ue = 1, k = 3, task = 0, ipar = False, s = 0.0, nest = 2865, wrk = array([], dtype=float64), iwrk = array([], dtype=int32), per = True

ValueError: Invalid inputs.
      add_note = <built-in method add_note of ValueError object>
      args = ('Invalid inputs.',)
      with_traceback = <built-in method with_traceback of ValueError object>