Orbital elements: Difference between revisions
m
→{{header|Sidef}}: use the Vector class added in Sidef 3.60
m (→version 2: changed whitespace.) |
m (→{{header|Sidef}}: use the Vector class added in Sidef 3.60) |
||
Line 1,173:
true_anomaly
) {
var (i, j, k) = (
)
Line 1,194 ⟶ 1,195:
rotate(\i, \j, argument_of_periapsis)
var l = (eccentricity == 1
var (c, s) = with(true_anomaly) { (.cos, .sin) }
Line 1,204 ⟶ 1,205:
var speed = muladd(i, rprime*c - r*s, j, rprime*s + r*c)
speed /= speed
speed *= sqrt(2/r - 1/semimajor_axis)
struct Result { position, speed }
Result(
}
Line 1,224 ⟶ 1,225:
<pre>
Arguments: [1, 1/10, 0, 355/678, 0, 0]:
Position :
Speed :
Arguments: [1, 1/10, 0.174532925199432957692369076848861271344287188854, 0.523598775598298873077107230546583814032861566563, 0.785398163397448309615660845819875721049292349844, 0]:
Position :
Speed :
</pre>
| |||