[(-(y_0 - y_c) / c, y_0, 0, 0), ((2 * c * dy_0 * dy_c + dy_0 * y_0 - dy_0 *
    y_c + dy_c * y_0 - dy_c * y_c) / (c * dy_0 + c * dy_c + 2 * y_0 - 2 *
    y_c), -(2 * c ** 3 * dy_0 * dy_c ** 2 - c ** 2 * dy_0 ** 2 * y_0 + 2 * 
    c ** 2 * dy_0 * dy_c * y_0 - 4 * c ** 2 * dy_0 * dy_c * y_c + c ** 2 * 
    dy_c ** 2 * y_0 - 2 * c ** 2 * dy_c ** 2 * y_c - 2 * c * dy_0 * y_0 ** 
    2 + 2 * c * dy_0 * y_c ** 2 - 4 * c * dy_c * y_0 * y_c + 4 * c * dy_c *
    y_c ** 2 - 2 * y_0 ** 3 + 2 * y_0 ** 2 * y_c + 2 * y_0 * y_c ** 2 - 2 *
    y_c ** 3) / (c * dy_0 + c * dy_c + 2 * y_0 - 2 * y_c) ** 2, -4 * (dy_0 -
    dy_c) * (c * dy_0 + y_0 - y_c) ** 2 * (c * dy_c + y_0 - y_c) ** 2 / (c *
    dy_0 + c * dy_c + 2 * y_0 - 2 * y_c) ** 3, -4 * (c * dy_0 + y_0 - y_c) **
    2 * (c * dy_c + y_0 - y_c) ** 4 / (c * dy_0 + c * dy_c + 2 * y_0 - 2 *
    y_c) ** 4)]
How do i solve a non linear equation in Sympy?