Matlab - Duda con optimización en scilab

deff('fx=ec(x)','fx=2*sin(x)-(x^2/10)');

xlow=0

xhigh=4

maxit=1000

es=0.000001

opcion=1

function [xopt,fopt,iter,ea]=goldmax(xlow,xhigh,maxit,fx)

    xl=xlow

    xh=xhigh

    iter=1

    d=(xh-xl)*((sqrt(5)-1)/2*(xh-xl))

    x1=xl+d

    x2=xh+d

    f1=ec(x1)

    f2=ec(x2)

    if f2<f1 then

        xopt=x1

        fopt=f1

    else

        xopt=x2

        fopt=f2

    end

    es=0.00001

    ea=es*10

    while ea>es & iter<maxit

        iter=iter+1

        d=d*R

        if f1>f2 then

            xl=x2

            x2=x1

            x1=xl+d

            f2=f1

            f1=ec(x1)

         else

             xh=x1

             x1=x2

             x2=xh-d

             f1=f2

             f2=ec(x2)

          end

          if f1>f2 then

              xopt=x1

              fopt=f1

          else

              xopt=x2

              fopt=f2

           end

           if xopt~=0

               ea=100*(1-R)*abs((xh-xl)/xopt)

           end

      end

endfunction

[xopt,fopt,iter,ea]=gold(xlow,xhigh,maxit,ec)

disp(xopt)

function [xopt,fx,iter,ea]=metodo_seccion_aurea1(xlow,xhigh,maxiter,fun)

% Hola a todos soy novata en este foro y espero ue alguien me pueda

% ayudar.Tengo que implementar una función para encontran un máximo

% empleando el método de la sección aurea, sin embargo al ejecutarlo

% no consigo que entre en ciclo while y me imprime xopt inicial de

% ese ciclo que no es el que me interesa.Ojala alguien me ayude a

% encontrar el error .Muchas gracias os adjunto el programa que

% para los parametros dados.muchas gracias 

% [xopt,fx,iter,ea]=metodo_seccion_aurea1(xlow,xhigh,maxiter,fun) 

%[xopt,fx,iter,ea]=metodo_seccion_aurea1(0,4,1000,'2*sin(x)-(x^2/10)')

%fun='2*sin(x)-(x^2/10)';

% xlow=0;

% xhigh=4;

% maxiter=1000;

es=0.000001;

R=(sqrt(5)-1)/2;

xl=xlow;

xu=xhigh;

iter=1;

d=R*(xu-xl);

x1=xl+d;

x2=xu-d;

f=inline(fun);

f1=f(x1);

f2=f(x2);

if f1>f2

    xopt=x1;

    fx=f1;

else

    xopt=x2;

    fx=f2;

end

while 1

    d=R*d;

    if f1>f2

        xl=x2;

        x2=x1;

        x1=xl+d;

        f2=f1;

        f1=f(x1);

    else

        xu=x1;

        x1=x2;

        x2=xu-d;

        f1=f2;

        f2=f(x2);

    end

    iter=iter+1;

    if f1>f2

        xopt=x1;

        fx=f1;

    else

        xopt=x2;

        fx=f2;

    end

    if xopt~=0

        ea=(1-R)*abs((xu-xl)/xopt)*100;

    end

    if ea<=es || iter>=maxiter

        break;

    end

end

EJECUCION

>> [xopt,fx,iter,ea]=metodo_seccion_aurea1(0,4,1000,'2*sin(x)-(x^2/10)')

xopt =

    1.4276

fx =

    1.7757

iter =

    40

ea =

  7.5679e-007