x=var('x')
     y=function('y',x)
     desolve(diff(y,x)+(2*x*e^x^2*cosh(y))/(e^x^2),y)