x=solve('x^2+2*x*sin(x-0.571)+2*cos(x-0.571)-(pi-0.571)^2+2')
x =
1.2317050037279494364890533718536481