func TestNode(t *testing.T) { j := julian.CalendarGregorianToJD(2065, 6, 24) var e pe.Elements pe.Mean(pe.Mercury, j, &e) if Ω := pe.Node(pe.Mercury, j); Ω != e.Node { t.Fatal(Ω, "!=", e.Node) } }
func TestInc(t *testing.T) { j := julian.CalendarGregorianToJD(2065, 6, 24) var e pe.Elements pe.Mean(pe.Mercury, j, &e) if i := pe.Inc(pe.Mercury, j); i != e.Inc { t.Fatal(i, "!=", e.Inc) } }
func ExampleEllipticAscending_venus() { // Example 39.c, p. 278 var k planetelements.Elements planetelements.Mean(planetelements.Venus, julian.CalendarGregorianToJD(1979, 1, 1), &k) t, _ := node.EllipticAscending(k.Axis, k.Ecc, k.Peri-k.Node, perihelion.Perihelion(perihelion.Venus, 1979)) y, m, d := julian.JDToCalendar(t) fmt.Printf("%d %s %.3f\n", y, time.Month(m), d) // Output: // 1978 November 27.409 }
func ExampleMean() { // Example 31.a, p. 211 j := julian.CalendarGregorianToJD(2065, 6, 24) var e pe.Elements pe.Mean(pe.Mercury, j, &e) fmt.Printf("L: %.6f\n", e.Lon.Deg()) fmt.Printf("a: %.9f\n", e.Axis) fmt.Printf("e: %.8f\n", e.Ecc) fmt.Printf("i: %.6f\n", e.Inc.Deg()) fmt.Printf("Ω: %.6f\n", e.Node.Deg()) fmt.Printf("ϖ: %.6f\n", e.Peri.Deg()) // Output: // L: 203.494701 // a: 0.387098310 // e: 0.20564510 // i: 7.006171 // Ω: 49.107650 // ϖ: 78.475382 }