我想用sympy解下面的表达式。我已经定义了所有的符号和函数。c、 x，r是函数，f是我要解的p的方程。我将有两个可能的p，我想画如下。p也在0和1之间，我没有把它放进去，也不知道如何定义它。x是c函数对k的导数，我把它的逆函数放在f中。

代码：

import sympy as sym
sym.init_printing()

p,q,a,c0,c1,z,k = sym.symbols('p,q,a,c0,c1,z,k')

c = sym.Eq(c0+0.5*c1*(1/k**2))

x = sym.diff(c,k)

R = sym.Eq(z*(a+x**(-1)))

f = sym.Eq((1-p+p*q)*(a+x**(-1))-c-(1-p)*R)

sym.solve([f],(p))