reset()
x = var('x')
y = function('y', x)
DE = desolve(diff(y,x)-(1/sin(y)), y)
DE.simplify()