func TestAnomalyDistance2(t *testing.T) {
	var e nearparabolic.Elements
	for _, d := range tdat2 {
		e.TimeP = base.J2000 + rand.Float64()*base.JulianCentury
		e.PDis = d.q
		e.Ecc = d.e
		ν, _, err := e.AnomalyDistance(e.TimeP + d.t)
		if (err == nil) != d.c {
			t.Errorf("%#v", d)
			continue
		}
		if math.Abs(ν.Deg()-d.ν) > math.Pow(10, float64(-d.p)) {
			t.Errorf("got ν = %.*f expected %.*f",
				d.p+1, ν.Deg(), d.p+1, d.ν)
		}
	}
}
func TestAnomalyDistance(t *testing.T) {
	var e nearparabolic.Elements
	for _, d := range tdat {
		e.TimeP = base.J2000 + rand.Float64()*base.JulianCentury
		e.PDis = d.q
		e.Ecc = d.e
		ν, r, err := e.AnomalyDistance(e.TimeP + d.t)
		if err != nil {
			t.Error(err)
			continue
		}
		if math.Abs(ν.Deg()-d.ν) > 1e-5 {
			t.Errorf("got ν = %.6f expected %.6f", ν.Deg(), d.ν)
		}
		if math.Abs(r-d.r) > 1e-6 {
			t.Errorf("got r = %.7f expected %.7f", r, d.r)
		}
	}
}