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

Double Star Calculator - Orbit

- Revision:   $Rev: 1216 $
- Report created:   2025-07-11 17:44:35

Orbital Elements for B1646AB
WDS: B1646AB WDS ID/Discoverer
P: 208.000 Orbital period (years)
a: 0.188 Semi-major axis (arcseconds)
i: 116.2 Inclination (degrees)
Ω: 15.9 Longitude of ascending node (degrees)
T: 2029.400 Epoch of periastron passage (year)
e: 0.6 Eccentricity
ω: 161.1 Longitude of periastron (degrees)
Start-Year: 1921.0 Year
End-Year: 2129.0 Year
Epoch: 2025.0 Besselian Epoch

Further data from the WDS Sixth Catalog of Orbits for B1646AB
Position: 090723.26-435718.3 Epoch-2000 right ascension and declination
WDS-ID: 09074-4357 WDS designation
Discoverer Code: B1646AB Discover designation and components
ADS number: . Aitken Double Star catalog
HD number: 78547 HD catalog number
Hipparcos number: . Hipparcos catalog number
V1: 9.5 Magnitude of the primary
V2: 9.9 Magnitude of the secondary
Equinox: Equinox
Date: 2019 Date of last observation
Orbit grade: 4 Orbit grade, ranging from 1 'definitive' to 5 'indeterminate'
Flags: Flags

Double Star Orbit for B1646AB
--> -->
 
 
ValueError		 Python 3.12.3: /usr/bin/python3
Fri Jul 11 17:44:35 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 = [-0.122986139985027, -0.1222026713153042, -0.12139191486624396, -0.1205540832244928, -0.11968938375944385, -0.11879801882587046, -0.11788018596180488, -0.11693607808241879, -0.11596588367073822, -0.11496978696610301, -0.11394796815135502, -0.11290060349613705, -0.11182786564323577, -0.11072992367026688, -0.10960694333418629, -0.10845908723816983, -0.10728651501263967, -0.106089383497493, -0.10486784692611748, -0.10362205711180522, ...], y_list = [0.2560265826234579, 0.2574616956937705, 0.25883931789189096, 0.260159208808716, 0.261421118248795, 0.262624785979084, 0.2637699414713765, 0.2648563036371288, 0.2658835805535795, 0.26685146918022307, 0.2677596550648831, 0.26860781207249995, 0.2693956019840026, 0.27012267426957665, 0.27078866569510546, 0.27139319998073697, 0.27193588743358416, 0.27241632456516823, 0.2728340936927159, 0.2731887625233801, ...]
 /var/www/cgi-bin/ds_orbit.py in Q_get_smoothed_orbit(x_list=[-0.122986139985027, -0.1222026713153042, -0.12139191486624396, -0.1205540832244928, -0.11968938375944385, -0.11879801882587046, -0.11788018596180488, -0.11693607808241879, -0.11596588367073822, -0.11496978696610301, -0.11394796815135502, -0.11290060349613705, -0.11182786564323577, -0.11072992367026688, -0.10960694333418629, -0.10845908723816983, -0.10728651501263967, -0.106089383497493, -0.10486784692611748, -0.10362205711180522, ...], y_list=[0.2560265826234579, 0.2574616956937705, 0.25883931789189096, 0.260159208808716, 0.261421118248795, 0.262624785979084, 0.2637699414713765, 0.2648563036371288, 0.2658835805535795, 0.26685146918022307, 0.2677596550648831, 0.26860781207249995, 0.2693956019840026, 0.27012267426957665, 0.27078866569510546, 0.27139319998073697, 0.27193588743358416, 0.27241632456516823, 0.2728340936927159, 0.2731887625233801, ...])
    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.07226715, -0.06924986, -0.06612336, -0...08330992, -0.08069683, -0.07798741, -0.07517857]), y = array([-0.00374145, -0.00880208, -0.01384873, -0...01648842, 0.01144562, 0.00638923, 0.00132475]), s undefined, per undefined
 /usr/lib/python3/dist-packages/scipy/interpolate/_fitpack_py.py in splprep(x=[array([-0.07226715, -0.06924986, -0.06612336, -0...08330992, -0.08069683, -0.07798741, -0.07517857]), array([-0.00374145, -0.00880208, -0.01384873, -0...01648842, 0.01144562, 0.00638923, 0.00132475])], 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.07226715, -0.06924986, -0.06612336, -0...08330992, -0.08069683, -0.07798741, -0.07517857]), array([-0.00374145, -0.00880208, -0.01384873, -0...01648842, 0.01144562, 0.00638923, 0.00132475])], 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([[-7.22671464e-02, -6.92498579e-02, -6.612...e-02, 6.38922533e-03, -3.74145323e-03]]), w=array([1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1..., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1.]), u=array([0. , 0.00704958, 0.01415272, 0.021... 0.97242049, 0.97921611, 0.98607988, 1. ]), ub=0, ue=1, k=3, task=0, s=0.0, t=array([], dtype=float64), full_output=0, nest=215, 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([[-7.22671464e-02, -6.92498579e-02, -6.612...e-02, 6.38922533e-03, -3.74145323e-03]]), w = array([1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1..., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1.]), u = array([0. , 0.00704958, 0.01415272, 0.021... 0.97242049, 0.97921611, 0.98607988, 1. ]), ub = 0, ue = 1, k = 3, task = 0, ipar = False, s = 0.0, nest = 215, 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>