let f t u =
    [| u.(1); (if u.(1) < 0. then
                 let u = abs_float u.(1) in
                 let u = u**(1. /. 3.) in
                 let u = -. u in
                 -.a *. u
               else -. a *. u.(1)**(1. /. 3.)) |] in