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:46:18

Orbital Elements for WSI87AD

WDS:	WSI87AD	WDS ID/Discoverer
P:	300.000	Orbital period (years)
a:	0.28	Semi-major axis (arcseconds)
i:	90.0	Inclination (degrees)
Ω:	271.5	Longitude of ascending node (degrees)
T:	1969.000	Epoch of periastron passage (year)
e:	0.8	Eccentricity
ω:	242.0	Longitude of periastron (degrees)
Start-Year:	1875.0	Year
End-Year:	2175.0	Year
Epoch:	2025.0	Besselian Epoch

Further data from the WDS Sixth Catalog of Orbits for WSI87AD

Position:	172445.86-591249.1	Epoch-2000 right ascension and declination
WDS-ID:	17248-5913	WDS designation
Discoverer Code:	WSI87AD	Discover designation and components
ADS number:	.	Aitken Double Star catalog
HD number:	157081	HD catalog number
Hipparcos number:	85216	Hipparcos catalog number
V1:	7.56	Magnitude of the primary
V2:	8.1	Magnitude of the secondary
Equinox:	2000	Equinox
Date:	2022	Date of last observation
Orbit grade:	4	Orbit grade, ranging from 1 'definitive' to 5 'indeterminate'
Flags:	n	Flags

Double Star Orbit for WSI87AD

--> -->

ValueError

Python 3.12.3: /usr/bin/python3
Fri Jul 11 18:46: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 = [0.11530389425045344, 0.11274807192614553, 0.11018058906187748, 0.10760157009477433, 0.10501114218256372, 0.10240943525914871, 0.09979658274943953, 0.09717272105925781, 0.0945379902607053, 0.09189253422426114, 0.08923650087447865, 0.08657004246151999, 0.08389331584975708, 0.0812064828247998, 0.07850971042041364, 0.07580317126693015, 0.0730870439629145, 0.07036151347199396, 0.06762677154695212, 0.06488301718338732, ...], y_list = [0.00301933873146422, 0.0029524121685351157, 0.0028851802636211616, 0.0028176462752231533, 0.002749813533083373, 0.0026816854396408177, 0.0026132654887497195, 0.0025445572523222267, 0.0024755643982754096, 0.00240629069399035, 0.0023367400130078013, 0.002266916342138567, 0.0021968237890207885, 0.0021264665901597557, 0.002055849119488552, 0.0019849758974914995, 0.001913851600936619, 0.0018424810732669457, 0.0017708693357058336, 0.0016990215991364108, ...]

/var/www/cgi-bin/ds_orbit.py in Q_get_smoothed_orbit(x_list=[0.11530389425045344, 0.11274807192614553, 0.11018058906187748, 0.10760157009477433, 0.10501114218256372, 0.10240943525914871, 0.09979658274943953, 0.09717272105925781, 0.0945379902607053, 0.09189253422426114, 0.08923650087447865, 0.08657004246151999, 0.08389331584975708, 0.0812064828247998, 0.07850971042041364, 0.07580317126693015, 0.0730870439629145, 0.07036151347199396, 0.06762677154695212, 0.06488301718338732, ...], y_list=[0.00301933873146422, 0.0029524121685351157, 0.0028851802636211616, 0.0028176462752231533, 0.002749813533083373, 0.0026816854396408177, 0.0026132654887497195, 0.0025445572523222267, 0.0024755643982754096, 0.00240629069399035, 0.0023367400130078013, 0.002266916342138567, 0.0021968237890207885, 0.0021264665901597557, 0.002055849119488552, 0.0019849758974914995, 0.001913851600936619, 0.0018424810732669457, 0.0017708693357058336, 0.0016990215991364108, ...])

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([-9.16567157e-05, -2.96055636e-03, -5.8287...55e-03, 5.64555579e-03, 2.77726843e-03]), y = array([-2.40011557e-06, -7.75248967e-05, -1.5263...19e-04, 1.47834081e-04, 7.27253334e-05]), s undefined, per undefined

/usr/lib/python3/dist-packages/scipy/interpolate/_fitpack_py.py in splprep(x=[array([-9.16567157e-05, -2.96055636e-03, -5.8287...55e-03, 5.64555579e-03, 2.77726843e-03]), array([-2.40011557e-06, -7.75248967e-05, -1.5263...19e-04, 1.47834081e-04, 7.27253334e-05])], 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([-9.16567157e-05, -2.96055636e-03, -5.8287...55e-03, 5.64555579e-03, 2.77726843e-03]), array([-2.40011557e-06, -7.75248967e-05, -1.5263...19e-04, 1.47834081e-04, 7.27253334e-05])], 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([[-9.16567157e-05, -2.96055636e-03, -5.828...e-04, 1.47834081e-04, -2.40011557e-06]]), 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.00167811, 0.00335579, 0.005...9118, 0.99496712, 0.99664413, 1. ]), ub=0, ue=1, k=3, task=0, s=0.0, t=array([], dtype=float64), full_output=0, nest=307, 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([[-9.16567157e-05, -2.96055636e-03, -5.828...e-04, 1.47834081e-04, -2.40011557e-06]]), 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.00167811, 0.00335579, 0.005...9118, 0.99496712, 0.99664413, 1. ]), ub = 0, ue = 1, k = 3, task = 0, ipar = False, s = 0.0, nest = 307, 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>