In [11]:
from numpy import *
from sympy import *

In [15]:
x,y = var('x y')

In [16]:
f = csc(x)-cot(x)
f.diff()


Out[16]:
cot(x)**2 + Derivative(csc(x), x) + 1

In [ ]:
x =