func ExampleApproxPosition() {
// Example 21.a, p. 132.
eq := &coord.Equatorial{
base.NewRA(10, 8, 22.3).Rad(),
base.NewAngle(false, 11, 58, 2).Rad(),
}
epochFrom := 2000.0
epochTo := 1978.0
mα := base.NewHourAngle(true, 0, 0, 0.0169)
mδ := base.NewAngle(false, 0, 0, 0.006)
precess.ApproxPosition(eq, eq, epochFrom, epochTo, mα, mδ)
fmt.Printf("%0.1d\n", base.NewFmtRA(eq.RA))
fmt.Printf("%+0d\n", base.NewFmtAngle(eq.Dec))
// Output:
// 10ʰ07ᵐ12ˢ.1
// +12°04′32″
}
func ExampleApproxPosition() {
// Example 21.a, p. 132.
eq := &coord.Equatorial{
unit.NewRA(10, 8, 22.3),
unit.NewAngle(' ', 11, 58, 2),
}
epochFrom := 2000.0
epochTo := 1978.0
mα := unit.HourAngleFromSec(-0.0169)
mδ := unit.AngleFromSec(0.006)
precess.ApproxPosition(eq, eq, epochFrom, epochTo, mα, mδ)
fmt.Printf("%0.1d\n", sexa.FmtRA(eq.RA))
fmt.Printf("%+0d\n", sexa.FmtAngle(eq.Dec))
// Output:
// 10ʰ07ᵐ12ˢ.1
// +12°04′32″
}
