let y_0 x = (1. -. (x ** 2.)) *. beta