func ExampleIlluminated_moon() {
// Example 48.a, p. 347.
k := base.Illuminated(69.0756 * math.Pi / 180)
fmt.Printf("k = %.4f\n", k)
// Output:
// k = 0.6786
}
func ExampleIlluminated_venus() {
// Example 41.a, p. 284.
k := base.Illuminated(math.Acos(.29312))
fmt.Printf("%.3f\n", k)
// Output:
// 0.647
}
func ExampleIlluminated_moon() {
// Example 48.a, p. 347.
k := base.Illuminated(unit.AngleFromDeg(69.0756))
fmt.Printf("k = %.4f\n", k)
// Output:
// k = 0.6786
}
func ExamplePhaseAngleEq2() {
const p = math.Pi / 180
i := moonillum.PhaseAngleEq2(134.6885*p, 13.7684*p, 20.6579*p, 8.6964*p)
k := base.Illuminated(i)
fmt.Printf("k = %.4f\n", k)
// Output:
// k = 0.6775
}
func ExamplePhaseAngle3() {
i := moonillum.PhaseAngle3(julian.CalendarGregorianToJD(1992, 4, 12))
k := base.Illuminated(i)
fmt.Printf("i = %.2f\n", i*180/math.Pi)
fmt.Printf("k = %.4f\n", k)
// Output:
// i = 68.88
// k = 0.6801
}
func TestPhaseAngleEcl2(t *testing.T) {
j := julian.CalendarGregorianToJD(1992, 4, 12)
λ, β, _ := moonposition.Position(j)
λ0 := solar.ApparentLongitude(base.J2000Century(j))
i := moonillum.PhaseAngleEcl2(λ, β, λ0)
k := base.Illuminated(i)
ref := .6775
if math.Abs(k-ref) > 1e-4 {
t.Errorf("k = %.4f", k)
}
}
func ExamplePhaseAngleEq2() {
i := moonillum.PhaseAngleEq2(
unit.RAFromDeg(134.6885),
unit.AngleFromDeg(13.7684),
unit.RAFromDeg(20.6579),
unit.AngleFromDeg(8.6964))
k := base.Illuminated(i)
fmt.Printf("k = %.4f\n", k)
// Output:
// k = 0.6775
}
