\frac{e \cdot \sin v}{1 + e \cdot \cos v}
\sin v \cdot \left(e + \cos v \cdot \left({e}^{3} - e \cdot e\right)\right)
(FPCore (e v) :precision binary64 (/ (* e (sin v)) (+ 1.0 (* e (cos v)))))
(FPCore (e v) :precision binary64 (* (sin v) (+ e (* (cos v) (- (pow e 3.0) (* e e))))))
double code(double e, double v) {
return (e * sin(v)) / (1.0 + (e * cos(v)));
}
double code(double e, double v) {
return sin(v) * (e + (cos(v) * (pow(e, 3.0) - (e * e))));
}



Bits error versus e



Bits error versus v
Results
Initial program 0.1
Simplified0.1
Taylor expanded in e around 0 0.5
Simplified0.5
Taylor expanded in v around 0 0.6
Final simplification0.6
herbie shell --seed 2021313
(FPCore (e v)
:name "Trigonometry A"
:precision binary64
:pre (and (<= 0.0 e) (<= e 1.0))
(/ (* e (sin v)) (+ 1.0 (* e (cos v)))))