Matlab - Ayuda matlab problema biseccion

function vm=mrastreo2(f,a,b,n)

% asi ya se puede evaluar al funcion

prt=linspace(a,b,(2^n)+1);

vm=[];

flag=0;

for i=1:n

if f(prt(i))*f(prt(i+1))<=0 % hay un cambio de signo

      vm=[vm;[prt(i) prt(i+1)]];

    flag=flag+1;

end

end

if flag==0

    display('no se encontraron cambios de digno')

    v=[]

else if n>10

    display('error')

end

end

>> syms x; f=x/2-cos(x)+1/2; a=-pi, b=pi; n=3,  vm=mrastreo2(f,a,b,n)

a =

   -3.1416

n =

     3

vm =

   -3.1416   -2.3562

   -2.3562   -1.5708