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

Double Star Calculator - Orbit

- Revision:   $Rev: 1216 $
- Report created:   2025-07-11 18:43:18

Orbital Elements for BU823AB
WDS: BU823AB WDS ID/Discoverer
P: 900.000 Orbital period (years)
a: 2.264 Semi-major axis (arcseconds)
i: 40.0 Inclination (degrees)
Ω: 97.7 Longitude of ascending node (degrees)
T: 1963.040 Epoch of periastron passage (year)
e: 0.596 Eccentricity
ω: 356.9 Longitude of periastron (degrees)
Start-Year: 1575.0 Year
End-Year: 2475.0 Year
Epoch: 2025.0 Besselian Epoch

Further data from the WDS Sixth Catalog of Orbits for BU823AB
Position: 170636.18+003907.9 Epoch-2000 right ascension and declination
WDS-ID: 17066+0039 WDS designation
Discoverer Code: BU823AB Discover designation and components
ADS number: 10341 Aitken Double Star catalog
HD number: 154621 HD catalog number
Hipparcos number: 83716 Hipparcos catalog number
V1: 8.73 Magnitude of the primary
V2: 9.81 Magnitude of the secondary
Equinox: Equinox
Date: Date of last observation
Orbit grade: 4 Orbit grade, ranging from 1 'definitive' to 5 'indeterminate'
Flags: Flags

Double Star Orbit for BU823AB
--> -->
 
 
ValueError		 Python 3.12.3: /usr/bin/python3
Fri Jul 11 18:43:18 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 = [3.4886777332636707, 3.485589019627618, 3.482456685101099, 3.47928070830376, 3.476061067430713, 3.472797740251765, 3.4694907041106253, 3.4661399359240908, 3.4627454121812273, 3.459307108942532, 3.455825001839082, 3.4522990660716646, 3.448729276409892, 3.445115607191312, 3.4414580323204857, 3.43775652526807, 3.434011059069872, 3.4302216063258895, 3.426388139199349, 3.4225106294157124, ...], y_list = [0.7014978524958763, 0.7071862002488698, 0.7128656978032689, 0.7185362616809409, 0.7241978080949824, 0.7298502529446268, 0.7354935118101319, 0.7411274999476262, 0.7467521322839009, 0.7523673234111572, 0.7579729875817218, 0.763569038702695, 0.7691553903305929, 0.7747319556658923, 0.7802986475475653, 0.7858553784475433, 0.7914020604651381, 0.7969386053214292, 0.8024649243535608, 0.8079809285090239, ...]
 /var/www/cgi-bin/ds_orbit.py in Q_get_smoothed_orbit(x_list=[3.4886777332636707, 3.485589019627618, 3.482456685101099, 3.47928070830376, 3.476061067430713, 3.472797740251765, 3.4694907041106253, 3.4661399359240908, 3.4627454121812273, 3.459307108942532, 3.455825001839082, 3.4522990660716646, 3.448729276409892, 3.445115607191312, 3.4414580323204857, 3.43775652526807, 3.434011059069872, 3.4302216063258895, 3.426388139199349, 3.4225106294157124, ...], y_list=[0.7014978524958763, 0.7071862002488698, 0.7128656978032689, 0.7185362616809409, 0.7241978080949824, 0.7298502529446268, 0.7354935118101319, 0.7411274999476262, 0.7467521322839009, 0.7523673234111572, 0.7579729875817218, 0.763569038702695, 0.7691553903305929, 0.7747319556658923, 0.7802986475475653, 0.7858553784475433, 0.7914020604651381, 0.7969386053214292, 0.8024649243535608, 0.8079809285090239, ...])
    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.91172013, -0.91189254, -0.91139288, -0...2, -0.90720097, -0.90936998, -0.91087804]), y = array([-0.01153099, -0.03567313, -0.05978897, -0...4, 0.06087574, 0.03676103, 0.01261962]), s undefined, per undefined
 /usr/lib/python3/dist-packages/scipy/interpolate/_fitpack_py.py in splprep(x=[array([-0.91172013, -0.91189254, -0.91139288, -0...2, -0.90720097, -0.90936998, -0.91087804]), array([-0.01153099, -0.03567313, -0.05978897, -0...4, 0.06087574, 0.03676103, 0.01261962])], 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.91172013, -0.91189254, -0.91139288, -0...2, -0.90720097, -0.90936998, -0.91087804]), array([-0.01153099, -0.03567313, -0.05978897, -0...4, 0.06087574, 0.03676103, 0.01261962])], 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.91172013, -0.91189254, -0.91139288, .... 0.06087574, 0.03676103, -0.01153099]]), w=array([1., 1., 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.00207236, 0.00414286, 0.006...9114, 0.99377149, 0.99584981, 1. ]), ub=0, ue=1, k=3, task=0, s=0.0, t=array([], dtype=float64), full_output=0, nest=907, 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.91172013, -0.91189254, -0.91139288, .... 0.06087574, 0.03676103, -0.01153099]]), w = array([1., 1., 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.00207236, 0.00414286, 0.006...9114, 0.99377149, 0.99584981, 1. ]), ub = 0, ue = 1, k = 3, task = 0, ipar = False, s = 0.0, nest = 907, 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>