x=0.7 f1=sym('pi/2*(cos(a))^2+x*sin(2*a)=0.2'); c=solve(f1,a