Matlab - Index exceeds matrix dimensions.

clc

clear all

x=-5:0.01:5;

fdx=sin(x);

plot(x,fdx,'color','r');

line([-5,5],[0,0],'color', 'k');

%genera la grafica de la funcion para asi encontrar Xi y Xu

xi=input('Valor xi:')

xu=input('Valor xu:')

%Evaluar la funcion en los valores

for i=1:30

fdxi=sin(xi)

fdxu=sin(xu)

if (fdxi*fdxu)<0

    xr=(xi+xu)/2;

    fdxr=sin(xr);

    if(fdxi*fdxr)<0

        xu=xr;

    else if(fdxi*fdxr)>0

            xi=xr;

        end

    end

end

if(fdxi*fdxr)==0

    break

end

hold on

plot(xr,0,'o')

end

plot(xr,0,'*','color','k');

xr=0; xr=NaN; n=10;

for i=1:n

    xr(i)=(xi+xu)/2

    e(i)=(abs(xr(i+1)-xr(i))/xr(i+1)*100)

end

xr(1) = (xi+xu)/2;

for i=1:n

    xr(i+1)=(xi+xu)/2

    e(i)=(abs(xr(i+1)-xr(i))/xr(i+1)*100)

end

clc

clear

close all

x=-5:0.01:5;

fdx=sin(x);

plot(x,fdx,'color','r');

 hold on

line([-5,5],[0,0],'color', 'k');

%genera la grafica de la funcion para asi encontrar xi y xu

%xi=input('Valor xi:');

%xu=input('Valor xu:');

xi=-1.4; xu=1.3;

%Evaluar la funcion en los valores

for i=1:30

    fdxi=sin(xi);

    fdxu=sin(xu);

    if fdxi*fdxu <0

             xr=(xi+xu)/2;

             fdxr=sin(xr);

            if fdxi*fdxr<0

                    xu=xr;

            elseif fdxi*fdxr>0

                    xi=xr;

            end

    end

    if fdxi*fdxr==0

            break

    end

    plot(xr,0,'o')

    Xr(i)=xr;

end

plot(xr,0,'*','color','k');

%xr=0; xr=NaN; n=10;

for i=1:length(Xr)-1

    e(i)=(abs(Xr(i+1)-Xr(i))/Xr(i+1)*100);

end

hold off