Matlab - Ayuda por favor es urgente

%

%-----------------------------------------------------

%"x" ==> ángulo (comprendido entre -1 y +1)

%"n" ==> N° de interaccioens (entre 0 y 85, 85 porque es el valor máximo que

%el programa reconoce sin hacer infinito el denominador de la función)

%-----------------------------------------------------

%

function [y]=arcsen(x,n)

  if x>=-1 && x<=1

    for n=0:n

        j=((factorial(2.*n))./(4.^n.*(factorial(n)).^2.*(2.*n+1))).*(x.^(2.*n+1));

        d(n+1)=j;

    end

        y=sum(d);

  else y=NaN;

  end

end

%

%-----------------------------------------------------

%"x" ==> ángulo (comprendido entre -1 y +1)

%"n" ==> N° de interaccioens (entre 0 y 85, 85 porque es el valor máximo que

%el programa reconoce sin hacer infinito el denominador de la función)

%-----------------------------------------------------

%

function [y]=arccos(x,n)

  y=(pi/2)-arcsen(x,n);

end

%

%-----------------------------------------------------

%"x" ==> ángulo (comprendido entre -1 y +1)

%"n" ==> N° de interaccioens (entre 0 y 85, 85 porque es el valor máximo que

%el programa reconoce sin hacer infinito el denominador de la función)

%-----------------------------------------------------

%

function [y]=arctan(x,n)

  if x>=-1 && x<=1

     for n=0:n

         j=((-1).^n./(2.*n+1)).*(x.^(2.*n+1));

         d(n+1)=j;

     end

         y=sum(d);

  else

      y=NaN;

  end

end