x=var('x')   #define a variable x
y=function('y',x)   #define y to be a function of that variable
t= diff(y,x)-(2*y)*cos(2*x)/sin(2*x)
desolve(t,[y,x],ics=(pi/4,3.8))