Programa biseccion matlab

function biseccion1

f1=input('ingrese la funcion ','s');

f=inline(f1);

a=input('ingrese el limite inferior a: ');

b=input('ingrese el limite superior b: ');

n=input('ingrese el numero de iteraciones n:  ');

f=inline(f1)

for i=1:1:n

    xr=(a+b)/2;

    if f(a)*f(xr)<0;

        b=xr

    else f(xr)*f(b)<0;

        a=xr

    end

end

disp('La aproximacion de la raiz es ')

disp(xr)

function biseccion2

f1=input('ingrese la funcion ','s');

a=input('ingrese el limite inferior a: ');

b=input('ingrese el limite superior b: ');

es=input('ingrese el error tolerable :  ');

f=inline(f1)

cont=1

ea=abs((b-a)/2)*100

while ea>es

    cont=cont+1;

    xr=(a+b)/2;

    if f(a)*f(xr)<0;

        b=xr

        ea=abs((xr-a)/xr)*100

    else f(xr)*f(b)<0;

        a=xr

      ea=abs((xr-b)/xr)*100

    end

end

disp('numero de iteraciones ')

disp(cont)

disp('La aproximacion de la raiz es ')

disp(xr)

function gauss

A=input('ingrese la matriz A= ');

B=input('ingrese la columna B= ');

n=length(A);

AB=[A,B];

X=zeros(n,1);

for j=1:1:n;

    if AB(j,j)==0

        disp('reordenar la matriz');

    end

    AB(j,:)=(1/AB(j,j))*AB(j,:);

    for i=[1:j-1,j+1:n];

        m=-AB(i,j);

        AB(i,:)=m*AB(j,:)+AB(i,:);

    end

end

X=AB(:,n+1)

disp('la solucion es')

disp(X)

function raphsonmultiple

f=input('ingrese la funcion f(x)= ','s');

g=input('ingrese la funcion g(x)= ','s');

x1=input('ingrese el valor inicial x1= ');

y1=input('ingrese el valor inicial y1= ');

n=input('ingrese el numero de iteraciones n= ');

syms x y

dfx=diff(f,x);

dfy=diff(f,y);

dgx=diff(g,x);

dgy=diff(g,y);

for i=1:1:n

    dfx1=subs(dfx,{x,y},{x1,y1});

...