pub fn astroid(x: f64, y: f64) -> f64
Solve k^4+2k^3-(x^2+y^2-1)k^2-2y^2k-y^2 = 0 for positive root k This solution is adapted from Geocentric::Reverse