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
}
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
}