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:41:03

Orbital Elements for BU79AB

WDS:	BU79AB	WDS ID/Discoverer
P:	302.000	Orbital period (years)
a:	1.32	Semi-major axis (arcseconds)
i:	123.6	Inclination (degrees)
Ω:	172.8	Longitude of ascending node (degrees)
T:	2112.000	Epoch of periastron passage (year)
e:	0.451	Eccentricity
ω:	303.0	Longitude of periastron (degrees)
Start-Year:	1874.0	Year
End-Year:	2176.0	Year
Epoch:	2025.0	Besselian Epoch

Further data from the WDS Sixth Catalog of Orbits for BU79AB

Position:	231737.61-013118.2	Epoch-2000 right ascension and declination
WDS-ID:	23176-0131	WDS designation
Discoverer Code:	BU79AB	Discover designation and components
ADS number:	16649	Aitken Double Star catalog
HD number:	219657	HD catalog number
Hipparcos number:	115012	Hipparcos catalog number
V1:	8.23	Magnitude of the primary
V2:	9.59	Magnitude of the secondary
Equinox:	2000	Equinox
Date:	2003	Date of last observation
Orbit grade:	4	Orbit grade, ranging from 1 'definitive' to 5 'indeterminate'
Flags:	n	Flags

Double Star Orbit for BU79AB

--> -->

ValueError

Python 3.12.3: /usr/bin/python3
Fri Jul 11 19:41:04 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>

956

957 # Smooth curve if necessary

=> 958 xa_smooth, ya_smooth = Q_get_smoothed_orbit(xa_list, ya_list)

959 plt.plot(xa_smooth, ya_smooth, label=f'Orbit of {wds}', color='#528B8B')

960

xa_smooth undefined, ya_smooth undefined, Q_get_smoothed_orbit = <function Q_get_smoothed_orbit>, xa_list = [0.8040028210883816, 0.8117156762596737, 0.8191638918154774, 0.8263509847880628, 0.8332804163537667, 0.8399555911502791, 0.8463798568344854, 0.8525565038511671, 0.8584887653860961, 0.8641798174799876, 0.8696327792823573, 0.8748507134266532, 0.8798366265101075, 0.8845934696635885, 0.889124139198399, 0.8934314773184229, 0.8975182728873508, 0.9013872622418798, 0.905041130042832, 0.9084825101570712, ...], ya_list = [0.401428445601899, 0.38301083657672436, 0.3644683291877286, 0.3458096153558042, 0.327043051300054, 0.3081766739841112, 0.28921821670504766, 0.2701751238696561, 0.2510545650011312, 0.23186344801726982, 0.21260843181952502, 0.19329593823029204, 0.17393216331410022, 0.15452308811647758, 0.1350744888525972, 0.11559194657611997, 0.09608085635694709, 0.07654643599518395, 0.05699373429696997, 0.037427638936523185, ...]

/var/www/cgi-bin/ds_orbit.py in Q_get_smoothed_orbit(x_list=[0.8040028210883816, 0.8117156762596737, 0.8191638918154774, 0.8263509847880628, 0.8332804163537667, 0.8399555911502791, 0.8463798568344854, 0.8525565038511671, 0.8584887653860961, 0.8641798174799876, 0.8696327792823573, 0.8748507134266532, 0.8798366265101075, 0.8845934696635885, 0.889124139198399, 0.8934314773184229, 0.8975182728873508, 0.9013872622418798, 0.905041130042832, 0.9084825101570712, ...], y_list=[0.401428445601899, 0.38301083657672436, 0.3644683291877286, 0.3458096153558042, 0.327043051300054, 0.3081766739841112, 0.28921821670504766, 0.2701751238696561, 0.2510545650011312, 0.23186344801726982, 0.21260843181952502, 0.19329593823029204, 0.17393216331410022, 0.15452308811647758, 0.1350744888525972, 0.11559194657611997, 0.09608085635694709, 0.07654643599518395, 0.05699373429696997, 0.037427638936523185, ...])

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([-4.30555624e-01, -4.38327106e-01, -4.4518....01648496e-01, -4.12224351e-01, -4.21858131e-01]), y = array([-0.00201353, -0.04362089, -0.08513802, -0...1, 0.12245754, 0.08112054, 0.03959769]), s undefined, per undefined

/usr/lib/python3/dist-packages/scipy/interpolate/_fitpack_py.py in splprep(x=[array([-4.30555624e-01, -4.38327106e-01, -4.4518....01648496e-01, -4.12224351e-01, -4.21858131e-01]), array([-0.00201353, -0.04362089, -0.08513802, -0...1, 0.12245754, 0.08112054, 0.03959769])], 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([-4.30555624e-01, -4.38327106e-01, -4.4518....01648496e-01, -4.12224351e-01, -4.21858131e-01]), array([-0.00201353, -0.04362089, -0.08513802, -0...1, 0.12245754, 0.08112054, 0.03959769])], 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([[-4.30555624e-01, -4.38327106e-01, -4.451...22457542e-01, 8.11205424e-02, -2.01352573e-03]]), 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.00686287, 0.0136857 , 0.020...6589, 0.97927864, 0.98619688, 1. ]), ub=0, ue=1, k=3, task=0, s=0.0, t=array([], dtype=float64), full_output=0, nest=309, 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([[-4.30555624e-01, -4.38327106e-01, -4.451...22457542e-01, 8.11205424e-02, -2.01352573e-03]]), 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.00686287, 0.0136857 , 0.020...6589, 0.97927864, 0.98619688, 1. ]), ub = 0, ue = 1, k = 3, task = 0, ipar = False, s = 0.0, nest = 309, 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>