fun=inline('cos(pi/2.*cos(x))-0.707*sin(x)','x')
x1 =fzero(fun,[0.5 1]),x2 =fzero(fun,[2 2.5])
x1 =
0.8893
x2 =
2.2523