x=var('x')
y=function('y',x)
deq1=diff(y,x)-(2*(y-1)*tanh(2*x))==0
r=desolve(deq1,[y,x],ics=(0,4))
show(r.simplify_full())