I did it, I used plane normal as basis and created a rotation matrix and a translation matrix using a point on the plane, transformed spline using inverse of it, and then just calculated roots of spline’s cubic polynomial for a single dimension.
took me a while to figure out, neither google nor mathematica did not help much.