Exemple #1
0
func ExamplePhysical() {
	// Example 42.a, p. 291
	e, err := pp.LoadPlanet(pp.Earth, "")
	if err != nil {
		fmt.Println(err)
		return
	}
	m, err := pp.LoadPlanet(pp.Mars, "")
	if err != nil {
		fmt.Println(err)
		return
	}
	DE, DS, ω, P, Q, d, k, q := mars.Physical(2448935.500683, e, m)
	fmt.Printf("DE = %+.2f\n", DE*180/math.Pi)
	fmt.Printf("DS = %+.2f\n", DS*180/math.Pi)
	fmt.Printf("ω = %.2f\n", ω*180/math.Pi)
	fmt.Printf("P = %.2f\n", P*180/math.Pi)
	fmt.Printf("Q = %.2f\n", Q*180/math.Pi)
	fmt.Printf("d = %.2f\n", d*180/math.Pi*60*60) // display as arc sec
	fmt.Printf("k = %.4f\n", k)
	fmt.Printf("q = %.2f\n", q*180/math.Pi*60*60) // display as arc sec
	// Output:
	// DE = +12.44
	// DS = -2.76
	// ω = 111.55
	// P = 347.64
	// Q = 279.91
	// d = 10.75
	// k = 0.9012
	// q = 1.06
}
Exemple #2
0
func ExamplePhysical() {
	// Example 42.a, p. 291
	e, err := pp.LoadPlanet(pp.Earth)
	if err != nil {
		fmt.Println(err)
		return
	}
	m, err := pp.LoadPlanet(pp.Mars)
	if err != nil {
		fmt.Println(err)
		return
	}
	DE, DS, ω, P, Q, d, q, k := mars.Physical(2448935.500683, e, m)
	fmt.Printf("DE = %+.2f\n", DE.Deg())
	fmt.Printf("DS = %+.2f\n", DS.Deg())
	fmt.Printf("ω = %.2f\n", ω.Deg())
	fmt.Printf("P = %.2f\n", P.Deg())
	fmt.Printf("Q = %.2f\n", Q.Deg())
	fmt.Printf("d = %.2f\n", d.Sec()) // display as arc sec
	fmt.Printf("k = %.4f\n", k)
	fmt.Printf("q = %.2f\n", q.Sec()) // display as arc sec
	// Output:
	// DE = +12.44
	// DS = -2.76
	// ω = 111.55
	// P = 347.64
	// Q = 279.91
	// d = 10.75
	// k = 0.9012
	// q = 1.06
}