% Este programa ilustra el Método de Falsa Posición.
clear
i=0;
p0=1;
p1=2;
q0=demof_(p0);
q1=demof_(p1);
p=p1-(q1*(p1-p0))/(q1-q0);
 f=demof_(p);
if f==0
   fprintf('La raíz es %12.5f.\n',p);
else
   while abs(p1-p0)>.000001
      p=p1-(q1*(p1-p0))/(q1-q0);
      pn=demof_(p)*demof_(p0);
        if pn > 0
        p0=p;
        else
        p1=p;
     end
     i=i+1;
     fprintf('#%i La función evaluada en %12.10f ',i,p);
     fprintf('es: %12.10f\n',demof_(p));
end
end   




